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Preface 


These lecture notes provide an introduction to some of the main concepts of topo¬ 
logical insulators, a branch of solid state physics that is developing at a fast pace. 
They are based on a one-semester course for MSc and PhD students at the Eotvos 
University, Budapest, which the authors have been giving since 2012. 

Our aim is to provide an understanding of the core topics of topological insulators 
- edge states, bulk topological invariants, bulk-boundary correspondence - with as 
simple mathematical tools as possible. We restricted our attention to one- and two- 
dimensional band insulators. We use noninteracting lattice models of topological 
insulators, and build these up gradually to arrive from the simplest one-dimensional 
case (the Su-Schrieffer-Heeger model for polyacetylene) to two-dimensional time- 
reversal invariant topological insulators (the Bernevig-Hughes-Zhang model for 
HgTe). In each case we introduce the model first, discuss its properties, and then 
generalize. The prerequisite for the reader is quantum mechanics and not much else: 
solid state physics background is provided as we go along. 

Since this is an introduction, rather than a broad overview, we try to be self- 
contained and give citations to the current literature only where it is absolutely nec¬ 
essary. For a broad overview, including pointers to the original papers and current 
topics, we refer the reader to review articles and books in the Introduction. 

Supporting material for these lecture notes in the form of ipython notebooks will 
be made available online. 

Despite our efforts, the book inevitably contains typos, errors, and less compre¬ 
hensible explanations. We would appreciate if you could inform us of any of those; 
please send your comments to janos . asboth@wigner . mta . hu. 
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Introduction 


The band theory of electric conduction was one of the early victories of quantum 
mechanics in the 1920s. It gave a simple explanation of how some crystalline mate¬ 
rials are electric insulators, even though electrons in them can hop from one atom to 
the next. In the bulk of a band insulator, the electrons occupy eigenstates that form 
energy bands. In a band insulator, there are no partially filled bands: completely 
filled bands are separated by an energy gap from completely empty bands, the gap 
represents the energy cost of mobilizing electrons. In contrast, materials with par¬ 
tially filled bands are conductors, where there are plane wave states available to 
transmit electrons across the bulk at arbitrarily low energy. Although we now know 
of situations where band theory is inadequate (e.g., for Mott insulators), it remains 
one of the cornerstones of solid states physics. 

The discovery of the Quantum Hall Effect (1980) has shown that the simple 
division into band insulators and metals is not the end of the story, not even in band 
theory. In the quantum Hall effect, a strong magnetic field confines the motion of 
electrons in the bulk, but the same field forces them into delocalized edge states on 
the surface. A two-dimensional metal in strong magnetic field is thus an insulator 
in the bulk, but conducts along the surface, via a discrete number of completely 
open edge state channels (in the language of the Landauer-Biittiker formalism). 
The number of edge state channels was linked to the Chern number, a topological 
invariant of the occupied bands (1982). 

Over the last twenty years, theoretical progress over artificial systems has shown 
that the external magnetic field is not necessary for an insulator to have robust con¬ 
ducting edge states: instead, the nontrivial topology of the occupied bands is the 
crucial ingredient. The name topological insulator was coined for such systems, 
and their study became a blossoming branch of solid state physics. Following the 
theoretical prediction (Bernevig, Hughes and Zhang, 2006 ID), electronic transport 
measurements confirmed that a thin layer of HgTe is a topological insulator (Konig 
et al, 2007 fllU l. Since that time, a host of materials have been shown to be three- 
dimensional topological insulators, and thin films and quantum wires shown to be 
two- and one-dimensional topological insulators 0. 
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The intense theoretical interest in topological insulators has led to signature re¬ 
sults, such as the “the periodic table of topological insulators” 1^ . which shows 
that similarly to phase transitions in statistical mechanics, it is the dimensionality 
and the basic symmetries of an insulator that decide whether it can be a topological 
insulator or not. Although it was derived by different ways of connecting topological 
insulators of various dimensions and symmetries (so-called dimensional reduction 
schemes), the mathematically rigorous proof of the periodic table is still missing. 

The field of topological insulators is very active, with many experimental chal¬ 
lenges and open theoretical problems, regarding the effects of electron-electron in¬ 
teraction, extra crystalline symmetries, coupling to the environment, etc. 


Literature 

To get a quick and broad overview of topological insulators, with citations for rel¬ 
evant research papers, we recommend the review papers nil ESI 13. For a more 
in-depth look, there are already a few textbooks on the subject (by Bemevig and 
Hughes in, by Shen f^ , and one edited by Franz and Molenkamp ifTOl ). To see 
the link between momentum-space topology and physics in a broader context, we 
direct the reader to a book by Volovik 041 . 

There are also introductory courses on topological insulators with a broad scope. 
We recommend the lectures by Charles Kane (the video recording of the version 
given at Veldhoven is freely available online), and the online EdX course on topol¬ 
ogy in condensed matter by a group of lecturers, with the corresponding material 
collected at topocondmat.org. 


These lecture notes 

Our aim with this set of lecture notes is to complement the literature cited above: 
we wish to provide a close look at some of the core concepts of topological insula¬ 
tors with as simple mathematical tools as possible. Using one-and two-dimensional 
noninteracting lattice models, we explain what edge states and what bulk topolog¬ 
ical invariants are, how the two are linked (this is known as the bulk-boundary 
correspondence), the meaning and impact of some of the fundamental symmetries. 

To keep things as simple as possible, throughout the course we use noninteract¬ 
ing models for solid state systems. These are described using single-particle lattice 
Hamiltonians, with the zero of the energy corresponding to the Fermi energy. We 
use natural units, with h=\ and length measured by the lattice constant. 
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Chapter 1 

The Su-Schrieffer-Heeger (SSH) model 


We take a hands-on approach and get to know the basic concepts of topological 
insulators via a concrete system: the Su-Schrieffer-Heeger (SSH) model describes 
spinless fermions hopping on a one-dimensional lattice with staggered hopping am¬ 
plitudes. Using the SSH model, we introduce the concepts of single-particle Hamil¬ 
tonian, the difference between bulk and boundary, chiral symmetry, adiabatic equiv¬ 
alence, topological invariants, and bulk-boundary correspondence. 



Fig. 1.1 Geometry of the SSH model. Filled (empty) circles are sites on sublattice A (B), each 
hosting a single state. They are grouped into unit cells: the n = 6th cell is circled by a dotted 
line. Hopping amplitudes are staggered: intracell hopping v (thin lines) is different from intercell 
hopping w (thick lines). The left and right edge regions are indicated by blue and red shaded 
background. 


1.1 The SSH Hamiltonian 

The Su-Schrieffer-Heeger (SSH) model describes electrons hopping on a chain 
(one-dimensional lattice), with staggered hopping amplitudes, as shown in Fig. [U] 
The chain consist of N unit cells, each unit cell hosting two sites, one on sublattice 
A, and one on sublattice B. Interactions between the electrons are neglected, and so 
the dynamics of each electron is described by a single-particle Hamiltonian, of the 
form 


N N-l 

H=v^ (^\m,B) {m,A\+h.c.) +w (^\m+l,A) {m,B\+h.c.). (1-1) 

m=l m=l 


1 
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Here \m,A) and \m,B), with m G {1,2,... ,A^}, denote the state of the chain where 
the electron is on unit cell m, in the site on sublattice A, respectively, B, and h.c. 
stands for Hermitian Conjugate (e.g., h.c. of \m,B) {m,A\ is \m,A) {m,B\). 

The spin degree of freedom is completely absent from the SSH model, since no 
term in the Hamiltonian acts on spin. Thus, the SSH model describes spin-polarized 
electrons, and when applying the model to a real physical system, e.g., polyacety¬ 
lene, we have to always take two copies of it. In this chapter we will just consider a 
single copy, and call the particles fermions, or electrons, or just particles. 

We are interested in the dynamics of fermions in and around the ground state of 
the SSH model at zero temperature and zero chemical potential, where all negative 
energy eigenstates of the Hamiltonian are singly occupied (because of the Pauli 
principle). As we will show later, due to the absence of onsite potential terms, there 
are N such occupied states. This situation - called half hlling - is characteristic of 
the simplest insulators such as polyacetylene, where each carbon atom brings one 
conduction electron, and so we hnd 1 particle (of each spin type) per unit cell. 

For simplicity, we take the hopping amplitudes to be real and nonnegative, v,w> 
0. If this was not the case, if they carried phases, v = | v| and w = |w| , with 

0H' S [0, n), these phases could always be gauged away. This is done by a redefini¬ 
tion of the basis states: \m,A) —>• ^ ^nd \m,B) —)■ 

The matrix for the Hamiltonian of the SSH model, Eq. ([U}, on a real-space 
basis, for a chain of A = 4 unit cells, reads 


/O V 0 0 0 0 0 0\ 
vOwOOOOO 
OwOvOOOO 
OOvOwOOO 
OOOwOvOO 
0000 vOwO 
OOOOOwOv 

\0 0 0 0 0 0 V 0/ 


( 1 . 2 ) 


External and internal degrees of freedom 

There is a practical representation of this Hamiltonian, which emphasizes the sep¬ 
aration of the external degrees of freedom (unit cell index m) from the internal de¬ 
grees of freedom (sublattice index). We can use a tensor product basis. 


1^, Oj) y 1 ^ 7)0 |Oj) € .^^^xtemal ^'^^nternal, (1-2) 


with m= 1,... ,N, and a G {A,B}. On this basis, with the Pauli matrices. 


do = 




(1.4) 


the Hamiltonian can be written 




1.2 Bulk Hamiltonian 
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H = v 


m=1 m= 1 \ 




(1.5) 


The intracell hopping shows up here as an intracell operator, while the intercell 
hopping as a hopping operator. 


1.2 Bulk Hamiltonian 

As every solid-state system, the long chain of the SSH model has a bulk and a 
boundary. The bulk is the long central part of the chain, the boundaries are the two 
ends, or “edges” of the chain, indicated by shading in Fig. o We first concentrate 
on the bulk, since, in the thermodynamic limit of N it is much larger than 
the boundaries, and it will determine the most important physical properties of the 
model. Although the treatment of the bulk using the Fourier transformation might 
seem like a routine step, we detail it here because different conventions are used in 
the literature. More on this in Appendix ??. 

The physics in the bulk, the long central part of the system, should not depend on 
how the edges are defined, and so for simplicity we set periodic (Born-von Karman) 
boundary conditions. This corresponds to closing the bulk part of the chain into a 
ring, with the bulk Hamiltonian defined as 


N 


^buik = 52 {v\in,B) {m,A\+w\{mmodN) + 1,A) {m,B\'j + h.c.. (1.6) 


m=\ 


We are looking for eigenstates of this Hamiltonian, 


HMMk))=E„ik)\%ik)), 


(1.7) 


with n G {1,... ,2N}. 

Bulk momentum-space Hamiltonian 

Due to the translation invariance of the bulk, Bloch’s theorem applies, and we look 
for the eigenstates in a plane wave form. We introduce the plane wave basis states 
only for the external degree of freedom. 



where the wavenumber k was chosen to take on values from the first Brillouin zone. 
The Bloch eigenstates read 
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\W„{k)) = \k)®\un{k))- \un{k))=an{k)\A)+bn{k)\B). (1.9) 

The vectors \un{k)) € ^6 eigenstates of the bulk momentum-space Hamil¬ 

tonian H{k) defined as 

^(^) = (^|4uik|^)= E (^,«|//buik|^,j3)-|a)(j3|; (1.10) 

H{k)\Un{k))=En{k)\un{k)). (1.11) 


Periodicity in wavenumber 


Although Eq. has a lot to do with the continous-variable Bloch theorem, 
^n,k(x) = jt(-v), this correspondence is not direct. In a discretization of the 

continuous-variable Bloch theorem, the internal degree of freedom would play the 
role of the coordinate within the unit cell, which is also transformed by the Fourier 
transform. Thus, the function cell-periodic, m„ 1) = u„ i^{x), but not pe¬ 

riodic in the Brillouin zone, M„j(-+ 2 re(^+ 1) u„ i^{x). Our Fourier transform acts 
only on the external degree of freedom, and as a result, we have periodicity in the 
Brillouin zone, 

H{k + 27t)=H{k)- \u„{k + 27t)) = \un{k)) . (1.12) 


This convention simplifies the formulas for the topological invariants immensely. 
Note, however, that the other convention, the discretization of the Bloch theorem, is 
also widely used in the literature. We compare the two approaches in Appendix ??. 

As an example, for the SSH model on a chain of A = 4 unit cells, the Schrodinger 
equation, Eq. using Eq. ([L^, translates to a matrix eigenvalue equation. 


^0v00000 w^ 
vOwOOOOO 
OwOvOOOO 
OOvOwOOO 
OOOwOvOO 
OOOOvOwO 
OOOOOwOv 
yw OOOOOvOy 


/ a{k)e‘^ \ 


/ a(k)e'^ \ 

b{k)e^^ 


b{k)e^ 



a{k)e^‘^ 

b{k)e^''^ 

a\k)e^'^ 

=m 

b\k)e^''^ 

a{k)e^^^ 

blky^’^ 


b{k)e^^^ 

a(k)e^‘’^ 



\b{k)e^‘y 


\b{li)e^‘’^) 


(1.13) 


The Schrodinger equation defining the matrix H{k) of the bulk momentum-space 
Hamiltonian reads 


/ 0 

\v-\-we’^ 0 J ’ 



H{k) = 


(1.14) 








1.2 Bulk Hamiltonian 
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1.2.1 The hopping is staggered to open a gap 


The dispersion relation of the bulk can be read off from Eq. ( |1.14| l, using the fact 
that H{kY =E{kYl 2 . This gives us 


m = 


v + e 


= + 2vw cos k. 


(1.15) 


We show this dispersion relation for five choices of the parameters in Fig. |1.2 



Fig. 1.2 Dispersion relations of the SSH model, Eq. jl.151 , for five settings of the hopping ampli¬ 
tudes: (a): V — 1,w = 0; (b): v = l,w — 0.6; (c): v = w = 1; (d): v = 0.6,w = 1; (e): v = 0,w = 1. 
In each case, the path of the endpoints of the vector d(k) representing the bulk momentum-space 
Hamiltonian, Eqs. jl.l7| l and l |1.181 , are also shown on the dx,dy plane, as the wavenumber is 
sweeped across the Brillouin zone, k = 0 —> 271. . 


As long as the hopping amplitudes are staggered, v 7 ^ w, (Figs. 1.2 (a),(b),(d),(e), 
there is an energy gap of 2A separating the lower, filled band, from the upper, empty 
band, with 


A = min*:£'(k) = |v —w|. (1-16) 

Without the staggering, i.e., if v = w, (Fig. |1.2| (c), the SSH model describes a con¬ 
ductor. In that case there are plane wave eigenstates of the bulk available with arbi¬ 
trarily small energy, which can transport electrons from one end of the chain to the 
other. 

The staggering of the hopping amplitudes occurs naturally in many solid state 
systems, e.g., polyacetylene, by what is known as the Peierls instability. A detailed 
analysis of this process neccesitates a model where the positions of the atoms are 
also dynamical ly . Nevertheless, we can understand this process intuitively just 
from the effects of a slight staggering on the dispersion relation. As the gap due 
to the staggering of the hopping amplitudes opens, the energy of occupied states is 
lowered, while unoccupied states move to higher energies. Thus, the staggering is 
energetically favourable. 
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1.2.2 Information beyond the dispersion relation 


Although the dispersion relation is useful to read off a number of physical properties 
of the bulk of the system (e.g., group velocities), there is also important information 
about the bulk that it does not reveal. Stationary states do not only have an energy 
and wavenumber eigenvalue, but also an internal structure, represented by the com¬ 
ponents of the corresponding vector \u„{k)) G internal- We now define a compact 
representation of this information for the SSH model. 

The bulk momentum-space Hamiltonian tt(k) of any two-band model (i.e., a 
model with 2 internal states per unit cell), reads 

N(/c) = dx(k)dx + dy{k)dy + d^(k)d^ = t/o(^)^o+ d(k)a. (1-17) 

For the SSH model, d{){k) = 0, and the real numbers dxj^z G the components of 
the k-dependent 3-dimensional vector d(k), read 

dx{k) =v+ wco^,k\ dy{k) = w s\Vik\ dz{k)—Q. (1-18) 


The internal structure of the eigenstates with momentum k is given by the direction 
in which the vector d{k) of Eq. ( |1.18[ l points (the energy is given by the magnitude 
of d(k); for details see Sect. 2.5 1 . 

As the wavenumber runs through the Brillouin zone, k = 0 —> 27Z, the path that 
the endpoint of the vector d(k) traces out is a closed circle of radius w on the dx,dy 
plane, centered at (v,0). For more general 2-band insulators, this path need not be a 
circle, but it needs to be a closed loop due to the periodicity of the bulk momentum- 
space Hamiltonian, Eq. p.l2| l, and it needs to avoid the origin, to describe an insu¬ 
lator. The topology of this loop can be characterized by an integer, the bulk winding 
number V. This counts the number of times the loop winds aroung the origin of the 
dx,dy plane. For example, in Fig. |1.2[ ^f),(g), we have v = 0, in Fig. 1.2 i),(j), we 
have V = 1, while in Fig. 1 1.2(h), the winding number v is undefined. 


1.3 Edge states 

Like any material, the SSH Hamiltonian does not only have a bulk part, but also 
boundaries (which we refer to as ends or edges). The distinction between bulk and 
edge is not sharply defined, it describes the behaviour of energy eigenstates in the 
thermodynamic limit. In the case we consider in these lecture notes, the bulk is trans¬ 
lation invariant, and then the we can distinguish edge states and bulk states by their 
localized/delocalized behaviour in the thermodynamic limit. We will begin with the 
fully dimerized limits, where the edge regions can be unambiguously defined. We 
then move away from these limits, and use a practical definition of edge states. 










1.3 Edge states 
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1.3.1 Fully dimerized limits 


The SSH model becomes particularly simple in the two fully dimerized cases: if the 
intercell hopping amplitude vanishes and the intracell hopping is set to 1, v = 1, w = 
0, or vice versa, v = 0,w = 1. In both cases the SSH chain falls apart to a sequence 


of disconnected dimers, as shown in Fig. 1.3 



Fig. 1.3 Fully dimerized limits of the SSH model, where the chain has fallen apart to disconnected 
dimers. In the trivial case (top, only intracell hopping, v = 1, w = 0), every energy eigenstate is an 
even or an odd superposition of two sites at the same unit cell. In the topological case, (bottom, 
only intercell hopping, v = 0,^ = 1), dimers are between neighboring unit cells, and there is 1 
isolated site per edge, that must contain one zero-energy eigenstate each, as there are no onsite 
potentials. . 


The bulk in the fully dimerized limits has flat bands 

In the fully dimerized limit, one can choose a set of energy eigenstates which are 
restricted to one dimer each. These consist of the even (energy £ = -|-1) and odd 
(energy £ = — 1) superpositions of the two sites forming a dimer. 

In the V = 1, w = 0 case, which we call trivial, we have 

v=l,w = 0: H{\m,A) zL\m,B)) = zL{\m,A) zt\m,B)). (1-19) 

The bulk momentum-space Hamiltonian is H(k) = dx, independent of the wavenum¬ 
ber k. 

In the V = 0, w = 1 case, which we call topological, each dimer is shared between 
two neighboring unit cells, 

v = 0,w=l: i/(|OT,£) ± |m-|-1,A)) = ±(|m,£) ± |m-|-1,A)), (1-20) 

for m = 1,...,A^— 1. The bulk momentum-space Hamiltonian now is H{k) = 
dx cos k + Gy sin k. 

In both fully dimerized limits, the energy eigenvalues are independent of the 
wavenumber, £(k) = 1. In this so-called flat-band limit, the group velocity is zero, 
which again shows that as the chain falls apart to dimers, a particle input into the 
bulk will not spread along the chain. 
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The edges in the fully dimerized limit can host zero energy states 

In the trivial case, v = l,w = 0, all energy eigenstates of the SSH chain are given 
by the formulas of the bulk, Eq. ( |1.19| ). A topological, fully dimerized SSH chain, 
with V = 0, w = 1, however, has more energy eigenstates than those listed Eq. ( |1.20[ ). 
Each end of the chain hosts a single eigenstate at zero energy, 

y = 0,w=l; H\l,A)=H\N,B)=0. (1.21) 

These eigenstates have support on one site only. Their energy is zero because onsite 
potentials are not allowed in the SSH model. These are the simplest examples of 
edge states. 


1.3.2 Moving away from the fully dimerized limit 

We now examine what happens to the edge states as we move away from the fully 
dimerized limit. To be specific, we examine how the spectrum of an open topological 
chain, v = 0,w = 1, of A = 10 unit cells changes, as we continuously turn on the 
intracell hopping amplitude v. The spectra, Eig. |1.4[ reveal that the energies of the 
edge states remain very close to zero energy. 

The wavefunctions of almost-zero-energy edge states have to be exponentially 
localized at the left/right edge, because the zero of energy is in the bulk band gap. A 
plot of the wavefunctions (which have only real components, since the Hamiltonian 
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Fig. 1.4 Energy spectrum and wave functions of a finite-sized SSH model. The number of unit 
cells is N = 10. (a) Energy spectrum of the system for intercell hopping amplitude vr = 1 as a 
function the intracell hopping amplitude v. v < 1 (v > 1) corresponds to the topological (trivial) 
phases, (b) and (c) shows the wave functions of the hybridized edge states, while (d) shows a 
generic bulk wave function. 
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is real), Fig. |1.4| reveals that the almost-zero-energy eigenstates are odd and even 
superpositions of states localized exponentially on the left and right edge. This is a 
result of the exponentially small overlap between the left and the right edge states. 
We will later show, in Sect. 1.5.2 that the edge-state energies are also controlled 
by this overlap, and are of the order E = , with a localization length ^ = 

l/log(v/w). 

There is an important property of the right/left edge states, which is only re¬ 
vealed by the plot of the wavefunctions. Fig. 1.4 These states have nonvanishing 
components only on the A/B sublattice. 

In the following, we show the generality of these properties, and show the link 
between the bulk winding number and the presence/absence of edge states, known 
as bulk-boundary correspondence. In the case of the SSH model, all this hinges on 
a property of the model known as chiral symmetry. 


1.4 Chiral symmetry 

In quantum mechanics, we say that a Flamiltonian H has a symmetry represented 
by a unitary operator U if 


UHU''=H. ( 1 . 22 ) 

In case of a symmetry, U and H can be diagonalized together, and therefore, H has 
no matrix elements between two eigenstates of U with different eigenvalues. This 
can be understood as a superselection rule: if we partition the Hilbert space into 
different sectors, i.e., eigenspaces of U, labeled by the corresponding eigenvalues, 
then the dynamics as defined by H can be regarded separately in each sector. 


No unitary symmetries 


A unitary symmetry can be simply made to disappear if we restrict ourselves to 
one sector of the Hilbert space. This is how we obtained the bulk momentum-space 
Hamiltonian, in Sect. 11.21 where the symmetry was the lattice translation operator 

l,A)i 


U = \m + 1,A) {m,A\ + \m 

were the quasimomenta k. 


1,B) {m.B\, and the labels of the superselection sectors 


A different type of symmetry 

The word “symmetry” is also used in a different sense in condensed matter physics. 
We say that a system with Hamiltonian H has chiral symmetry, if 

= -H, 


(1.23) 
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with an operator F that is not only unitary, but fulfils some other criteria as well. 
Notice the extra minus sign on the right hand side. This has important consequences, 
which we come to later, but first discuss the criteria on the symmetry operator. 

First, the chiral symmetry operator has to be unitary and Hermitian, = F, 
which can be written succintly as 

r^r = f^ = l. (1.24) 

The reason for this requirement is that if the operator F^ was nontrivial, it would 
represent a unitary symmetry, since 

FFHff = -tut =H. (1.25) 

This could still leave room for the chiral symmetry operator to square to a state- 
independent phase, F^ = e“^. However, this can be got rid of by a redefinition of the 
chiral symmetry, F —> . 

Second, it is also required that the sublattice operator F be local. The system 
is assumed to consist of unit cells, and matrix elements of F between sites from 
different unit cells should vanish. In the SSH chain, this means that for m ^ m', we 
have {m,a\r \m',a') = 0, for any a, a' G {A,B). To keep things simple, we can 
demand that the sublattice operator act in the same way in each unit cell (although 
this is not strictly necessary), its action represented by a unitary operator 7 acting 
on the internal Hilbert space of one unit cell, i.e., 

N 

F = 70 70 ... 07 = 07 , (1.26) 

m=l 

where N is the number of unit cells. 

A third requirement, which is often not explicitly stated, is that the chiral sym¬ 
metry has to be robust. To understand what we mean by that, first note that in solid 
state physics, we often deal with Hamiltonians with many local parameters that vary 
in a controlled or uncontrolled way. An example is the SSH model, where the values 
of the hopping amplitudes could be subject to spatial disorder. We gather all such 
parameters in a formal vector, and call it G S. Here S is the set of all realizations 
of disorder that we investigate. Instead of talking about the symmetries of a Hamil¬ 
tonian H, we should rather refer to symmetries of a set of Hamiltonians {//((*)}, for 
all ^ G S. This set has chiral symmetry represented by F if 

G S ; FH(4)F = -H, (1.27) 

with the symmetry operator F independent of the parameters ^. This is the robust¬ 
ness of the chiral symmetry. 



1.4 Chiral symmetry 
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1.4.1 Consequences of chiral symmetry for energy eigenstates 

We now come to the consequences of chiral symmetry, which are very different 
from those of conventional symmetries, due to the extra minus sign in its dehnition, 
Eq. ([T^. 


Sublattice symmetry 

Chiral symmetry is also called sublattice symmetry. Given the chiral symmetry op¬ 
erator r, we can define orthogonal sublattice projectors and Pg, as 

PA^\{l + t)- (1.28) 

where I represents the unity operator on the Hilbert space of the system. Note that 
Pa+Pb = I, and PaPb = 0. The dehning relation of sublattice symmetry, Eq. ( |1.23[ ), 
can be written in an equivalent form by requiring that the Hamiltonian induces no 
transitions from any site on one sublattice to any site on the same sublattice, 

PaHPa = PbHPb = 0; H = PaHPb + PbHPa ■ (1.29) 

In fact, using the projectors Pa and Pg is an alternative and equivalent way of dehn¬ 
ing chiral symmetry. 


Symmetric spectrum 

The spectrum of a chiral symmetric Hamiltonian is symmetric. Eor any state with 
energy E, there is a chiral symmetric partner with energy —E. This is simply seen, 

H|v/„) =£„|r„) Hr\\i/„) = -tH\\iin) = -tE„\\i/„) = -E„r\\i/„). 

(1.30) 

This carries different implications for nonzero energy eigenstates and zero energy 
eigenstates. 

Eor E„ 7 ^ 0, the states | y/,,) and E | y/,,) are eigenstates with different energy, and, 
therefore, have to be orthogonal. This implies that every nonzero energy eigenstate 
of H has equal support on both sublattices, 

If£„7^0; 0 = {Wn\r Wn) = {Wn\PAWn) - {Wn\PBWn) ■ (1-31) 

Eor E„ = 0, zero energy eigenstates can be chosen to have support on only one 
sublattice. This is because 


Ifi/|V4) = 0: HpAlB\Wn)=H{\Wn)±P'\Wn)) =0. 


(1.32) 
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These projected zero-energy eigenstates are eigenstates of f, and therefore are chi¬ 
ral symmetric partners of themselves. 


1.4.2 Sublattice projectors and chiral symmetry of the SSH model 

The Hamiltonian of the SSH model, Eq. ([U), is bipartite: the Hamiltonian includes 
no transitions between sites with the same sublattice index. The projectors to the 
sublattices read 


N N 

PA=Y.\m,A){n,A\-, Pg = \m,B) {n,B\. (1.33) 

m=l m=l 

Chiral symmetry is represented by the sublattice operator E^, that multiplies all com¬ 
ponents of a wavefunction on sublattice B by (-1), 

L=Pa-Pb. (1.34) 

Note that this operator has the properties required of the chiral symmetry operator 
above: it is unitary, Hermitian, and local. 

The chiral symmetry of the SSH model is a restatement of the fact that the Hamil¬ 
tonian is bipartite, 

PaPiPa = PbHPb = 0; = -H- (1-35) 

This relation holds because H only contains terms that are multiples of \m,A) {in',B\, 
or of \in,B) {in',A\ with in,in' G Z. Upon multiplication from the left and the right 
by E^-, such a term picks up a single factor of —1 (because of the multiplication 
from the left or because of the multiplication from the right). Note that this relation, 
equivalent to an anticommutation of N and E^, holds whether or not the hopping 
amplitudes depend on position: therefore, the chiral symmetry represented by E^ 
has the required property of robustness. 


1.4.3 Consequence of chiral symmetry: Bulk winding number for 
the SSH model 

The path of the endpoint of d{k), as the wavenumber goes through the Brillouin 
zone, k = 0 ^ 2n, is a closed path on the dx,dy plane. This path has to avoid the 
origin, d = 0: if there was a k at which d(k) = 0, the gap would close at this k, and 
we would not be talking about an insulator. Because of chiral symmetry, the vector 
d(k) is restricted to lie on the d^dy plane. 
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a^H{k)az = 0 => d^ = 0. (1.36) 

This is then a closed, directed loop on the plane, and thus has a well defined integer 
winding number about the origin. 


Winding number as the multiplicity of solutions 


The simplest way to obtain the winding number graphically is counting the num¬ 
ber of times d(A:) intersects a curve that goes from the origin of the dx,dy plane to 
infinity. 

1. Since d{k) is a directed curve, it has a left side and a right side. Paint the left side 
blue, the right side red, as shown in Fig. |1.5| (a). 

2. Take a directed curve ^ going from 0 to infinity. We can call this the “line of 
sight to infinity”, although it need not be a straight line. A simple choice is the 


half-infinite line, dy — Q,dx> 0. Two other choices are shown in Fig. 1.5 

3. Identify the intersections of d(A:) with .if. 

4. Each intersection has a signature: this is -fl if the line of sight meets it from the 
blue side, — 1 for the red side. 

5. The winding number v is the sum of the signatures. 

To show that the winding number v defined above is a topological invariant, we 
need to consider how it can change under continuous deformations of .if or of d(A:). 
Due to the deformations the intersections of .if and d(A:) can move, but this does 
not change v. They can also appear or disappear, at points where .if and d(k) touch. 
However, they can only appear or disappear pairwise: a red and a blue intersection 
together, which does not change V. As an example, the two choices of the line of 


sight if in Fig. 1.5 (a), have 1 or 3 intersections, but the winding number is -|-1, for 
either of them. 


Winding number as an integral 

The winding number can also be written as a compact formula using the unit vector 
d, defined as 


Fig. 1.5 The endpoints of the 
vector d(^) as k goes across 
the Brillouin zone (red or blue 
closed circles). 


A dy 
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d = 


d 

Wl 


(1.37) 


Thi s is the result of projecting the curve of d(A:) to the unit circle, as shown in 
Fig. 1.5 (b). The vector d{k) is well defined for all k because d{k) ^ 0. 

You can check easily that the winding number v is given by 


V = (^d(k) X dk. (1.38) 

To calculate v directly from the bulk momentum-space Hamiltonian, note that it 
is off-diagonal (in the basis of eigenstates of the chiral symmetry operator a^), 

m) = ; m = d,{k) - idy{k). (1.39) 

The winding number of d(k) can be written as an integral, using the complex loga¬ 
rithm function, log(|/i| = log \h\ -f iaigh. It is easy to check that 

V=;^/ dk^\ogh{k). (1.40) 

zTTi J —71 dk 


Here during the calculation of the integral, the branch cut for the logarithm is always 
shifted so that the derivative is always well defined. The above integral is always 
real, since \h{k = —7r)| = \h{k = 7t)\. 


Winding number of the SSH model 

For the SSH model, the winding number is either 0 or 1, depending on the parame¬ 
ters. In the trivial case, when the intracell hopping dominates the intercell hopping, 

V > w, the winding number is V = 0. In the topological case, when w > v, we have 

V = 1. 

To change the winding number v of the SSH model, we need to either a) pull the 
path of d(k) through the origin in the dx,dy plane, or (b) lift it out of the plane and 
put it back on the plane at a different position. This is illustrated in Fig. |1.6| Method 
(a) requires closing the bulk gap. Method (b) requires breaking chiral symmetry. 


1.5 Number of edge states as topological invariant 

We now introduce the notion of adiabatic deformation of insulating Hamiltonians. 
An insulating Hamiltonian is adiabatically deformed if 

• its parameters are changed continuously, 

• the important symmetries of the system are maintained. 



1.5 Number of edge states as topological invariant 
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• the bulk gap around E = 0 remains open. 

The deformation is a fictitious process, and does not take place in time. However, 
if we do think of it as a process in real time, the adiabatic theorem ns tells us, 
that, starting from the many-body ground state (separated from excited states by 
the energy gap), and performing the deformation slowly enough, we end up in the 
ground state, at least as far as the bulk of the system is concerned. At the edges of 
a system, changes can occur, and there is a subtle point to be made about adiabatic 
deformations being slow, but not too slow, that the edges should still be considered 
separately. We will come back to this point in Chapt. 


Adiabatic equivalence of Hamiltonians 

Two insulating Hamiltonians are said to be adiabatically equivalent or adiabatically 
connected if there is an adiabatic deformation connecting them, that respects the 
important symmetries. For example, in the phase diagram Fig. |1.7| of the SSH model, 
the two Hamiltonians corresponding to the two black points in the topological phase 
(w > v) are adiabatically connected, as one can draw a path between them which 
does not cross the gapless topological-trival phase boundary w = v. 


Topological invariant 

We call an integer number characterizing an insulating Hamiltonian a topological 
invariant, or adiabatic invariant, if it cannot change under adiabatic deformations. 
Note that the use of adiabatic deformations implies two properties of the topological 



Fig. 1.6 The endpoints of the vector d(^) as k goes across the Brillouin zone (red or blue closed 
circles), for various parameter settings in the SSH model. In (a), intracell hopping v is set to 0.5, 
and intercell hopping is gradually increased from 0 to 1 (outermost red circle). In the process, 
the bulk gap was closed and reopened, as the origin (black point) is included in some of the blue 
circles. The winding number changed from 0 to 1. In (b), we gradually decrease intracell hopping v 
to 0 so that the winding number again changes to 1. In the process, the bulk gap closes and reopens. 
In (c), keeping w = I, we increase the intracell hopping from 0.5 to 2.5, but avoid closing the bulk 
gap by introducing a sublattice potential, //jubiattice = We do this by tuning a parameter 0 from 
0 to 7t, and setting v = 1.5 — cos0, and u = sin0. At the end of the process, 0 = 0, there is no 
sublattice potential, so chiral symmetry is restored, but the winding number is 0 again. 
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invariant; 1) it is only well defined in the thermodynamic limit, 2) it depends on the 
symmetries that need to be respected. An example for a topological invariant is the 
winding number v of the SSH model. 

We know that two insulating Hamiltonians are not adiabatically equivalent if their 
topological invariants differ. Consider as an example two Hamiltonians correspond¬ 
ing to two points on different sides of the phase boundary in Fig. |1.7| of the SSH 
model. One might think that although there is no continous path connecting them in 
the phase diagram, continuously modifying the bulk Hamiltonian by the addition of 
extra terms can lead to a connection between them. However, their winding numbers 
differ, and since winding numbers cannot change under adiabatic deformation, we 
know that they are not adiabatically equivalent. 


Number of edge states as a topological invariant 


We have seen in Sect. 1.3.2 that the number of edge states at one end of the SSH 
model was an integer that did not change under a specific type of adiabatic defor¬ 
mation. We now generalize this example. 

Consider energy eigenstates at the left end of a gapped chiral symmetric one¬ 
dimensional Hamiltonian in the thermodynamic limit, i.e., with length N ^ in 
an energy window from —e < E < e, with e in the bulk gap. There can be nonzero 
energy edge states in this energy window, and zero energy edge states as well. Each 
nonzero energy state has to have a chiral symmetric partner, with the state and its 
partner occupying the same unit cells (the chiral symmetry operator is a local uni¬ 
tary). The number of zero energy states is finite (because of the gap in the bulk), and 
they can be restricted to a single sublattice each. There are Na zero energy states on 
sublattice A, and Nb states on sublattice B. 


Fig. 1.7 Phase diagram of 
the SSH model. The wind¬ 
ing number of the bulk 
momentum-space Hamil¬ 
tonian H{k) can be V = 0, 
if V > w, or V = 1, if V < w. 
This defines the trivial (gray) 
and the topological phase 
(white). The boundary sep¬ 
arating these phases (black 
solid line), corresponds to 
V = tv, where the bulk gap 
closes at some k. Two Hamil¬ 
tonians in the same phase are 
adiabatically connected. 
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Consider the effect of an adiabatic deformation of the Hamiltonian, indexed by 
some continuous parameter tf : 0 —?► 1, on the number Na—Nb- The Hamiltonian 
respect chiral symmetry, and its bulk energy gap exceeds 2e, for all values of d. 

The deformation can create zero energy states by bringing a nonzero energy edge 
state = 0)) to zero energy, Eo{d) = 0 for > t/' but not for d < d'. In that case, 
the chiral symmetric partner of 1%), which is F |'ff)((i)) up to a phase factor, has to 
move simultaneously to zero energy. The newly created zero energy edge states are 
Pa \^o{d')) andPg !%(£/')), which occupy sublattice A andB, respectively. Thus, the 
number Na — Nb is unchanged. 

The deformation can also bring a zero energy state Iff)) to energy £ > 0 at some 
d — d'. However, it must also create a chiral symmetric partner with energy £ < 0 at 
the same d'. This is the time reverse of the process of the previous paragraph: here, 
both Na and Nb must decrease by 1, and, again, Na — Nb is unchanged. 

The deformation can move nonzero energy states in or out of the — e < £ < e 
energy window. This obviously has no effect on the number Na — Nb- 

Due to the deformation, the wavefunction of a zero energy eigenstate can change 
so that it extends deeper and deeper into the bulk. However, because of the gap 
condition, zero energy states have to have wavefunctions that decay exponentially 
towards the bulk, and so this process cannot move them away from the edge. Thus, 
Na and Nb cannot be changed this way. 

The arguments above show that Na — Nb, the net number of edge states on sub¬ 
lattice A at the left edge, is a topological invariant. 


Bulk-boundary correspondence in the SSH model 

We have introduced two topological invariants for the SSH model: the winding num¬ 
ber V, of Eq. ( |1.38| l, and the net number of edge states, Na — Nb, of this section. The 
hrst one was obtained from the bulk Hamiltonian only, the second by looking at the 
low energy sector of the left edge. In the trivial case of the SSH model, v > w, both 
are 0; in the topological case, v < w, both are 1. This shows that we can use the bulk 
topological invariant (the winding number) to make simple robust predictions about 
the low-energy physics at the edge. This is a simple example for the bulk-boundary 
correspondence, a recurrent theme in the theory of topological insulators, which 
will reappear in various models in the forthcoming chapters. 


1.5.1 Bound states at domain walls 


Edge states do not only occur at the ends of an open chain, but also at interfaces 
between different insulating domains of the same chain. This can be understood via 


the fully dimerized limit, Eig. 1.8 There are two types of domain walls here: those 
containing single isolated sites, which host 0 energy states on a single sublattice (no 
onsite potentials are allowed), and those containing trimers. On a trimer, the odd 
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—3 m = 6 m = 7 


Fig. 1.8 A long, fully dimerized SSH chain with 3 domains. The boundaries between the domains, 
the “domain walls”, host zero energy eigenstates (yellow shading). These can be localized on a 
single site, as for the domain wall at « = 3, or on a superposition of sites, as the odd superposition 
of the ends of the trimer shared between the n = 6 and « = 7 unit cells. . 


superposition of the two end sites form a zero energy eigenstate. In the the example 
of Fig. fLS] this is 

W(|6,B)-|7,B))=0. (1.41) 

Note that, just as the edge states at the ends of the chain, these zero energy states at 
the interfaces have wavefunctions that take nonzero values on one sublattice only. 

From a perfect dimerized phase without domains it is only possible to germinate 
an even number of interfaces. This means that if one encounters a domain wall with 
a localized state on one sublattice then there will be another domain wall somewhere 
in the system - possibly at the system’s edge - with a localized state on the opposite 
sublattice. 

Consider a domain wall in an SSH system that is not in the fully dimerized limit. 
The wavefunctions of the edge states at the domain walls will penetrate to some 
small depth into the bulk, with exponentially decaying evanescent tails. For two 
domain walls at a distance of M unit cells, the two edge states on the walls will 
hybridize, form “bonding” and “anti-bonding” states. At half filling, of these only 
the negative energy eigenstate will be occupied. This state hosts a single electron, 
however, its wavefunction is localized with equal weight on the two domain walls. 
Hence each domain wall, when well separated from other domain walls and the ends 
of the chain, will carry half an electronic charge. This effect is sometimes referred 
to as “fractionalization” of the charge. 


1.5.2 Exact calculation of edge states 

The zero energy edge states of the SSH model can also be calculated exactly, even 
in the absence of translational invariance. Take an SSH model on N unit cells, with 
complex intracell and intercell hopping amplitudes, 

N N-i 

H ^ Y, {vm\m,B) {m,A\+h.c.) + Y {wm\m+l,A) {m,B\+h.c.). (1.42) 

m=l m=l 

We are looking for a zero energy eigenstate of this Hamiltonian, 
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H ^ {am\m,A)+bm\m,B)) =0. 

m=\ 

This gives us 2N equations for the amplitudes and bm, which read 


(1.43) 


m= ; v^am + Wmam+i =0; + Vm+i^m+i = 0; (1.44a) 

boundaries: VA^aA^ = 0; vibi=Q. (1.44b) 

The first set of equations is solved by 

m = 2,... ,N : a™ = 

m = 1,... ,A^-1 : bm = 


"i-l _y . 

n— «i; 

./=! 

(1.45) 

N-\ 

n 

(1.46) 


However, we also have to fulfil Eqs. ( |1.44b| l, which give 

bi =aN = 0 . 


(1.47) 


These equations together say that, in the generic case, there is no zero energy eigen¬ 
state, — bj^i — 0. 

Although there is no exactly zero energy state, Eqs. ( |1.45| ), ( |1.46| l and ( |1.47| l 
admit two approximate solutions in the thermodynamic limit, N if the average 
intercell hopping is stronger than the intracell hopping. More precisely, we define 
the “bulk average values”. 


log|v| = 


1 


iV-l 


N- 


Y ElogI 

^ m=l 


log|w| = 


1 


N-l 


N- 


Y ElogI 

^ m=l 


(1.48) 


Eqs. ( |1.45| ) and ( |1.46| l translate to 

with the localization length 


= (1.49) 

ki 


1 


log|w|-log|v| 


(1.50) 


If in the thermodynamic limit, the bulk average values, Eqs. ( |1.48| l make sense, and 
^ > 0, we have two approximate zero energy solutions. 


N 


|T) ='Y^am \m,A ); 

m=\ 


1^) = E 

m=\ 


(1.51) 
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with the coefficients Um and bm chosen according to Eqs. ( |1.45| l and ( |1.46| l, and ai, 
respectively, b{^, used to fix the norm of \L), respectively, \R). 


Hybridization of edge states 


The two states \L) and \R) hybridize under H to an exponentially small amount, and 
this induces a small energy splitting. We can obtain an estimate for the splitting, and 
the energy eigenstates, to a good approximation using adiabatic elimination of the 
other eigenstates. In this approximation, the central quantity is the overlap 




v^b^ 




{R\H\L) = aie 

with some (j) G [0,27r). The energy eigenstates are approximated as 

| L )+ e ''^'/2 \ R ) 


| 0 +) = 

| 0 -) = 


V2 

|L)-e''(’/2 \R) 

^ V2 




£_ = - 


aie 


(1.52) 

(1.53) 

(1.54) 


The energy of the hybridized states thus is exponentially small in the system size. 


Problems 

1.1. Higher winding numbers 

The SSH model is one-dimensional in space, and has a two-dimensional internal 
Hilbert space. Construct a lattice model that has these properties of the SSH model, 
but which has a bulk winding number of 2. Generalize the construction for an arbi¬ 
trary integer bulk winding number. 

1.2. Complex-valued hopping amplitudes 

Generalize the SSH model in the following way. Assume that the hopping ampli¬ 
tudes V = and w = \w\e“l’"' are complex, and include a third complex-valued 

hopping amplitude z = \z\e“l’' between the states \m,A) and \m+ l,B) for every m. 
Provide a specific example where the tuning of one of the phases changes the bulk 
winding number. 

1.3. A possible generalization to two dimensions 

Consider a two dimensional generalization of the SSH model. Take parallel copies 
of the SSH chain and couple them without breaking chiral symmetry. What will 
happen with the edge states? 














Chapter 2 

Berry phase, Chern number 


To describe the theory of topological band insulators we will use the language of 
adiabatic phases. In this chapter we review the basic concepts: the Berry phase, the 
Berry curvature, and the Chern number. We further describe the relation between the 
Berry phase and adiabatic dynamics in quantum mechanics. Finally, we illustrate 
these concepts using the two-level system as a simple example. 

For pedagogical introductions, we refer the reader to Berry’s original paper 1^, 
and papers from the Americal Journal of Physics ifTSlflTl . For the application to solid 
state physics, we will mostly build on Resta’s lecture note ll26l . and Niu’s review 
paperll3^. 


2.1 Discrete case 

The subject of adiabatic phases is strongly related to adiabatic quantum dynamics, 
when a Hamiltonian is slowly changed in time, and the time evolution of the quan¬ 
tum state follows the instantaneous eigenstate of the Hamiltonian. In that context, 
as time is a continuous variable and the time-dependent Schrodinger equation is a 
differential equation, the adiabatic phase and the related concepts are expressed us¬ 
ing differential operators and integrals. We will arrive to that point later during this 
chapter; however, we start the discussion using the language of discrete quantum 
states. Besides the conceptual simplicity, this language also offers an efficient tool 
for the numerical evaluation of Chern number, which is an important topological 
invariant in the context of two-dimensional electron systems. 


2.1.1 Relative phase of two nonorthogonal quantum states 

In quantum mechanics, the state of a physical system is represented by an equiva¬ 
lence class of vectors in a Hilbert space: a multiplication by a complex phase factor 
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does not change the physical content. A gauge transformation is precisely such a 
multiplication: 


If') withae [0,2;r). (2.1) 

In that sense, the phase of a vector jf') does not represent physical information. We 
can try to dehne the relative phase 712 of two nonoithogonal states and as 

7i2 = -arg('f'i|'f'2), (2.2) 


where arg(z) denotes the phase of the complex number z, with the specihcation that 
arg(z) G (—?t]. Clearly, the relative phase 712 fulhls 


m\m' 


(2.3) 


However, the relative phase is not invariant under a local gauge transformation, 

\'Fj) ^ g-irn+iiai-ai)^ (2.4) 


2.1.2 Berry phase 

Take N >3 states in a Hilbert space, order them in a loop, and ask about the phase 
around the loop. As we show below, the answer - the Berry phase - is gauge invari¬ 
ant. For states j’f;), with j = 1,2,...,A, and for the ordered list L = (1,2,...,A) 
which dehne the loop, shown in Fig. |2.1[ the Berry phase is dehned as 

jL = -arge-‘(y‘2+l^3+...+rvi) = -arg(('f'i | -ft) Cft | f'a) ■ ■ ■ (fV | f'l)) • (2.5) 

To show the gauge invariance of the Berry phase, it can be rewritten as 

Yl = - argTr (l-Fi) {W, \ j-ft) Cftl • • • \^n) (fVl) ■ (2.6) 

Here, we expressed the Berry phase Yl using projectors that are themselves gauge 
invariant. 

Even though the Berry phase is not the expectation value of some operator, it is 
a gauge invariant quantity, and as such, it can have a direct physical signihcance. 
We will hnd such a signihcance, but hrst, we want to gain more intuition about its 
behaviour. 



2.1 Discrete case 
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Fig. 2.1 Berry phase, Berry flux and Berry curvature for discrete quantum states, (a) The Berry 
phase /z, for the loop L consisting of N = 3 states is defined from the relative phases 712, yii, Ya- 
(b) The Berry phase of a loop defined on a lattice of states can be expressed as the sum of the Berry 
phases Fyi and ^2,1 of the plaquettes enclosed by the loop. The plaquette Berry phase is also 
called Berry flux. 


2.1.3 Berry flux 

Consider a Hilbert space of quantum states, and a finite two-dimensional square lat¬ 
tice with points labelled by n, m G Z, 1 <n<N, and 1 <m< M. Assign a quantum 
state 1^5, ,,,) from the Hilbert space to each lattice site. Say you want to know the 
Berry phase of the loop L around this set, 


Yl = -argexp 


C N-\ M-\ 

y^. 7(n,l),(n+l,l) “f ^ Y{N,m),(N,m+\) 

n=\ m=l 

N-\ M-1 

+ 52 7(«+l,M),(n,M) + 52 7(l,m+l).(l.m) 
n=\ m=\ 


(2.7) 


as shown in Fig. |2.1| Although the Berry phase is a gauge invariant quantity, cal¬ 
culating it according to the recipe above involves multiplying together many gauge 
dependent complex numbers. The alternative route, via Eq. ( |231 l, involves multli- 
plying gauge independent matrices, and then taking the trace. 

There is a way to break the calculation of the Berry phase of the loop down to 
a product of gauge independent complex numbers. To each plaquette (elementary 
square) on the grid, with indexing the lower left corner, we dehne the Berry 
flux F„ ,fi of the plaqette using the sum of the relative phases around its boundary. 


Fnm — -argexp 

,(«+!,m) 7(«+l,m),(n+l,m+l) 

7(n,m+l),(n.m))] ; (^■8) 

for n = 1,... and m = 1,... ,M. Note that the Berry flux is itself a Berry phase 
and is therefore gauge invariant. Alternatively, we can also write 


Ffjm — ^g ( {^n,m \ | 

{%+l.m+l I \ %,n,) ) , (2.9) 




24 


2 Berry phase, Chern number 


Now consider the product of all plaquette phase factors e 




n n ^ 

n=\ m=\ 


= exp 


n=\ m=\ 


— exp 


«=1 m=l 


' 7(n+l.m),(rt+l,m+l) “f 7(«+l,m+l),(«,m+l) “f 7(«,m+l),(n,m)) 


( 2 . 10 ) 


Each internal edge of the lattice is shared between two plaquettes, and therefore 
occurs twice in the product. However, since we fixed the orientation of the plaquette 
phases, these two contributions will always be complex conjugates of each other, 
and cancel each other. Therefore the exponent in the right-hand-side of Eq. ( |2.10| l 
simplihes to the exponent appearing in Eq. (|2.7|l, implying 


exp 






n=l m=l 


( 2 . 11 ) 


This result is reminsicent of the Stokes theorem connecting the integral of the curl 
of a vector field on an open surface and the line integral of the vector held along 
the boundary of the surface. In Eq. ( |2.1 l| l, the sum of the relative phases, i.e., the 
Berry phase yi, plays the role of the line integral, whereas the double sum of the 
Berry Huxes plays the role of the surface integral. There is an important difference 
with respect to the Stokes theorem, namely, the equality of the total Berry Hux and 
the Berry phase is not guaranteed: Eq. ( |2.11| i only tells that they are either equal or 
have a difference of 2n times an integer. 


2.1.4 Chern number 

Consider states in a Hilbert space ari'anged on a grid as above, with n,m G Z, 

1 <n<N, and I <m<M, but now imagine this grid to be on the surface of a torus. 
We use the same definition for the Berry flux per plaquette as in ( |2.9| l, but now with 
n mod N -f 1 in place of n -f 1 and m mod M -f 1 in place of m -|-1. 

The product of the Berry flux phase factors of all plaquettes is now 1, 

M N 

= ( 2 . 12 ) 

m=\ n=\ 

The same derivation can be applied as for Eq. p.l l| i above, but now every edge is 
an internal edge, and so all contributions to the product cancel. 

The Chern number Q associated to our structure is defined via the sum of the 
Berry fluxes of all the plaquettes forming the closed torus surface: 

Q = 

y,yy, 


(2.13) 
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The fact that the Chern number Q is defined via the gauge invariant BerTy fluxes 
ensures that Q itself is gauge invariant. Furthermore, taking the arg of Eq. ( |2.12| i 
proves that the Chern number Q is an integer. 

It is worthwhile to look a little deeper into the discrete formula for the Chern 
number. We can define modified Berry fluxes F„m as 

(2.14) 

Since each edge is shared between two neighboring plaquettes, the sum of the mod¬ 
ified Berry fluxes over all plaquettes vanishes, 

M N 

( 2 . 15 ) 

m=l n=\ 

If — TT < Fnm < Tt, then we have F„m = Fnm- However, Fnm can be outside the range 
then as the logarithm is taken in Eq. Fnm is taken back into [—TT, n) 
by adding a (positive or negative) integer multiple of 2n. In that case, we say the 
plaquette nm contains a number Qnm € Z of vortices, with 

= (2.16) 
IK 

We have found a simple picture for the Chern number; The Chern number Q, that 
is, the sum of the Berry fluxes of all the plaquettes of a closed surface, is the number 
of vortices on the surface. 


= (2.17) 

nm nm 

Although we proved it here for the special case of a torus, the derivation is easily 
generalized to all orientable closed surfaces. We focused on the torus, because this 
construction can be used as a very efficient numerical recipe to discretize and calcu¬ 
late the (continuum) Chern number of a 2-dimensional insulator jTJl, to be defined 
in Sect. 12.2.^ 


2.2 Continuum case 

We now assume that instead of a discrete set of states, {we have a continuum, 
l'E(R)), where the R’s are elements of some D-dimensional parameter space 
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2.2.1 Berry connection 


We take a smooth directed path i.e., a curve in the parameter space 




(2.18) 


We assume that all components of |'f'(R)) are smooth, at least in an open neighbor¬ 
hood of the the curve The relative phase between two neighbouring states on the 
curve corresponding to the parameters R and R + c/R, is 


(T(R)|T(R + r/R)) . 

\{'F{R)\'F{R + dR)W 


47=/('T(R)|VR|'f'(R)).r/R, (2.19) 


obtained to first order in dR 0. The quantity multiplying dR on the right-hand 
side defines the Berry connection. 


A(R) = i ('f'(R) I VR'f'(R)) = -Im ('f'(R) | VR'f'(R)). (2.20) 

Here | VR'f'(R)) is defined by requiring for every Hilbert space vector |^>), that 

(0| VR'f'(R))=VR(0|'T(R)). (2.21) 

The second equality in Eq. ( |2.20| ) follows from the conservation of the norm, 
VR('f'(R)|'T(R))=0. 

We have seen in the discrete case that the relative phase of two states is not gauge 
invariant; neither is the Berry connection. Under a gauge transformation, it changes 
as 

l'f'(R)) ^ l'f'(R)) : A(R) ^ A(R) - VRa(R). (2.22) 


2.2.2 Berry phase 

Consider a closed directed curve in parameter space. The Berry phase along the 
curve is defined as 


7(‘^) = -argexp 


A dR 


(2.23) 


The Berry phase of a closed directed curve is gauge invariant, since it can be in¬ 
terpreted as a limiting case of the discrete Berry phase, via Eqs. ( |2.20| i, ( |2.19| l, and 
and the latter has been shown to be gauge invariant. 
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2.2.3 Berry curvature 

As in the discrete case above, we would like to express the gauge invariant Berry 
phase as a surface integral of a gauge invariant quantity. This quantity is the Berry 
curvature. Similarly to the discrete case, we consider a two-dimensional parameter 
space, and for simplicity denote the parameters as x and y. We take a simply con¬ 
nected region in this two-dimensional parameter space, with the oriented bound¬ 
ary curve of this surface denoted by and consider the continuum Berry phase 
corresponding to the boundary. 


Smoothness of the manifold of states 


Before relating the Berry phase to the Berry curvature, an important note on the 
manifold |T'(R)) of considered states is in order. From now on, we consider a man¬ 
ifold of states, living in our two-dimensional parameter space, that is smooth, in 
the sense that the map R |T'(R)) (T'(R)| is smooth. Importantly, this condition 
does not necessarily imply that that the function R |'f'(R)), also referred to as a 
gauge describing our manifold, is smooth. (For further discussion and examples, see 
Sect. 2.5.1 ) Nevertheless, even if the gauge R |'f'(R)) is not smooth in a point Rq 
of the parameter space, one can always find an alternative gauge |'F'(R)) which is 
(i) locally smooth, that is, smooth in the point Rq, and (ii) locally generates the same 
map as |'f'(R)), that is, for which |'F'(R)) ('F'(R)| = |'f'(R)) (’f'(R)| in an infinites¬ 
imal neighborhood of Rq. Let us formulate an intuitive argument supporting the 
latter claim using quantum-mechanical perturbation theory. Take the Hamiltonian 
H(R) = — l'F(R)) ('f'(R)l, which can be substituted in the infinitesimal neighbor¬ 
hood of Rq with H(Ro -t- AR) = H(Ro) -f AR- (VH)(Ro). According to first-order 
perturbation theory, the ground state of the latter is given by 


D 

|'F'(Ro + AR)) = l'f'(Ro))- £ |f'«(Ro)) ('f'„(Ro)|AR-(Vi7)(Ro)|'F(Ro)), 

n=2 

(2.24) 

where the states |T(,(Ro)) (n — 2,3,... ,D), together with |'F(Ro)), form a basis of 
the Hilbert space. On the one hand, Eq. p.24| l defines a function that is smooth in 
Rq, hence the condition (i) above is satisfied. On the other hand, as |'F'(Ro -f AR)) 
is the ground state of H(Ro 4- AR), condition (ii) is also satisfied. 


Berry phase and Berry curvature 

Now return to our original goal and try to express the Berry phase as a surface 
integral of a gauge invariant quantity; to this end, we start by relating the Berry 
phase to its discrete counterpart; 
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, lim Jd^, 


(2.25) 


where we discretize the parameter space using a square grid of steps Ax, Ay, and 
express the integral as the discrete Berry phase of a loop approximating d,^, in 
the limit of an infinitesimally fine grid. Then, from Eq. p.25| l and the Stokes-type 
theorem in Eq. (|2.1 l|l, we obtain 


exp 


A c/R 




= lim e 'E""'-'"™, 

Ax,Ay^0 


(2.26) 


where the nm sum goes for the plaquettes forming the open surface Further¬ 
more, let us take a guage |'f''(R)) and the corresponding Berry connection A' that is 
smooth in the plaquette nm; this could be |'f'(R)) and A if that was already smooth. 
Then, due to the gauge invariance of the Berry flux we have 

e-'P'"" = e-‘FL, (2.27) 


where is the Berry flux corresponding to the locally smooth gauge. Furthermore, 
in the limit of an infinitely fine grid it holds that 


Km = (x„ + Ax+A'y (^Xn+uym + 4y 

f \ . .if ^y\. 

I 2 I {^n-iym 2 / 

Taylor expansion of the Berry connection around R„,„ = (^x„ + + 

order yields 


(2.28) 

to first 


Km = [dxA'yiR.^m) - dyA',{R„„,)] AxAy. (2.29) 


Thereby, with the definition of the Berry curvature as 

B= lim (2.30) 

Ax.Ay^O AxAy 

we obtain a quantity that is gauge invariant, as it is defined via the gauge invariant 
Berry flux, and is related to the Berry connection via 


B = d.MyiRnm) - dyA!^[Rnm). (2.31) 

We can rephrase Eq. p.29| l as follows: the Berry flux for the nm plaquette is ex¬ 
pressed as the product of the Berry curvature on the plaquette and the surface area 
of the plaquette. 

Substituting Eqs. ( |2.27| i and ( |2.29| l into Eq. ( |2.26| l yields 
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exp 


A dR 


d,^ 


= exp 


B{x,y)dxdy 


(2.32) 


which is the continuum version of the result ( |2.11| ). Equation ( |2.32| l can also be 
rephrased as 


Y{d^) = 


(2.33) 


A special case where the usual Stokes theorem works 


A shortcut towards a stronger result than Eq. ( |2.32| l is offered in the special case 
when l'E(R)) is smooth in the neighborhood of the open surface Then, a direct 
application of the two-dimensional Stokes theorem implies 



dR = 



— dyAx)dxdy = / Bdxdy 


(2.34) 


Summarizing Eqs. ( |2.32| i and ( |2.34[ ), we can say that line integral of the Berry con¬ 
nection equals the surface integral of the Berry curvature if the set of states |'E(R)) 
is smooth in the neighborhood of but they might differ with an integer multiple 
of 2;r otherwise. 


The case of the three-dimensional parameter space 

Let us briefly discuss also the case of a 3D parameter space. This will be particularly 
useful in the context of two-level systems. Starting with the case when the gauge 
l'E(R)) on the two-dimensional open surface embedded in the 3D parameter 
space is smooth in the neighborhood of we can directly apply the 3D Stokes 
theorem to convert the line integral of A to the surface integral of the curl of A to 
obtain 


</ A dR= f R dS, (2.35) 

where the Berry curvature is defined as the vector field B(R) via 

B(R) = VrxA(R), (2.36) 

which is gauge invariant as in the two-dimensional case. Even if |'E(R)) is not 
smooth on the relation 

Y{d^) = - = - arge^'-^-^® ''® (2.37) 

holds, similarly to the two-dimensional result Eq. ( |2.32| l. 

Note furthermore that the Berry phase Y{d^) is not only gauge invariant, but 
also invariant against continuous deformations of the two-dimensional surface ^ 
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embedded in 3D, as long as the Berry curvature is smooth everywhere along the 
way. 

We also remark that although we used the three-dimensional notation here, but 
the above results can be generalized for any dimensionality of the parameter space. 

The notation A and B for the Berry connection and Berry curvature suggest that 
they are much like the vector potential and the magnetic held. This is a useful anal¬ 
ogy, for instance, VrB = 0, from the dehnition ( |2.361 l. Nevertheless, it is not true 
that in every problem where the Berry curvature is nonzero, there is a physical mag¬ 
netic held. 


2.2.4 Chern number 

In the discrete case, we dehned the Chern number as a sum of Berry huxes for a 
square lattice living on a torus. Here, we take a continuum parameter space that has 
the topology of a torus. The motivation is that certain physical parameter spaces in 
fact have this torus topology, and the corresponding Chern number does have physi¬ 
cal signihcance. One example will be the Brillouin zone of a two-dimensional lattice 
representing a solid crystalline material. The Brillouin zone has a torus topology, as 
the momentum vectors {kx,ky), {kx + 27t,ky), and {kx,ky + 2n) are equivalent. 

Quite naturally, in the continuum dehnition of the Chern number, the sum of 
Berry huxes is replaced by the surface integral of the Berry curvature; 



(2.38) 


As this can be interpreted as a continuum limit of the discrete Chern number, it 
inherits the properties of the latter; the continuum Chern number is a gauge invariant 
integer. 

For future reference, let us lay down the notation to be used for calculating the 
Chern numbers of electronic energy bands in two-dimensional crystals. Consider 
a square lattice for simplicity, which has a square-shaped Brillouin zone as well. 
Our parameter space ^ is the two-dimensional Brillouin zone now, which has a 
torus topology as discussed above. The parameters are the Cartesian components 
kx^ky G [—TTjTt) of the momentum vector k. The electronic energy bands and the 
corresponding electron wavefunctions can be obtained from the bulk momentum- 
space Hamiltonian H{kx,ky). The latter dehnes the Schrodinger equation 


H(k)|M„(k)) =£'„(k)|M„(k)), 


(2.39) 


where n = 1,2,... is the band index, which has as many possible values as the 
dimension of the Hilbert space of the internal degree freedom of our lattice model. 
Note that dehning the Berry connection, the Berry curvature and the Chern number 
for the «th band is possible only if that band is separated from the other bands by 
hnite energy gaps. The Berry connection of the nth band, in line with the general 
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definition (|2.20|l, reads 


(k) = i {u„ (k) I dkj \u„ (k)) , for j = x,y. 


(2.40) 


The Chern number of the nth band, in correspondence with Eqs. ( |2.38| l and ( |2.31[ ), 
reads 



(2.41) 


Certain approximations of the band-structure theory of electrons provide low¬ 
dimensional momentum-space Hamiltonians that can be diagonalized analytically, 
allowing for an analytical derivation of the Chern numbers of the electronic bands. 
More often, however, the electronic wave functions are obtained from numerical 
techniques on a finite-resolution grid of {kx,ky) points in the Brillouin zone. In that 
case, the Chern number of a chosen band can still be effectively evaluated using the 
discrete version of its definition ( |2.13| l. 

The Chern number of a band of an insulator is a topological invariant in the fol¬ 
lowing sense. One can imagine that the Hamiltonian describing the electrons on 
the lattice is deformed adiabatically, that is, continuously and with the energy gaps 
separating the nth band from the other bands kept open. In this case, the Berry cur¬ 
vature varies continuously, and therefore its integral for the Brillouin zone, which is 
the Chern number, cannot change as the value of the latter is restricted to integers. If 
the deformation of the crystal Hamiltonian is such that some energy gaps separating 
the nth band from a neighboring band is closed and reopened, that is, the deforma¬ 
tion of the Hamiltonian is not adiabatic, then the Chern number might change. In 
this sense, the Chern number is a similar topological invariant for two-dimensional 
lattice models as the winding number is for the one-dimensional SSH model. 


2.3 Berry phase and adiabatic dynamics 

In most physical situations of interest, the set of states whose geometric features 
(Berry phases) we are interested in are eigenstates of some Hamiltonian H. Take a 
physical system with D real parameters that are gathered into a formal vector R = 
{Ri ,R2, ■ ■ ■ ,Rd)- The Hamiltonian is a smooth function H{R) of the parameters, at 
least in the region of interest. We order the eigenstates of the Hamiltonian according 
to the energies ^^(R), 


H{R)\n{R))=En{R)\n{R)}. 


(2.42) 


We call the set of eigenstates |n(R)) the snapshot basis. 

The definition of the snapshot basis involves gauge fixing, i.e., specifying the 
otherwise arbitrary phase prefactor for every |n(R)). This can be a tricky issue: 
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even if in theory a gauge exists where all elements of the snapshot basis are smooth 
functions of the parameters, this gauge might be very challenging to construct. 

We consider the following problem. We assume that the system is initialized with 
R = Rq and in an eigenstate |n(Ro)) that is in the discrete part of the spectrum, i.e., 
£'„(R) -£'„_i(R) and£'„+i(R) -£'„(R) are nonzero. At time f = 0 we thus have 

R(f = 0) = Ro; I v/(f = 0)) = |n(Ro)). (2.43) 

Now assume that during the time t —0 T the parameters R are slowly changed: 
R becomes R(f), and the values of R(f) dehne a continuous directed curve Also, 
assume that |n(R)) is smooth along the curve The state of the system evolves 
according to the time-dependent Schrodinger equation: 

ij^\xif{t))=H{R{t))\xif{t)). (2.44) 

Further, assume that R is varied in such a way that at all times the energy gaps 
around the state |n(R(f))) remain hnite. We can then choose the rate of variation of 
R(f) along the path ^ to be slow enough compared to the frequencies corresponding 
to the energy gap, so the adiabatic approximation holds In that case, the system 
remains in the energy eigenstate |n(R(f))), only picking up a phase. We are now 
going to find this phase. 

By virtue of the adiabatic approximation, we take as Ansatz 

\xj/{t)) = e'A(0e-'7oA,(R(^')K |„(R(f))). (2.45) 

For better readability, in the following we often drop the t argument where this leads 
to no confusion. The time derivative of Eq. ( |2.45| l reads 

i±\xif(t))=e<r'^e-<!oMm)dr' |„(r))+£„(r) |„(r)). 

(2.46) 

To show what we mean by | ^n(R(f))), we write it out explicitly in terms of a hxed 
basis, that of the eigenstates at R = Rq: 

|n(R))=^c™(R)|m(Ro)); (2.47) 

m 

fiTl /iJi. 

\in{R{t))) = — . |Vr«(R)) = — I:Vrc„,(R) |m(Ro)). (2.48) 

We insert the Ansatz ( |2.45| l into the right hand side of the Schrodinger equa¬ 
tion ( |2.44| ), use the snapshot eigenvalue relation ( |2.42| i, simplify and reorder the 
Schrodinger equation, and obtain 
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Multiplying from the left by (n(R)|, and using Eq. ( |2.48| l, we obtain 

= i (n(R(0) I ®«(R(0)) = («(R) I Vr«(R)) ■ (2-50) 

We have found that for the directed curve ^ in parameter space, traced out by 
R(f), there is an adiabatic phase 7n('^), which reads 

7„i^) = Ji{n{R)\W^n(R))dR. (2.51) 

-if 


A related result is obtained after a similar derivation, if the parameter space of the R 
points is omitted and the snapshot basis |n(f)) is parametrized directly by the time 
variable. Then, the adiabatic phase is 

Yti{t) = j i{n{t')\ dtin{t')')dt'. (2.52) 

Equation ( |2.51| l allows us to formulate the key message of this section as the 
following. Consider the case of an adiabatic and cyclic change of the Hamiltonian, 
that is, when the curve ^ is closed, implying R(r) = Rq. In this case, the adiabatic 
phase reads 


y„{^) = ji{n{R)\V^n{R))dR. 

-f 


(2.53) 


Therefore, the adiabatic phase picked up by the state during a cyclic adiabatic 
change of the Hamiltonian is equivalent to the Berry phase corresponding to the 
closed oriented curve representing the Hamiltonian’s path in the parameter space. 

Two further remarks are in order. Eirst, on the face of it, our derivation seems to 
do too much. It seems that we have produced an exact solution of the Schrodinger 
equation. Where did we use the adiabatic approximation? In fact, Eq. ( |2.50| l does not 
imply Eq. ( |2.49| l. Eor the more complete derivation, showing how the nonadiabatic 
terms appear, see ns. 

The second remark concerns the measurability of the Berry phase. The usual way 
to experimentally detect phases is by an interferometric setup. This means coher¬ 
ently splitting the wavefunction of the system into two parts, taking them through 
two adiabatic trips in parameter space, via R(f) and R'(f), and bringing the parts 
back together. The interference only comes from the overlap between the states; it 
is maximal if |n(R(7’))) = |n(R'(7’))), which is typically ensured if R(7’) = R'(7’). 
The difference in the adiabatic phases y„ and is the adiabatic phase associ¬ 
ated with the closed loop which is the path obtained by going forward along 
f = 0 —>■ r ; R(f), then coming back along f = T —>• 0 ; R'(f). 
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2.4 Berry’s formulas for the Berry curvature 

Berry providedll&l two practical formulas for the Berry curvature. Here we present 
them in a form corresponding to a three-dimensional parameter space. To obtain 
the two-dimensional case, where the Berry curvature B is a scalar, one can identify 
the latter with the component of the three-dimensional case treated below; for 
generalization to higher than 3 dimensions, see the discussion in Berry’s paperl^. 
First, 


Bj = -Im Ejkt dk {n \ din) = -Im Sjki {d^n \ din) -f 0, (2.54) 

where the second term is 0 because d^di — d[dk but Sju = —£jik- 

To obtain Berry’s second formula, inserting a resolution of identity in the snap¬ 
shot basis in the above equation, we obtain 

= -Im Y, (^« I «') X («' I > (2-55) 

n'^n 


where the parameter set R is suppressed for brevity. The term with n' = n is omitted 
from the sum, as it is zero, since because of the conservation of the norm, (Vn | n) = 
— {n I Vn). To calculate {n' \ Vn), start from the dehnition of the eigenstate |n), act 
on both sides with V, and then project unto |n'); 

H\n)=E„\n); (2.56) 

{VH)\n) +H\Vn) = {VE„)\n) +En\Vn)- (2.57) 

(n'l VH|n)-t-(n'|i/|Vn) =0-f£„(n'| Vn) . (2.58) 

Act with El towards the left in Eq. ( |2.58| l, rearrange, substitute into ( |2.55| l, and you 
should obtain the second form of the Berry curvature, which is manifestly gauge 
invariant; 


b(«) 


-Im Y 

n'^n 


(n|VH|n') X (n'|Vi/|n) 

[En-En')^ 


(2.59) 


This shows that the monopole sources of the Berry curvature, if they exist, are the 
points of degeneracy. 

A direct consequence of Eq. ( |2.59| l, is that the sum of the Berry curvatures of all 
eigenstates of a Hamiltonian is zero. If all the spectrum of H(R) is discrete along a 
closed curve then one can add up the Berry phases of all the energy eigenstates. 
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y- (w| VRj/|«') X («'| VrH|w) 
« n'^n {En-En'Y 



+ {n' \ VrH \n) X {n\ VrH 


(2.60) 


The last equation holds because a x b = — b x a for any two vectors a,b. 

2.5 Example: the two-level system 

So far, most of the discussion on the Berry phase and the related concepts have been 
kept rather general. In this section, we illustrate these concepts via the simplest 
nontrivial example, that is, the two-level system. 

2.5.1 No continuous global gauge 

Consider a Hamiltonian describing a two-level system: 


H(d) = dxdx + dyOy + = d • ff, 


(2.61) 


with d = {dx,dy,dY) € K^\{0}. Here, the vector d plays the role of the parameter 
R in of preceding sections, and the parameter space is the punctured 3D Euclidean 
space K^\{0}, to avoid the degenerate case of the energy spectrum. Note the ab¬ 
sence of a term proportional to Oq: this would play no role in adiabatic phases. Be¬ 
cause of the anticommutation relations of the Pauli matrices, the Hamiltonian above 
squares to a multiple of the identity operator, H{dY = d^Ob. Thus, the eigenvalues 
of H{d) have to have absolute value |d|. 

A practical graphical representation of H(d) is the Bloch sphere, shown in 
Fig.|2.2| The spherical angles 9 G [0,7r) and (p G [0,27r) are defined as 




(2.62) 


We denote the two eigenstates of the Hamiltonian i/(d) by |-|-d) and |—d), with 


H(d) |±d) = ±|d| |±d> • 


(2.63) 


These eigenstates depend on the direction of the 3-dimensional vector d, but not on 
its length. The eigenstate with £ = -f |d| of the corresponding Hamiltonian is: 
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(a) 


A d. 


ib) 



d. 


'X 



'X 


Fig. 2.2 The Bloch sphere. A generic traceless gapped two-level Hamiltonian is a linear combi¬ 
nation of Pauli matrices, H{d) = d - This can be identified with a point in R^\{0}. The eigenen- 
ergies are given by the distance of the point from the origin, the eigenstates depend only on the 
direction of the vector d, i.e., on the angles 6 and (p, as defined in subfigure (a) and in Eq. l |2.62^ 
The Berry phase of a closed curve is half the area enclosed by the curve when it is projected 
onto the surface of the Bloch sphere. 



(2.64) 


while the eigenstate with £ = — |d| is |— d) = |-l--d)- The choice of the phase 

factors a and j3 above corresponds to hxing a gauge. We will now review a few 
gauge choices. 

Consider hxing a{9,(p) = 0 for all 9,(p. This is a very symmetric choice, in this 
way in formula ( |2.64| i, we hnd 9/2 and (p/2. There is problem, however, as you 
can see if you consider a full circle in parameter space; at any hxed value of 9, let 
(p = 0 ^ 2n. We should come back to the same Hilbert space vector, and we do, but 
we also pick up a phase of — 1. We can either say that this choice of gauge led to a 
discontinuity at ^ = 0, or that our representation is not single-valued. We now look 
at some attempts at hxing these problems, to hnd a gauge that is both continuous 
and single valued. 

As a hrst attempt, let us hx a = ^/2; denoting this gauge by subscript S, we have 



(2.65) 


The phase prefactor now gives an additional factor of — 1 as we make the circle in (p 
at hxed 9, and so it seems we have a continuous, single valued representation. There 
are two tricky points, however; the North Pole, 0 = 0, and the South Pole, 9 — n. 
At the North Pole, |(0,0, l))y = (1,0) no problems. This gauge is problematic at the 
South Pole, however (which explains the choice of subscript); there, |(0,0, —1))^ = 
(0,e''*’), the value of the wavefunction depends on which direction we approach the 
South Pole from. 
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We can try to solve the problem at the South Pole by choosing a = — ^/2, which 
gives us 


l+d)Ar — 


/e-''Pcos0/2\ 
sin 0/2 J ' 


( 2 . 66 ) 


As you can probably already see, this representation runs into trouble at the North 
Pole: 1(0,0, l))^=(e-'^0). 

We can try to overcome the problems at the poles by taking linear combinations 
of |+d)s and l+d)Ar5 with prefactors that vanish at the South and North Poles, re¬ 
spectively. A family of options is: 


0 0 

l+d);^ =e'^sin-|+d)5-Ecos-|+d)^ 

_/ cos |(cos I+sin \ 

\^sin|e"P(cos | -|-sin J ' 


(2.67) 

( 2 . 68 ) 


This is single valued everywhere, solves the problems at the Poles. However, it has 
its own problems: somewhere on the Equator, at 0 = 7r/2, ^ j ± ;r, its norm 
disappears. 

It is not all that surprising that we could not find a well-behaved gauge: there is 
none. By the end of this chapter, it should be clear, why. 


2.5.2 Calculating the Berry curvature and the Berry phase 

Consider the two-level system as defined in the previous section. Take a closed curve 
^ in the parameter space IR^\{0}. We are going to calculate the Berry phase j- of 
the |—d) eigenstate on this curve: 

= ^A(d)c/d, (2.69) 

with the Berry vector potential defined as 

A(d) = /(-d|Vd|-d). (2.70) 

The calculation becomes straightforward if we use the Berry curvature, 

B(d)=VdxA(d); 

7 -(^)= [ Rid)dS^, 

Jy 


(2.71) 

(2.72) 
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where is any surface whose boundary is the loop (Alternatively, it is a worth¬ 
while exercise to calculate the Berry phase directly in a fixed gauge, e.g., one of the 
three gauges of the previous chapter.) 

Specifically, we make use of Berry’s gauge invariant formulation ( |2.59| l of the 
Berry curvature, derived in the last chapter. In the case of the generic two-level 
Hamiltonian (|2.61|l, Eq. (|2.59|l gives 


B±(d) = 


(2.73) 


with 


VdH = a. (2.74) 

To evaluate ( |2.73| l, we choose the quantization axis parallel to d, thus the eigenstates 
simply read 



The matrix elements can now be computed as 

and similarly, 

(-|cTvl+) = /; 

(-|cT.|+)=0. 


So the cross product of the vectors reads 


(-|a|+)x(+|a|-) 



This gives us for the Berry curvature. 


B±(d) 


|d| 2d2 ■ 



(2.75) 


(2.76) 


(2.77) 

(2.78) 


(2.79) 


(2.80) 


We can recognize in this the field of a pointlike monopole source in the origin. 
Alluding to the analog between the Berry curvature and the magnetic field of elec¬ 
trodynamics (both are derived from a “vector potential”) we can refer to this field, 
as a “magnetic monopole”. Note however that this monopole exists in the abstract 
space of the vectors d and not in real space. 

The Berry phase of the closed loop in parameter space, according to Eq. ( |2.72[ ), 
is the flux of the monopole held through a surface whose boundary is ^. It is 















2.5 Example: the two-level system 


39 


easy to convince yourself that this is half of the solid angle subtended by the curve, 




(2.81) 


In other words, the Berry phase is half of the area enclosed by the image of 


projected onto the surface of the unit sphere, as illustrated in Fig. 2.2 


What about the Berry phase of the other energy eigenstate? From Eq. ( |2.73| ), the 
corresponding Berry curvature B+ is obtained by inverting the order of the factors 
in the cross product: this flips the sign of the cross product. Therefore the Berry 
phases of the ground and excited state fulfil the relation 


7+(^) = -7-('^)- (2.82) 

One can see the same result on the Bloch sphere. Since (+1 —) = 0, the point corre¬ 
sponding to |—) is antipodal to the point corresponding to j-f-). Therefore, the curve 
traced by the |—) on the Bloch sphere is the inverted image of the curve traced 
by |+). These two curves have the same orientation, therefore the same area, with 
opposite signs. 


2.5.3 Two-band lattice models and their Chern numbers 

The simplest case where a Chern number can arise is a two-band system. Consider 
a particle with two internal states, hopping on a two-dimensional lattice. The two 
internal states can be the spin of the conduction electron, but can also be some 
sublattice index of a spin polarized electron. In the translation invariant bulk, the 
wave vector k = {kx,ky) is a good quantum number, and the Hamiltonian reads 

H(k)=d(k)ff, (2.83) 

with the function d(k) mapping the each point of the Brillouin Zone to a 3D vector. 
Since the Brillouin zone is a torus, the endpoints of the vectors d(k) map out a 
deformed torus in K^\{0}. This torus is a directed surface: its inside can be painted 
red, its outside, blue. 

The Chern number of |—) (using the notation of Sect. of \ui (k))) is the flux 
of B_ (d) through this torus. We have seen above that B_ (d) is the magnetic held of 
a monopole at the origin d = 0. If the origin is on the inside of the torus, this flux is 
H-1. If it is outside of the torus, it is 0. If the torus is turned inside out, and contains 
the origin, the flux is -1. The torus can also intersect itself, and therefore contain the 
origin any number of times. 

One way to count the number of times is to take any line from the origin to infin¬ 
ity, and count the number of times it intersects the torus, with a H-1 for intersecting 
from the inside, and a -1 for intersecting from the outside. The sum is independent 
of the shape of the line, as long as it goes all the way from the origin to infinity. 
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Problems 


2.1. Discrete Berry phase and Bloch vectors 

Take an ordered set of three arbitrary, normalized states of a two-level system. Eval¬ 
uate the corresponding discrete Berry phase. Each state is represented by a vector 
on the Bloch sphere. Show analytically that if two of the vectors coincide, then the 
discrete Berry phase vanishes. 

2.2. Two-level system and the Berry connection 

Consider the two-level system defined in Eq. ( |2.61| l, and describe the excited en¬ 
ergy eigenstates using the gauge l+d)^ defined in Eq. ( |2.65| l. Using this gauge, 
evaluate and visualize the corresponding Berry connection vector field A(d). Is it 
well-defined in every point of the parameter space? Complete the same tasks using 
the gauge |-|-d)Ar defined in Eq. p.66| l. 

2.3. Massive Dirac Hamiltonian 

Consider the two-dimensional massive Dirac Hamiltonian H{kjc, ky) = ma^ -f -f 
kyOy, where m C M is a constant and the parameter space is 9 {kx,ky). (a) Take 
a circular loop with radius K in the parameter space, centered around the origin. 
Calculate the Berry phase associated to this loop and the ground-state manifold of 
the Hamiltonian; 7_ {m, k) =7. (b) Calculate the Berry connection B {kx, ky) for the 
ground-state manifold, (c) Integrate the Berry connection for the whole parameter 
space. How does the result depend on ml 


2.4. Absence of a continuous glohal gauge 

In Sect. 2.5.1 we have shown example gauges for the two-level system that were 
not globally smooth on the parameter space. Prove that such globally smooth gauge 
does not exist. 


2.5. Chern number of two-band models 

Consider a two-band lattice model with the Hamiltonian H{k) = d(k) • a. Express 
the Chern number of the lower-energy band in terms of d(k)/|d(k)|. 









Chapter 3 

Polarization and Berry phase 


The bulk polarization of a band insulator is a tricky concept. Polarization of a neu¬ 
tral molecule is easily defined using the difference in centers of the negative and 
positive charges constituting the system. When we try to apply this simple concept 
to the periodic bulk of a band insulator (assuming for simplicity that the positive 
atom cores are immobile and localized), we meet complications. The center of the 
negative charges should be calculated from electrons in the fully occupied valence 
bands. However, all energy eigenstates in the valence band are delocalized over the 
bulk, and so the center of charge of each electron in such a state is ill defined. Nev¬ 
ertheless, insulators are polarizable, and respond to an external electric field by a 
rearrangement of charges, which corresponds to a (tiny) current in the bulk. Thus, 
there should be a way to define a bulk polarization. 

In this chapter we show how a bulk polarization can be defined for band insula¬ 
tors using the so-called modern theory of polarization ll27l 1^ 1^ . The contribution 
of the electrons to the polarization is a property of the many-body electron state, a 
Slater determinant of the energy eigenstates from the fully occupied valence bands. 
The central idea is to rewrite the same Slater determinant using a different orthonor¬ 
mal basis, one composed of localized states, the so-called Wannier states. The con¬ 
tribution of each electron in a Wannier state to the center of charge can then be easily 
assessed, and then added up. 

We discuss the simplest interesting case, that of a one-dimensional two-band 
insulator with one occupied and one empty band, and leave the multiband case for 
later. We show that the center of charge of the Wannier states can be identified with 
the Berry phase of the occupied band over the Brillouin zone, also known as the Zak 
phase 1^ . 

For a more complete and very pedagogical introduction to the Berry phase in 
electron wavefunctions, we refer the reader to a set of lecture notes by Resta ll26l . 
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3 Polarization and BeiTy phase 


The Rice-Mele model 


The toy model we use in this chapter is the Rice-Mele model, which is the SSH 
model of Chapt.[T]with an extra staggered onsite potential. The Hamiltonian for the 
Rice-Mele model on a chain of N unit cells reads 


N N-1 

H = v ^ [\m,B) (m,A| +h.c.) -f w ^ [\m+ 1,A) {m,B\ +h.c.) 

m=\ m=l 


N 

-f M ^ ( |m,A) (m,A| — \m,B) {m,B\ ), 

m=l 


(3.1) 


with the staggered onsite potential u, the intracell hopping amplitude v, and intercell 
hopping amplitude w all assumed to be real. The matrix of the Hamiltonian for the 
Rice-Mele model on a chain of = 4 sites reads 


/mvOOOOOOX 
V—MwOOOOO 
0 w u vOOOO 
00 V— MwOOO 
OOOwMvOO ’ 
00 0 0 V —u V 0 
OOOOOwuv 
\0 0 0 0 0 0 V -uj 


(3.2) 


3.1 Wannier states in the Rice-Mele model 

The bulk energy eigenstates of a band insulator are delocalized over the whole sys¬ 
tem. We use as an example the bulk Hamiltonian of the Rice-Mele model, i.e., the 
model on a ring of N unit cells. As in the case of the SSH model. Sect. |1.2| the 
energy eigenstates are the plane wave Bloch states, 

\'B{k)) = \k)^\u{k)), (3.3) 


with 


\k) = ^ f 


' 2.71 

for k G {5k,25k, ■ ■ ■ ,N5k} with ^k = (3.4) 


We omit the index 1 from the eigenstate for simplicity. The \u{k)) are eigenstates of 
the bulk momentum-space Hamiltonian, 


H{k) = 


u v + we 


-ik 


v + we‘ 


ik 


—u 


(3.5) 
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with eigenvalue E {k). 

The Bloch states \^{k)) are spread over the whole chain. They span the occupied 
subspace, defined by the projector 

L (3.6) 

k£BZ 

The phase of each Bloch eigenstate \^{k)) can be set at will. A change of these 
phases, a gauge transformation, 

|M(^))^e'“«|M(^)); (3.7) 

gives an equally good set of Bloch states, with an arbitrary set of phases a (^) € K 
for k — 5k, 25k, • • • > 2.n. Using this freedom it is in principle possible to ensure that in 
the thermodynamic limit of N ^ the components of l’f'(^)) are smooth, contin¬ 
uous functions of k. However, this gauge might not be easy to obtain by numerical 
methods. We therefore prefer, if possible, to work with gauge-independent quanti¬ 
ties, like the projector to the occupied subspace defined in Eq. (|3.6|l. 


Defining properties of Wannier states 


The Wannier states \w{j)) G J^xtemai 0 =^nternaU with J = 1,... ,N, are defined by 
the following properties: 


(w(/) j w(J)) = 5jij Orthonormal set (3.8a) 

N 


52 |w(y)) (w( 7 )| = P Span the occupied subspace (3.8b) 

7=1 

{m+l\w{j+1)) = {m\w{j)) Related by translation (3.8c) 

lim (w(A/2)| (x —A/2)^ |w(A/2)) < oo Localization (3.8d) 

N^oo 

with the addition in Eq. p.8c| l defined modulo N. Requirement p.8d[ ), that of local¬ 
ization, uses the position operator, 

N 


X = ^ m{\m,A) {m,A\ + \m,B) {m,B\), (3.9) 

m=l 


and refers to a property of |w(y)) in the thermodynamic limit of N ^ that is not 
easy to define precisely. In this one-dimensional case it can be turned into an even 
stricter requirement of exponential localization, (w(J) | m) {m \ w{j)) < 
for some finite localization length ^ G K. 
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Wannier states are inverse Fourier transforms of the Bloch eigenstates 

Because of Bloch’s theorem, all energy eigenstates have a plane wave form not only 
in the canonical basis with p =A,B, but in the Wannier basis as well, 



(3.10) 


with some phase factors a{k). To convince yourself of this, consider the components 
of the right-hand-side in the basis of Bloch eigenstates. The right-hand-side is an 
eigenstate of the lattice translation operator S, with eigenvalue e^‘^, and therefore, 
orthogonal to all of the Bloch eigenstates \^{k’)) with k’ 7 ^ k. It is also orthogonal 
to all positive energy eigenstates, since it is in the occupied subspace. Thus, the only 
state left is |'f'(^)). 

From Eq. p.l0[ ), an inverse Fourier transformation gives us a practical Ansatz 
for Wannier states. 


1 


(3.11) 


There is still a large amount of freedom left by this form, since the gauge function 
a{k) is unconstrained. This freedom can be used to construct Wannier states as lo¬ 
calized as possible. If, e.g., a smooth gauge is found, where in the N ^ limit, 
the components of are analytic functions of k, we have exponential lo¬ 

calization of the Wannier states due to properties of the Fourier transform. (More 
generally, if a discontinuity appears first in the 1 th derivative of \^{k)), the compo¬ 
nents of the Wannier state |w(y)) will decay as .) 

Wannier centers can be identified with the Berry phase 

We first assume that we have found a continuous gauge. The center of the Wannier 
state |w( 0 )) can be calculated, using 



We find that the center of the Wannier state |w(y)) is 
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{w{j)\x\wij)) = ^ [ dk{u{k)\di,u{k))+j. (3.13) 

ZTT J —71 

The second term in this equation shows that the centers of the Wannier states are 
equally spaced, at a distance of one unit cell from each other. The first term, which 
is the Berry phase (divided by 2n) of the occupied band across the Brillouin zone, 
cf. Eq. ( |2.53| ), corresponds to a uniform displacement of each Wannier state by the 
same amount. 

We define the bulk electric polarization to be the Berry phase of the occupied 
band across the Brillouin zone, the first term in Eq. ( |3.13| l, 

^’electric = / dk {u{k) \ di,u{k)). (3.14) 

Although the way we derived this above is intuitive, it remains to be shown that this 
is a consistent definition. Erom Chapt.|^ it is clear that a gauge transformation can 
only change the bulk electric polarization by an integer. We will show explicitly in 
Chapt. I^that the change of this polarization in a quasi-adiabatic process correctly 
reproduces the bulk current. 


3.1.1 Wannier states using the projected position operator 

A numerically stable, gauge invariant way to find a tightly localized set of Wannier 
states is using the unitary position operator ESl . 

X = e'^>‘\ (3.15) 

This operator is useful, because it fully respects the periodic boundary conditions 
of the ring. The eigensystem of X consists of eigenstates localized in cell m with 
eigenvalue Thus, we can associate the expectation value of the position in 

state I'f') with the phase of the expectation value of X, 

{x) = ^lm\og{'¥\X\'¥). (3.16) 

2n 

The real part of the logarithm carries information about the degree of localization ESl 

m. 

In order to obtain the Wannier states, we restrict the unitary position operator to 
the filled bands, defining 


Xp = PxP. (3.17) 

We will show below that in the thermodynamic limit of A —> oo, the eigenstates of 
the projected position operator Xp form Wannier states. 

To simplify the operator Xp, consider 
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where = \\fk'=k+5k, and 0 otherwise. Using this, we have 


Xp^Y.\^{k')){W{k')\X\W{k)){W{k)\ 


k'k 


Y,{u{k + dk)\u{k))-\W{k + 5k)) {^{k) I. (3.19) 


k 


We can find the eigenvalues of Xp using a direct consequence of Eq. p.l9[ ), 
namely, that raising Xp to the Mh power gives an operator proportional to the unity 
in the occupied subspace. 



(3.20) 


We will refer to the constant of proportionality, W S C, given by 

W = {u{2n) I u{2k - 4)) •.. ■ • (m(24) I u{5k)) {u{5k) \ u{2n)) , (3.21) 

as the Wilson loop. Note that W is very similar to a discrete Berry phase, apart from 
the fact that |W| < 1 (although limA^^«, |W| = 1). The spectrum of eigenvalues of Xp 
is therefore composed of the Mh roots of W, 




These eigenvalues have the same magnitude |A„| = \/|W) < 1, and phases in the 
interval [0,27r), spaced by 5^- Because {w{j)\Xp\w{j)) = {w{j)\X\w{j)), the mag¬ 
nitude tells us about the localization properties of the Wannier states, and the phases 
can be interpreted as position expectation values. 

We now check whether eigenstates of Xp fulfil the properties required of Wannier 
states, Eq. ( |3.8| l. The relation {Xp)^ = WP above shows that eigenvectors ofXp span 
the occupied subspace. The eigenstates are related by translation, since 


S'^XpS = e‘^^Xp- 
Xp\>P) = 

XpS\'P) = 


(3.23) 

(3.24) 

(3.25) 


There is a problem with the orthogonality of the eigenstates though. We leave the 
proof of localization as an exercise for the reader. 

The projected unitary position operator Xp is a normal operator only in the ther¬ 
modynamic limit of N Eor finite N, it is not normal, i.e., it does not commute 
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with its adjoint, and as a result, its eigenstates do not form an orthonormal basis. 
This can be seen as a discretization error. 


3.2 Inversion symmetry and polarization 

For single-component, continuous-variable wavefunctions ^{r), inversion about the 
origin (also known as parity) has the effect 'f'(r) —> n^{r) — •?(—r). Two important 
properties of the unitary operator ft representing inversion follow: 77^ = 1, and 
He'*'' = A Hamiltonian has inversion symmetry if IJHW = H. 

When generalizing the inversion operator to lattice models of solid state physics 
with internal degrees of freedom, we have to keep two things in mind. 

First, in a finite sample, the edges are bound to break inversion symmetry about 
the origin (except for very fine-tuned sample preparation). We therefore only care 
about inversion symmetry in the bulk, and require that it take \k) —)■ \ —k). 

Second, each unit cell of the lattice models we consider also has its internal 
Hilbert space, which can be affected by inversion. This includes spin components 
(untouched by inversion) and orbital type variables (affected by inversion) as well. 
In general, we represent the action of inversion on the internal Hilbert space by a 
unitary operator n independent of the unit cell. 

The inversion operator is represented on the bulk Hamiltonian of a lattice model 
by an operator 77, which acts on ,^ntemai as 7t, 

771^) (g) |m) = |—^) (g) TT |m) ; (3.26) 

ft = = lintemal- (3.27) 

The action of the inversion operator on the bulk momentum-space Hamiltonian can 
be read off using its definition, 

nH{k)n-^ = n {k\H^,uik\k) ft-^ = {-klnH^kfi^^ \-k) = ftH{-k)ft\ (3.28) 

A lattice model has inversion symmetry in the bulk, if there exists a unitary and 
Hermitian ft acting on the internal space, such that 

ftH{-k)ft = H{k). (3.29) 

If all occupied bands can be adiabatically separated in energy, so we can focus on 
one band, with wavefunction \u{k)), inversion symmetry has a simple consequence. 
The eigenstates at —k and k are related by 

H{k) \u{k)) = E{k)\u{k)) H{-k)ft\u{k)) = E{k)ft\u{k )); 

^ \u{-k))=e‘‘^^’^'^ft\u{k)). 


(3.30) 

(3.31) 
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For the wavenumbers k = Q and k = K, the so-called time-reversal invariant mo¬ 
menta, this says that they have states with a definite parity, 

l«(0)) =Fo|m(0)); \u{n)) = pTt\u{n)), (3.32) 

with po = ±l; P;r = ±l- (3.33) 


3.2.1 Quantization of the Wilson loop due to inversion symmetry 

We now rewrite the Wilson loop W of a band of an inversion-symmetric one¬ 
dimensional Hamiltonian, assuming we have a discretization into a number 2M of 
A:-states, labeled by j = —M as 

I f|M(27t + ;4)) if;<0 ; 

' ||M(y5;t))) otherwise 

We use Eq. ( |3.31[ ), which takes the form 

\u^j)=e“!’i7t\uj). 

The Wilson loop W of a band of an inversion symmetric one-dimensional insu¬ 
lator can only take on the values ±1. We show this, using M = 3 as an example. 


(3.34) 


(3.35) 


W = {um I M 2 ) (m2 I Ml) (mi I mo) (mo I M_i) (m_i I M_2) (m-2 I Mm) 

= (mm I M2) (M 2 I Ml) (mi I Mo) (Mo|e''^'7t:|Mi) 

(mi| jfe^"^‘e''^7t:|M2) (M 2 I 

= (mi I mo) (mqI |mi) (M 2 I it |mm) (mm | M 2 ) = popn (3.36) 

W = {um I mm-i) ... (mi I Mo) (mo I M_i) ... (m_m+i | um) 

= {um I mm-i) ... (mi I Mo) (mo| e'‘^'it\ui) 

(mi 1|m2) {u2\7te-“^e'‘^^7t\u3)... (mm+i | |mm) 

= (mi I Mo) (moI a |mi) (mm-i| ^ |mm) (mm I mm-i) = ±1 (3.37) 

The statement about the Wilson loop can be translated to the bulk polarization, 
using Eq. ( |3.14[ ). Each band of an inversion symmetric one-dimensional insulator 
contributes to the bulk polarization 0 or 1/2. 
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Problems 


3.1. Inversion symmetry of the SSH model 

Does the SSH model have inversion symmetry? If it has, then provide the corre¬ 
sponding local operator acting in the internal Hilbert space. 

3.2. Ei gensta tes of the projected position operator are localized 

it is claimed that the eigenstates of the projected position operator Xp 


3.1.1 


In Sect. 

form a Wannier set. One necessary condition for that statement to be true is that the 
eigenstates are localized, see Eq. (|3.8d|l. Prove this. 






Chapter 4 

Adiabatic charge pumping, Rice-Mele model 


We now apply the Berry phase and the Chern number to show that the periodically 
and slowly changing the parameters of a one-dimensional solid, it is possible to 
pump particles in it. The number of particles (charge) pumped is an integer per cy¬ 
cle, that is given by a Chern number. Along the way we will introduce important 
concepts of edge state branches of the dispersion relation, and bulk-boundary cor¬ 
respondence. Since we are working towards understanding time-independent topo¬ 
logical insulators, this Chapter might seem like a detour. However, bulk-boundary 
correspondence of 2-dimensional Chern insulators, at the heart of the theory of topo¬ 
logical insulators, is best understood via a mapping to an adiabatic charge pump. 
The concrete system we use in this Chapter is the simplest adiabatic charge pump, 
the time-dependent version of the Rice-Mele (RM) model, 

N N-l 

H{t) = v{t) ^ ( \m,B) {m,A\ -fh.c.) +w{t) ^ ( jm-f 1,A) (m,B| -|-h.c.) 

m=l m=l 

N 

+ u{t)Y^ (|m,A)(m,A|-|m,B)(m,B|), (4.1) 

m=\ 

with the staggered onsite potential u, intracell hopping amplitude v, and intercell 
hopping amplitude w all assumed to be real and periodic functions of time t. In this 
Chapter, we are going to see how, by properly choosing the time sequences, we can 
ensure that particles are pumped along the chain. 


4.1 Charge pumping in a control freak way 

The most straightforward way to operate a charge pump in the Rice-Mele model is 
to make sure that the system falls apart at all times to disconnected dimers. This 
will happen if at any time either the intercell hopping amplitude w, or the intracell 
hopping amplitude v vanishes. We can then use the staggered onsite potential to 
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nudge the lower energy eigenstate to the right. If during the whole cycle we keep 
a finite energy difference between the two eigenstates, we can do the cycle slowly 
enough to prevent excitation. 


Adiabatic shifting of charge on a dimer 

As a first step towards the charge pumping protocol, consider a single dimer, i.e., 
N = 1. Using the adiabatic limit introduced in the last chapter, we can shift charge 
from one site to the other. The Hamiltonian reads 

H{t) = u{t)a^ + v{t)ajc, (4.2) 

with no hopping allowed at the beginning and end of the cycle, at f = 0, we have 
M = 1; V = 0; and atf = r,M = —1; v = 0;. 

We initialize the system in the ground state, which at time f = 0 corresponds to 
|A), a particle on site A. Then we switch on the hopping, which allows the particle to 
spill over to site B, and once it has done that, we switch the hopping off. To ensure 
that the particle spills over, we raise the onsite potential at A and lower it at B. A 
practical choice is 

u{t) = cos{7tt/T); v{t) = sinint/T), (4.3) 

whereby the energy gap is at any time 2. According to the adiabatic theorem, if //(f) 
is varied slowly enough, we will have shifted the charge to B at the end of the cycle. 


Putting together the control freak sequence 


Once we know how to shift a particle from \m,A) to \m,B), we can use that to shift 
the particle further from \m,B) to |m+l,A). For simplicity, we take a sequence 
constructed from linear ramps of the amplitudes, using the function / : [0,1) —> M: 


ifx< 1/8 
if l/8<x<3/8 
if3/8 <x< 1/2 
otherwise 

One period of the pump sequence, for 0 < f < T, reads 

M(f)=/(f/r)-/(f/r + i); 
v(f)=2/(f/r+i); 


fix) = 



(4.4) 


(4.5a) 

(4.5b) 

(4.5c) 
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This period, shown in Fig. |4.1| a), is assumed to then be repeated. Note that we 
shifted the beginning time of the sequence: now at times t/T = n the Hamil¬ 
tonian is the trivial SSH model, tjT = n + 1/4, disconnected monomers, at times 
t/T = n -f 1/2, it is the nontrivial SSH model. 

The time-dependent bulk momentum-space Hamiltonian reads 

H{k,t) = d{k,t)d = (v(f) + w{t) cos k)dx + w{t) sinkdy + u{t)d^, (4.6) 

which can be represented graphically as the path of the vector d{k,t) as the quasi¬ 
momentum goes through the Brillouin zone, k : 0 —> 2n, for various hxed values of 
time f, as in Fig.|4.1|(b). 



time t/T cell index m 

Fig. 4.1 The control freak pump sequence in the Rice-Mele model. The sequence is defined via 
Eqs. H-St and (a) Time dependence of the hopping amplitudes v, w and the sublattice po¬ 
tential u. (b) The surface formed by the vector d(^,f) corresponding to the bulk momentum-space 
Hamiltonian. The topology of the surface is a torus, but its parts corresponding to f 6 [0,0.25]r 
and t 6 [0.75, l]r are infinitely thin and appear as a line due to the vanishing value of w in these 
time intervals, (c) Instantaneous spectrum of the Hamiltonian H(l) of an open chain of W = 10 
sites. Red (blue) points represent states that are localized in the rightmost (leftmost) unit cells and 
have energies between -1 and 1. 
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Visualizing the motion of energy eigenstates 


We can visualize the effects of the control freak pumps sequence in the Rice-Mele 
model by tracing the trajectories of the energy eigenstates. At any time t, each in¬ 
stantaneous energy eigenstate can be chosen confined to a single dimer: either on 
a single unit cell, or shared between two cells. In both cases, we can associate a 
position with the energy eigenstates: the expectation value of the position operator 
X defined as per Eq. 


The trajectories of energy eigenstates in the position-energy space. Fig. 4.2 show 


that the charge pump sequence works rather like a conveyor belt for the eigenstates. 
We engineered the sequence as a unitary operation that pushes all negative energy 
states in the bulk to the right at the rate of one unit cell per cycle (a current of 
one particle per cycle). These orthogonal states, one by one, are pushed into the 
right end region, which has only room for one energy eigenstate. Eigenstates cannot 
pile up in the right end region: if they did, this would violate unitarity of the time 
evolution operator U{t) = where T stands for time ordering, since 

initially orthogonal states would acquire finite overlap. So, states pumped to the 
right edge have to go somewhere, and the only direction they can go is back towards 
the bulk. This on the other hand is only possible, if they acquire enough energy to be 
in the upper band, since all states in the lower band in the bulk are pushed towards 
the right. Moreover, in order to carry these states away from the right edge, and 
make room for those coming from the bulk, the pump sequence has to push upper 
band bulk states towards the left. 

To summarize, the control freak pump sequence is characterized by three state¬ 
ments. The protocol 

• in the bulk, pushes all £ < 0 eigenstates rightwards, by 1 unit cell per cycle, 

• at the right end, pushes 1 eigenstate per cycle with £ < 0 to £ > 0, 

• in the bulk, pushes all £ > 0 eigenstates leftwards, by 1 unit cell per cycle. 

If any one of these statements holds, the other two must also hold as a consequence. 



Fig. 4.2 An instantaneous energy eigenstate of the control freak pump sequence, as it is pumped 
through the system. At negative (positive) energy, it is pumped towards the right (left), upon reach¬ 
ing the right (left) end, it is pumped into the upper (lower) band. . 
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Edge states in the instantaneous spectrum 


We can see the charge pump at work indirectly - via its effect on the edge states 
- using the instantaneous spectrum, the eigenvalues of H{t) of the open chain. An 
example is shown for the control freak pump sequence of the Rice-Mele model on 
a chain of 10 unit cells (20 sites) in Fig. |4.1| Due to the special choice of the control 
freak sequence, the bulk consists of A — 1-fold degenerate states (the bands are flat). 
More importantly, there is an energy gap separating the bands, which is open around 
£ = 0 at all times. However, there are branches of the spectrum crossing this energy 
gap, which must represent edge states. 

To assign “right” or “left” labels to edge states in the instantaneous spectrum, it is 
necessary to examine the corresponding wavefunctions. In case of the control freak 
pump sequence, right (left) edge state wavefunctions are localized on the m = N 
{ m= \) unit cells, and the corresponding energy values are highlighted in green 


(red). The edge state branches in the dispersion in Fig. 4.1 clearly show that 1 state 
per cycle is pushed up in energy at the right edge. 


4.2 Moving away from the control freak limit 

We will now argue that the number of particles pumped by a cycle of a periodic 
adiabatic modulation of an insulating chain is an integer, even if the control freak 
attitude is relaxed. In the generic case, the energy eigenstates are delocalized over 
the whole bulk, and so we will need new tools to keep track the charge pumping 
process. The robust quantization of charge pumping was shown by Thouless, who 
calculated the bulk current directly; we defer this calculation to the next Chapter, 
and here argue using adiabatic deformations. 

As an example for a generic periodically modulated insulator, we take the Rice- 
Mele model, but we relax the control freak attitude. We consider a smooth modula¬ 
tion sequence. 


u(t) = sinf2f; 

(4.7a) 

v(f) = v-|-cosf2f; 

(4.7b) 

w{t) = 1, 

(4.7c) 


where the sequence is fixed by choosing the average value of the intracell hopping, 
V. With V = 1, this sequence can be obtained by an adiabatic deformation of the 
control freak sequence. We show the smooth pump sequence and its representation 
in the d space for v = 1 in Fig. |4.3| 
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time t/T 


cell index m 
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Fig. 4.3 The smooth pump sequence of the Rice-Mele model for v = 1. The hopping amplitudes 
and the sublattice potential (a) are varied smoothly as a function of time. The vector d(^, f) cor¬ 
responding to the bulk momentum-space Hamiltonian (b) traces out a torus in the 3-dimensional 
space. Instantaneous spectrum of the Hamiltonian H{t) on an open chain of A = 10 sites (c) re¬ 
veals that during a cycle, one state crosses over to the upper band on the right edge, and one to the 
lower band on the left edge (dark red/light blue highlights energies of edge states, whose wave- 
functions have than 60% weight on the rightmost/leftmost 2 unit cells). The wavefunctions of the 
edge states (d,e) are exponentially localized to one edge and have support overwhelmingly on one 
sublattice each. In contrast a typical bulk state (f) has a delocalized wavefunction with support on 
both sublattices. . 


Edge states in the instantaneous spectrum 


Consider the spectrum of the instantaneous energies on an open chain, with an ex¬ 
ample for N = 20 unit cells shown in Fig. |4.3| Since this charge sequence was ob¬ 
tained by adiabatic deformation of the control freak sequence above, each branch in 


the dispersion relation is deformed continuously from a branch in Fig. 4.1 


The edge states are no longer confined to a single unit cell, as in the control 
freak case. However, as long as their energy lies deep in the bulk band gap, they 
have wavefunctions that decay exponentially towards the bulk, and so they can be 
unambiguously assigned to the left or the right end. (In case of a degeneracy between 
edge states at the right and left end, we might find a wavefunction with components 
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on both ends. In that case, however, restriction of that state to the left/right end 
results in two seperate eigenstates, to a precision that is exponentially high in the 
bulk length). In Fig. we used the same simple criterion as in Chapt.[T]to define edge 
states: 


N 


I"F) is on the right edge ^ (\{^\m,A)\^+ \{^\m,B)\^j > 0.6', (4.8) 


m=N—l 



(4.9) 


m=l 


As in the control freak case, there is a branch of energy eigenstates crossing over 
from E < 0 to £ > 0 at the right edge, and from £ > 0 to £ < 0, at the left. 

We can define the edge spectrum to consist of edge state branches of the disper¬ 
sion relation, that are clearly assigned to the right end. More precisely, we take 
two limiting energies, e_ and £+, deep in the bulk gap, and only consider en¬ 
ergy eigenstates of the open chain with eigenvalues £„(f) between these limits, 
e_ < £„(f) < £+, with eigenstates localized at the right edge. Each edge state branch 
can begin (1) at f = 0 , as a continuation of another (or the same) edge state branch 
ending att = T, or (2) at £ = £_, or (3) at £ = £+. Each edge state branch can end (1) 
alt — T, to then continue in another (or the same) edge state branch at f = 0, or (2) at 
£ = £_, or (3) at £ = £+. That is 9 possible types of edge state branches. Taking into 
account that the edge state spectmm, like the total spectrum, has to be periodic in f, 
the number of edge state branches entering the energy range £_<£<£+ during a 
cycle is equal to the number of branches leaving it. 

The net number of edge states pumped in energy is a topological invariant 

We now define an integer Q, that counts the number of edge states pumped up in 
energy across at the right edge. Although this quantity is not easily represented by 
a closed formula, it is straightforward to read it off from the dispersion relation 
of an open system. We restrict our attention to the neighbourhood of an energy 
£ deep in the bulk gap around E — 0, such that at all points where E^ = e, the 
derivative dE„/dt does not vanish. Then every edge state energy branch entering 
this neighborhood crosses £ = £ either towards £ > £ or towards £ < £. During 
one cycle, we define for the states at the right edge 

N+ = number of times £ = £ is crossed from £ < £ to £ > £; (4.10) 

N- = number of times £ = £ is crossed from £ > £ to £ < £; (4.11) 

Q = N+ — N- = net number of edge states pumped up in energy . (4.12) 


Note that within the gap, Q is independent of the choice of £. If we found a value 
Qq al E — £o, but a different Q\ f Qq al E = Ei > Eq, this would require a net 
number Qq — Qi of edge state branches at the right edge to enter the energy region 
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Eo < £ < £1 during a cycle but never exit it. Since both £o and £i are deep in the 
bulk gap, away from the bulk bands, this is not possible. 

The net number of edge states pumped up inside the gap on the right edge, Q, is a 
topological invariant: its value cannot change under continuous deformations of the 
Hamiltonian H(t) that preserve the bulk gap. This so-called topological protection 
is straightforward to prove, by considering processes that might change this number. 


We do this using Fig. 4.4 



Fig. 4.4 Adiabatic deformations of dispersion relations of edge states on one edge, in an energy 
window that is deep inside the bulk gap. Starting from a system with 3 copropagating edge states 
(a), an edge state’s dispersion relation can develop a “bump”, (b)-(c). This can change the number 
of edge states at a given energy (intersections of the branches with the horizontal line corresponding 
to the energy), but always by introducing new edge states pairwise, with opposite directions of 
propagation. Thus the signed sum of edge states remains unchanged. Alternatively, two edge states 
can develop a crossing, that because of possible coupling between the edge states turns into an 
avoided crossing (d),(f). This cannot open a gap between branches of the dispersion relation (e), as 
this would mean that the branches become multivalued functions of the wavenumber kr (indicating 
a discontinuity mE{kx), which is not possible for a system with short-range hoppings). Therefore, 
the signed sum of edge states is also unchanged by this process. One might think the signed sum 
of edge states can change if an edge state’s direction of propagation changes under the adiabatic 
deformation, as in (g). However, this is also not possible, as it would also make a branch of the 
dispersion relation multivalued. Deformation of the Hamiltonian can also form a new edge state 
dispersion branch, as in (a)-(h)-(i), but because of periodic boundary conditions along kx, this 
cannot change the signed sum of the number of edge states. 


The number of times edge state branches intersect E = 0 can change because 
new intersection points appear. These can form because an edge state branch is 
deformed, and as a result, it gradually develops a “bump”, local maximum, and the 
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local maximum gets displaced from £ < 0 to £ > 0. For a schematic example, see 
Fig. 4.4 (a)-(c). Alternatively, the dispersion relation branch of the edge state can 


also form a local minimum, gradually displaced from £ > 0 to £ < 0. In both cases, 
the number of intersections of the edge band with the £ = 0 line grows by 2, but the 
two new intersections must have opposite pump directions. Therefore, both N+ and 
N- increase by 1, but their difference, Q = N+ —N-, stays the same. 

New intersection points can also arise because a new edge state branch forms. 
As long as the bulk gap stays open, though, this new edge state band has to be a 
deformed version of one of the bulk bands, as shown in Fig. |4.4| (a)-(h)-(i). Because 
the periodic boundary conditions must hold in the Brillouin zone, the dispersion 
relation of the new edge state has to “come from” a bulk band and go back to the 
same bulk band, or it can be detached from the bulk band, and be entirely inside the 
gap. In both cases, the above argument applies, and it has to intersect the £ = 0 line 
an even number of times, with no change of Q. 

The number of times edge state branches intersect £ = 0 can also decrease if 
two edge state branches develop an energy gap. However, to open an energy gap, 
the edge states have to be pumped in opposite directions. For states pumped in the 
same direction, energy crossing between them can become an avoided crossing, but 
no gap can be opened, as this would violate the single-valuedness of a dispersion 
relation branch, as illustrated in Fig. 4.4 (d)-(f). This same argument shows why 
it is not possible for an edge state to change its direction of propagation under an 
adiabatic deformation without developing a local maximum or minimum (which 
cases we already considered above). As shown in Fig. |4.4| (g), this would entail that 
at some stage during the deformation the edge state branch was not single valued. 


4.3 Tracking the charges with Wannier states 


Electrons in a solid are often described via Bloch states delocalized over the whole 
lattice. As we have seen in Sect. |3.1 [ though, one can represent a certain energy band 
with a set of Wannier states, which inherit the spatial structure (discrete translational 
invariance) of the lattice, and are well localized. Therefore, it seems possible to 
visualize the adiabatic pumping process by following the adiabatic motion of the 
Wannier functions as the parameters of the lattice Hamiltonian are varied in time. 

In fact, we will describe the adiabatic evolution of both the position and energy 
expectation values of the Wannier functions. By this, the toolbox for analyzing the 
adiabatic pumping procedure for control-freak-type pumping is extended to arbi¬ 
trary pumping sequences. 
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4.3.1 Plot the Wannier centers 

According to the result ( |3.13| l, the Wannier center positions of a certain band, in 
units of the lattice constant, are given by the Berry phase of that band divided by 2n. 
Hence, to follow the motion of the Wannier centers during the pumping procedure, 
we calculate the Berry phase of the given band for each moment of time. 

The numerically computed Wannier-center positions obtained for the smoothly 
modulated Rice-Mele sequence [defined via Eqs. ( |4.6| l and ( |4.7| i]. , with v = 1, on 
a finite lattice, are shown in Fig. |4.5^ . Solid (dotted) lines correspond to Wannier 
states of the valence (conduction) band. The results show that during a complete 
cycle, each Wannier center of the valence (conduction) band moves to the right 
(left) with a single lattice constant. Figure [43]7 shows time evolution of the position 
and energy expectation value of a single Wannier center in the valence/conduction 
band, over a few complete cycles. In complete analogy with control-freak pumping, 
these results suggest that the considered pumping sequence operates as a conveyor 
belt: it transports valence-band electrons from left to right, with a speed of one 
lattice constant per cycle, and would transport conduction-band electrons, if they 
were present, from right to left, with the same pace. 

An important result, which is not specific to the considered pumping cycle, arises 
from the above considerations. As pumping is cyclic, the Berry phase of a given 
band at f = 0 is equivalent to that at f = T. Therefore, the displacement of the Wan¬ 
nier center during a complete cycle is an integer. 



position {x) 


position {x) 


Fig. 4.5 Time evolution of Wannier centers and Wannier energies in a smoothly modulated topo¬ 
logical Rice-Mele pumping sequence [defined via Eqs. and (4.51 ]. The parameter of the se¬ 
quence is V = 1, corresponding to a Chem number of 1. In both subfigures, a solid (dotted) line 
corresponds to the valence (conduction) band, (a) Time evolution of the Wannier centers of the 
bands. During a cycle, each Wannier center of the valence (conduction) band moves to the right 
(left) with a single lattice constant, (b) Time evolution of the position and energy expectation value 
of a single Wannier center in the valence/conduction band, over a few complete cycles. 
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4.3.2 Number of pumped particles is the Chern number 

Our above results for the smoothly modulated Rice-Mele cycle suggest the inter¬ 
pretation that during a complete cycle, each electron in the filled valence band is 
displaced to the right by a single lattice constant. From this interpetation, it follows 
that the number particles pumped from left to right, through an arbitrary cross sec¬ 
tion of the lattice, during a complete cycle, is one. Generalizing this consideration 
for arbitary one-dimensional lattice models and pumping cycles, it suggests that the 
number of pumped particles is an integer. 

We now shown that this integer is the Chern number associated to the valence 
band, that is, to the ground-state manifold of the time-dependent bulk momentum- 
space Hamiltonian H{k,t). To prove this is, we first write the Wannier-center dis¬ 
placement Axq t for the complete cycle by splitting up the cycle [0, T] to small 
pieces At: 


n — \ 

Axqj = lim Y Axt-^ti+At, (4.13) 

f=0 

where At = T/n and f, = iAt. Then, we express the infinitesimal displacements with 
the Berry phases, 

i 

Ax{ti,ti + At) = — dk[{u„{ti + At) I di,Un{ti + At)) - {u„{ti) \ dkU„{ti))\^A4-) 
2,71 J — 71 

where the k argument is suppressed for brevity. Since the k = —n and k = n values 
are equivalent, the above integral can be considered as a line integral of the Berry 
connection to the closed boundary line dRi of the infinitesimally narrow rectangle 
{k,t) G Ri = [—TT, tt) X [tiAi + At\, that is. 


Ax(ti,ti + At)^ — i aW-c/R. (4.15) 

271 JdRi 

Using the fact that we can choose a gauge that is locally smooth on that rectangle 
Ri, we obtain 


Ax(ti,ti + At) = — [ B^"'>dkdt, (4.16) 

27t JdRi 

where is the Berry curvature associated to the nth eigenstate manifold of 
Together with Eq. ( |4.13| ), this result ensures that the Wannier-center displacement is 
the Chern number; 


4x0, = — B^"Ukdt. (4.17) 

27t J-K 

Note that even though we have not performed an explicit calculation of the 
valence-band Chern number of the smooth Rice-Mele pump cycle with v = 1, by 
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looking at the motion of the corresponding Wannier centers we can conclude that 
the Chern number is 1. 


4.3.3 Tuning the pump using the average intracell hopping 
amplitude v 


So far, the discussed results were obtained for the v = 1 special case of the smoothly 
modulated Rice-Mele pumping cycle. Now we ask the question: can the number of 
pumped particles be changed by tuning the parameter v? To show that the answer is 


affirmative and the pump has such a tunability, on Fig. 4.6 1 we plot the instantaneous 
energy spectrum corresponding to v = — 1. The spectrum reveals that this sequence, 
similarly to the v = 1 case, does pump a single particle per cycle. However, the 
direction of pumping is opposite in the two cases: Fig. |4.6^ shows that during a 
cycle, one edge state on the left (light blue) crosses over from the valence band to 
the conduction band, revealing that the particles are pumped from right to left in the 
valence band. 



Fig. 4.6 The smooth pump sequence of the Rice-Mele model for v = — 1, revealing reversed 
pumping with respect to the v = 1 case, (a) Instantaneous spectrum of the Hamiltonian H{t) on an 
open chain of A = 10 sites. During a cycle, one edge state on the left (light blue) crosses over from 
the valence to the conduction band, revealing that there is a single particle per cycle is pumped 
from right to left in the valence band. Wavefunctions of the edge states (b,c) as well as a typical 
bulk state (d) are also shown. 


4.3.4 Robustness against disorder 

So far, we have the following picture of an adiabatic pump in a long open chain. If 
we take a cross section at the middle of the chain, a single particle will be pumped 
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through that, from left to right, during a complete cycle. This implies that at the end 
of the cycle, the number of particles on the right side of the cross section has grown 
by one, therefore the final state reached by the electron system is different from the 
original, ground state. This implies that during the course of the cycle, a valence- 
band energy eigenstate deformed into a conduction-band state, and that occured on 
the right edge of the chain; the opposite happens on the left edge. 

Does that qualitative behavior change if we introduce disorder in the edge re¬ 
gions of the open chain? No: as long as the bulk of the chain remains regular, the 
pump works at the middle of the chain, and therefore the above conclusion about the 
exchange of a pair of states between the valence and conduction bands still holds. 
On the other hand, introducing disorder in the bulk seems to complicate the above- 
developed description of pumping in terms of Wannier-center motion, and therefore 
might change number of edge states and the qualitative nature of the instantaneous 
energy spectrum. 


Problem s 

4.1. Further control-freak pump cycles 

Construct a control-freak pump cycle where the spectrum of the bulk does not 
change during the entire cycle, and the pumped charge is (a) zero (b) one. 




Chapter 5 

Current operator and particle pumping 


In the previous chapter, we described quantized adiabatic pumping of particles in a 
one-dimensional lattice in an intuitive and visual fashion, using the concepts of the 
control-freak pumping cycle and the time evolution of the Wannier centers. Here, 
we provide a more formal description of the same effect. For simplicity, we consider 
two-band insulator lattice models with a completely filled lower band, which are 
described by a periodically time-dependent bulk momentum-space Hamiltonian of 
the form 


H{k,t)=A{k,t)-a, (5.1) 

where d{k,t) is a dimensionless three-dimensional vector fulfilling d{k,t) > 1, and 
a is the vector of Pauli matrices. This Hamiltonian is periodic both in momentum 
and in time, H{k + 27Z,t) = H{k,t T) — H[k,t), where T is the period of the time 
dependence of the Hamiltonian. The minimal energy gap between the two eigen¬ 
states of the Hamiltonian is 2. Furthermore, the frequency characterising the peri¬ 
odicity of the Hamiltonian is = 27r/7’. We call the periodically time-dependent 
Hamiltonian quasi-adiabatic, if ^ 1, and the adiabatic limit is defined as > 0, 
that is, T ^ 

For example, d can be chosen as 


( v-|-cosf2f -|-cosk\ 

sink j , (5.2) 

sinf2f j 

corresponding to the smoothly modulated Rice-Mele model, see Eq. ( |4.7| i and Eq. 

( |431 ). 

We will denote the eigenstate of H{k,t) with a lower (higher) energy eigenvalue 
as \ui{k,t)) (\u2{k,t)))- With this notation, we can express the central result of this 
chapter: the momentum- and time-resolved current carried by the electrons of the 
filled band equals the Berry curvature associated to that band. As a consequence, 
the number ^ of particles pumped through an arbitrary cross section of an infinite 
one-dimensional crystal during a complete adiabatic cycle is the momentum- and 
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time integral of the Berry curvature, that is, 

^ =—i— [ dt j dk(dk{u\{k^t)\dtU\{k,t)) — dt {u\{k,t)\dkUi{k,t))). {53) 
2 k Jo Jbz 

This is the Chern number associated to the ground-state manifold of H{k,t). As the 
latter is an integer, the number of pumped particles is quantized. This result was 
discovered by David Thouless lf33l . 

We derive Eq. ( |5.3| l via the following steps. In Sect. |5.1[ we consider a generic 
time-dependent lattice Hamiltonian, we express the number of particles moving 
through a cross section of the lattice using the curi'ent operator and the time-evolving 


energy eigenstates of the lattice. Then, in Sect. 5.2.2 we provide a description 
of the time-evolving energy eigenstates of the lattice in the case of periodic and 
quasi-adiabatic time dependence of the lattice Hamiltonian. This allows us to ex¬ 
press the number of pumped particles for quasi-adiabatic time dependence. Finally, 
in Sect. |5.3[ building on the latter result for quasi-adiabatic pumping, we take the 
adiabatic limit and thereby establish the connection between the current, the Berry 
curvature, the number of pumped particles, and the Chem number. 


5.1 Particle current at a cross section of the lattice 

Our aim here is to express the number of particles pumped through a cross sec¬ 
tion of the lattice, assuming that the time evolution of the Bloch states due to the 
time-dependence of the Hamiltonian is known. As intermediate steps toward this 
end, we derive the real-space current operator and the diagonal matrix elements 
of the momentum-space current operator, and establish an important relation be¬ 
tween those diagonal matrix elements and the momentum-space Hamiltonian. For 
concreteness, we first discuss these using the example of the Rice-Mele model in¬ 
troduced in the preceding chapter. It is straightforward to generalize the results for 
lattice models with an generic internal degree of freedom; the generalized results are 
also given below. Finally, we use the relation between the current and the Hamilto¬ 
nian to express the number of pumped particles with the time-evolving states and 
the Hamiltonian. 


5.1.1 Current operator in the Rice-Mele model 

We consider the Rice-Mele model with A ^ 1 unit cells and periodic boundary 
conditions. The real-space bulk Hamiltonian A^uik has the almost the same form as 
Eq. with the difference that the sum corresponding to intercell hopping runs 
up to N, and in accordance with the periodic boundary condition, the unit cell index 
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m should be understood as (m modA^). The bulk momentum-space Hamiltonian 
H{k) of the model is given in Eq. (|4.6|l. 


(f -I-1) ' 


(« + i) 


N 


segment S 


Fig. 5.1 A segment enclosed by two cross sections in the one-dimensional SSH model. The seg¬ 
ment S is defined as the part of the chain between the pth and {q + l)th unit cells. The current 
operators corresponding to the two cross sections can be established by considering the temporal 
change of the number of particles in the segment. The dashed line represents the periodic boundary 
condition. 


Influx of particles into a segment of the crystal 

We aim at establishing the operator representing the particle current flowing through 
a cross section of the one-dimensional crystal. We take a cross section between the 
pth and (p-f l)th unit cells, and denote the corresponding current operator as jp+ij 2 - 
To find the current operator, we first consider a segment S of the crystal, stretching 
between (and including) the {p + l)th and ^th unit cells, where q> p+\. The 
number of particles in that segment S is represented by the operator 

^s=Y^ Y, \m,a){m,a\. (5.4) 

meS ae{A,B} 

Now, Ehrenfest’s theorem ensures that the time-evolution of the number of particles 
embedded in the segment S follows 

dr{£s)t = -i{[^s,H{t)]), (5.5) 

Hence, we identify the operator describing the influx of particles into the segment 
as 


js{t) = -i[.^s,H{t)]. 

A straigthforward calculation shows that Eq. ( |5.6| l implies 

Js{t) = -iw{t){\p+l,A){p,B\-\p,B) (p+l,A| 
+ {q+\,A\ — \q+\,A) {q,B\) 


(5.6) 


(5.7) 
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Remarkably, the operator js{t) is time dependent, if the hopping amplitude w{t) is 
time dependent. 


Current operator at a cross section of the crystal 

Clearly, the terms in Eq. (|53 can be separated into two groups: the hrst two terms 
are hopping operators bridging the cross section p + 1/2, and the last two terms are 
bridging the cross section q+\/2. Thereby we dehne 

Jm+l/lit) = -iw{t) (|m+ 1,A) {m,B\ - \m,B) (m+ 1,A|), (5.8) 

and use this dehnition to express js as 

is = jp+\/2- jq+\/2- (5-9) 

This relation allows us to interpret /m+i/2(0 the current operator describing par¬ 
ticle flow, from left to right, across the cross section m+l/2. 


Relation of the current operator to the Hamiltonian and to the group velocity 

Later we will need the momentum-diagonal matrix elements of the curi'ent operator, 
which are defined as 


}m+l/2{k,t) = {k\ jm+i/2(f) \k) . (5.10) 

For the Rice-Mele model under consideration, these can be expressed using Eqs. 
(prsl) and (|5:8]l as 

From a comparison of this result and Eq. ( |4.6| l, we see that the momentum-diagonal 
matrix elements of the current operator are related to the momentum-space Hamil¬ 
tonian as 


Jm+i/2{k,t) = ^dkH{k,t). (5.12) 

This is the central result of this section. Even though we have derived it only for the 
case of the Rice-Mele model, it is a generic result. A generalization is outlined in 
the next section. 

In the case of a lattice without an internal degree of freedom, it is easy to see that 
Eq. ( |5.12| i establishes the equivalence /m+i/ 2 (^T) = Vk{t)/N between the current 
and the (instantaneous) group velocity Vk{t)\ in this case, the bulk momentum-space 
Hamiltonian equals the dispersion relation and its momentum deriva- 
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live Vj^it) = is defined as the group velocity of the energy eigenstates. This 

correspondence generalizes to lattices with an internal degree of freedom as well. 
The current carried by an instantaneous energy eigenstate \'y„{k,t)) = |^) 0 
of such a lattice is 


(M«l/m+l/2 \Un) = {Un\ [dkH] \u„) 


= — {u„ 

N ^ 






— ^ I \ ^n) (ttrt I T Efi I + Efi \uf^ I] |m«)} 


l/Tr;. 117:1/ I \1 ^kEfl 

— ^ [OkEn + E„ai^ \Ufi \Un)) — ^ , 


(5.13) 


where the arguments {k,t) are suppressed for brevity. 


5.1.2 Current operator in a generic one-dimensional lattice model 

Here we prove Eq. ( |5.12| l in a more general setting. Previously, we focused on 
the Rice-Mele model that has only two bands and hopping only between nearest- 
neighbor cells. Consider now a general one-dimensional lattice model with Nb bands 
and hnite-range hopping with range 1 < r <C A^; here, r = 1 corresponds to hopping 
between neighbouring unit cells only. As before, we take a long chain with A ^ 1 
unit cells, and assume periodic boundary conditions. 

The real-space Hamiltonian has the form 

N Nt, 

H{t)= Y. L f^ma.m'a'(t)\m,a){m'a'\, (5.14) 

m,w'=l a,(x'=l 

where m and m' are unit cell indices, and a,a' G {1,2,...,Nb} correspond to the 
internal degree of freedom within the unit cell. Due to the hnite-range-hopping as¬ 
sumption, the Hamiltonian can also be written as 

N r Nf, 

H{t)=YY E f^m+i,a-m,a'it)\m + i,a){ma'\. (5.15) 

m=\ i=—r a,a'=l 

Again, a unit cell index m should be understood as (m mod A). Also, due to the 
discrete translational invariance, we have aS implying 

N r Nfy 

^(0 = E E E (0 \m + i, a) {mo! 

m=l i=—ra.a’=\ 


(5.16) 
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The bulk momentum-space Hamiltonian then reads 

H{k,t) = {k\H{t) \k) = ^ l«) (a'| • (5.17) 

m=—r 


Next, we establish the operator representing the particle current flowing through 


a cross section of the one-dimensional crystal, the same way we did in Sect. 5.1.1 
We take a cross section between the pth and {p-\-\ )th unit cells, and denote the cor¬ 
responding current operator as jp+\/ 2 - To find the current operator, we first consider 
a segment S of the crystal, stretching between (and including) the {p -f l)th and ^th 
unit cells, where q — p> r. The number of particles embedded in that long segment 
is represented by the operator 


Nt 

= E L l'”“) 

meSa=\ 


bna \. 


(5.18) 


As discussed in the preceding section, we identify the operator describing the influx 
of particles into the wire segment as 

js{t) = -i[^s.H{t)]. (5.19) 


From this, a straigthforward calculation shows that 

75(0 ^^m€S^m'^S^a,a'=l \in(x} {vtl CC | 

^ ) (^Ct|] . (5.20) 

Note that Eq. ( |5.20[ ) testifies that the operator 75 (f) is constructed only from those 
hopping matrix elements of the Hamiltonian that bridge either the p -f 1/2 or the 
1/2 cross sections of the crystal, i.e., one of the two cross sections that terminate 
the segment under consideration. This is ensured by the condition that the segment 
is at least as long as the range r of hopping. A further consequence of this is that the 
terms in Eq. ( |5.20[ ) can be separated into two groups; one containing the hopping 
matrix elements bridging the cross section p + 1 / 2 , and one with those bridging the 
cross section q+l/2. The former reads 

l^a) (5.21) 

^m'a '^ . 

Using this as a definition for any cross section m +1/2, we conclude that Eq. ( |5.9| l 
holds without any change in this generalized case as well. This conclusion allows 

us to interpret 7 „,+i /2 the current operator describing particle flow, from left to 

right, across the cross section n + 1 / 2 . 

After defining the momentum-diagonal matrix elements of the current operator 
exactly the same way as in Eq. (|5.10|), the relation between the current operator and 
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the Hamiltonian has exactly the same form as in Eq. ( |5.12| l. This can be proven 
straightforwardly using Eqs. (|5.10|l, (|5.21|i, and the A:-derivative of Eq. (|5.17|i. 


5.1.3 Number of pumped particles 


Let us return to particle pumping in insulating two-band models. As the electrons 
are assumed to be non-interacting, and the time-dependent lattice Hamiltonian has 
a discrete translational invariance for all times, the many-electron state ^>(f) is a 
Slater determinant of Bloch-type single-particle states I’Ll (A:,f)) = \k) (g) \ui{k,t)). 
Here, we adopted the notation introduced in Sect. |1.2[ with n = 1 referring to the 


hlled band. There are also two additions with respect to the notation of Sect. 1.2 


hrst, we explicitly denote the time dependence of the state; second, we added a tilde 
here to denote that the state j"?! {k,t)') is not an instantaneous lower-band eigenstate 
\^i{k,t)) of the Hamiltonian, but is slightly different from that due to the quasi- 
adiabatic driving. 

The number of particles pumped through the cross section m-f 1/2 within the 
time interval t G [0, T] is the time-integrated current, that is. 






(5.22) 


where 7^+ 1 / 2(0 is the many-particle generalisation of the current operator dehned 
in Eq. ( |5.21| i, or, for the special case of the Rice-Mele model, in Eq. ( |5.8| l), and 
0{t) is the many-electron Slater determinant formed by the hlled Bloch-type single¬ 
particle states, introduced in the preceding paragraph. Equation ( |5.22| l can be con¬ 
verted to an expression with single-particle states: 


keBZ 


(5.23) 


which is related to the momentum-diagonal matrix elements of the current operator 
as 


^ = / dt Y. («i(00l/m+i/2(00l«i(00)- (5-24) 

keBZ 

Einally, this is rewritten using the current-Hamiltonian relation Eq. ( |5.12| ) as 

L {ai(k,t)\dkH{k,t)\ui{k,t)). (5.25) 

^ ■'O keBZ 

To evaluate this in the case of adiabatic, periodically time-dependent Hamiltonian, 
we hrst need to understand how the two-level wave functions \ui{k,t)) evolve in 
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time in the quasi-adiabatic case; then we can insert those in Eq. ( |5.25| l, and take the 
adiabatic limit. 


5.2 Time evolution governed by a quasi-adiabatic Hamiltonian 

Our goal here is to describe the time evolution of Bloch-type electronic energy 
eigenstates. Nevertheless, as the pumping dynamics preserves the wavenumber k, 
the task simplifies to describe the dynamics of distinct two-level systems, labelled 
by the wavenumber k. Therefore, in this section we discuss the dynamics of a single 
two-level system, hence the wavenumber k does not appear in the formulas. We will 
restore k when evaluating the number of pumped particles in the next section. To de¬ 
scribe the time evolution of the electronic states subject to quasi-adiabatic driving, 
it is convenient to use the so-called parallel-transport gauge or parallel-transport 
time parametrization, which we introduce below. Then, our goal is reached by per¬ 
forming perturbation theory in the small frequency ^ 1 characterizing the quasi- 
adiabatic driving. 


5.2.1 The parallel-transport time parametrization 


As mentioned earlier, the instantaneous energy eigenstates of the bulk momentum- 
space Hamiltonian H{k,t) are denoted as \u„(k,t)). Here, in order to simplify the 
derivations, we will use a special time parametrization (gauge) for these eigenstates, 
which is called the parallel-transport time parametrization or parallel-transport 
gauge. As mentioned above, we suppress the momentum k. 

We will call the smooth time parametrization \u„{t)) of the instantaneous nth 
eigenstate of the Hamiltonian H{t) a parallel-transport time parametrization, if for 
any time point t and any band n, it holds that 

{u„{t)\dt\un{t)) =0. (5.26) 


Using a time parametrization with this property will simplify the upcoming calcu¬ 
lations of this section. 


In Sect. 2.3 we used smooth parametrizations that were defined via a parameter 
space, and therefore were cyclic. For any time parametrization \u'„{t)) having those 
properties, we can construct a parallel-transport time parametrizaton |M„(f)) via the 
definition 


\un{t))=e‘^"^’^\u'„{t)), (5.27) 

where 7 „(f) is the adiabatic phase associated to the adiabatic time evolution of the 
initial state = 0)), governed by our adiabatically varying Hamiltonian H{t): 
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= dt' . (5.28) 

The fact that \un{t)) indeed fulfils Eq. ( |5.26| ) can be checked by performing the time 
derivation and the scalar product on the left hand side of the latter. 

As an interpretation of Eq. ( |5.27| l, we can say that a parallel-transport time 
parametrization is an adiabatically time-evolving state without the dynamical phase 
factor. Furthermore, as the Berry phase factor e'7n(^) is, in general, different from 1, 
the parallel-transport time parametrization (|5.27|l is, in general, not cyclic. 


5.2.2 Quasi-adiabatic evolution 

Here, following Thouless l|33l, we describe the quasi-adiabatic time evolution using 
stationary states, also known as Floquet states, that are characteristic of periodi¬ 
cally driven quantum systems. Again, we focus on two-level systems as introduced 
in Eq. (|53, and suppress the wave number k in our notation. The central result 
of this section is Eq. ( |5.42| l, which expresses how the instantaneous ground state 
mixes weakly with the instantaneous excited state due to the quasi-adiabatic time 
dependence of the Hamiltonian. In the next section, this result is used to evaluate 
the particle current and the number of pumped particles. 

As an example, we can consider the state corresponding to the wavenumber k = 0 
in the smoothly modulated Rice-Mele model with v = 1, see Eq. ( |5.2| i: 

( 2 -|-cosf2f \ 

0 . (5.29) 

sin(f2f) j 


Stationary states of periodically driven dynamics 

The stationary states are special solutions of the periodically time-dependent Schrodinger 
equation, which are essentially periodic with period T ; that is, which fulfill | (//(f + T)) = 
e^“l’ I V7(f)) for any t, with (j) being a f-independent real number. The number of such 
nonequivalent solutions equals the dimension of the Hilbert space of the quantum 
system, i.e., there are two of them for the case we consider. Here we describe station¬ 
ary states in the quasi-adiabatic case, when the time evolution of the Hamiltonian 
is slow compared to the energy gap between the instantaneous energy eigenvalues: 

^ 1. This condition suggest that the deviation from the adiabatic dynamics is 
small, and therefore each stationary state is in the close vicinity of either the instan¬ 
taneous ground state or the instantaneous excited state. Thereby, we will label the 
stationary states with the band index n, and denote them as \u„{t)). 

Since, after all, we wish to describe pumping in a lattice with a filled lower band 
and an empty upper band, we mostly care about the stationary state corresponding to 
the lower band, |mi (f)), and therefore want to solve the time-dependent Schrodinger 
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equation 


— idt\ui{t)) + H{t) |Mi(f)) =0. 


(5.30) 


Making use of the parallel-transport gauge 


We characterize the time evolution of the wave function |i5i (f)) by a time-dependent 
linear combination of the instantaneous energy eigenstates; 

|mi( f)) = ai(f)) + \u2{t )), (5.31) 


Recall that we are using the parallel-transport time parametrization, having the prop¬ 
erties ( |5.26| ) and ( |5.27| ). Therefore, in the adiabatic limit 0, we already now that 
fli(f) = 1 and 02 ( 1 ) = 0. Here, we are mostly interested in the quasi-adiabatic case 
defined via ^ 1, and then it is expected that ai (f) ^ 1 and 02 ( 1 ) ^ I. 

Before making use of that consideration in the form of perturbation theory in 
Q, we convert the time-dependent Schrodinger equation ( |5.30| l to two differential 
equations for the two unknown functions ai(f) and a 2 (f). We insert |Mi(f)) of Eq. 
(|5.31|l to the time-dependent Schrodinger equation (|5.30|l, yielding 


-id, \ui{t'))+a2it)e-'^odt'E2{t') |^2(f'))' 

+Hit) \u,{t))+a2{t)e-‘^odt'E2U') |„2(f))' 


(5.32) 


= 0 


After evaluating the time derivatives, the left hand side consists of 8 terms. Using the 
instantaneous eigenvalue relations H{t) ]««(?)) = E„{t) |M„(f)), two pairs of terms 
annihilate each other, and only 4 terms remain; 


dl(t) \u,{t))+a,{t)d, |«i(f)) |m2 (0) 

+ a2{t)e-'^^‘‘''’^^^''>d,\u2{t)) =Q, (5.33) 


where £(?) =E 2 {t)—E\{t) =2d{t). 

Projecting Eq. ( |5.33| l onto (mi (f)| and (M 2 (f)|, respectively, and making use of the 
parallel-transport gauge, yields 

diit)+a2{t)e-‘fo‘‘''’^^^'Huiit)\d,\u2{t)) =0, (5.34) 

a^{t){u2{t)\d,\uiit))+d2{t)e-‘^odr'E(,') (5 35 ) 

The latter result can be rewritten as 

d 2 {t) = -ai{t) {U 2 {t)\d, 


(5.36) 
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Making use of the quasi-adiabatic condition 

Note that the quasi-adiabatic condition has not been invoked so far; this is the next 
step. As mentioned above, the quasi-adiabatic nature of the Hamiltonian suggests 
that one of the two stationary states will be close to the instantaneous ground state, 
suggesting ai{t) ~ 1 and a 2 {t) ^ Q. Furthermore, we know that (mi( f)| dt \u 2 {t)) ^ 
£2, since variations in \u„{t)) become slower as the adiabatic limit is approached. 
The latter relation is explicitly demonstrated by the example in Eq. ( |5.29| l: if we use 
l“i(0) = (—sin(0/2),cos(0/2))^ and \u 2 it)) = (cos(0/2),sin(0/2)) with 0 = 
arctanfulfilling the parallel-gauge criterion, we find {ui{t)\dt \u 2 {t)) = 

l+2cosnf 

10+8cosl2?- 

As we are interested in the quasi-adiabatic case ^ 1, we drop those terms from 
( |5.34| l and ( |5.36| l that are at least second order in £2. This results in 

di(f)=0, (5.37) 

d2{t) = (5.38) 

If we assume ai(f = 0) = 1 +o{£2), then the first equation guarantees that a\{t) = 
1 -|-o(f2). Then this allows for a further simplification of Eq. ( |5.38| ): 

d2{t) = -{u2{t)\dMt))e'^l^‘‘^'^^*'l (5.39) 


Solution of the equation of motion 


The remaining task is to solve Eq. ( |5.39| l for 02 (0- Instead of doing this in a con¬ 
structive fashion, we give the solution a 2 {t) and prove that it indeed fulfills Eq. 
( |5.39| l up to the desired order. The solution reads 

Et 

Eirst, let us check if it solves the differential equation (|5.39|i: 


^ . {d, {U2{t)\ d, \ui{t))) j ■ (dtE,) {U2it)\ d, \ui{t)) J 

' E, Ef 




(5.41) 


The first two terms on the right hand side scale as £2^, whereas the third one scales 
as £2. Hence we conclude that in the quasi-adiabatic case, Eq. ( |5.40| l is the solu¬ 
tion of Eq. ( |5.39| l we were after. The corresponding solution of the time-dependent 
Schrodinger equation ( |5.30| ) is constructed using Eqs. ( |5.31| ), a\{t) = 1 and ( |5.38| ), 
and reads 






(5.42) 
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In words, Eq. ( |5.42| l assures that the stationary state has most of its weight in the 
instantaneous ground state |Mi(f)), with a small, ~ <C 1 admixture of the instan¬ 
taneous excited state \u 2 {t)). Interestingly, even though this small admixture van¬ 
ishes in the adiabatic limit —> 0, the corresponding contribution to the number of 
pumped particles can give a finite contribution, as the cycle period T goes to infinity 
in the adiabatic limit. This will be shown explicitly in the next section. 

Finally, we show that this state |Mi(f)) is indeed stationary. That is proven if we 
can prove that |mi( 7’)) is equal to |mi( 0)) up to a phase factor. This arises as the 
consequence of the following fact. If the Berry phase associated to the state \u\) is 
7 , that is, if |mi( 7’)) = e'^|Mi(0)), then 


[dt\ui{t))]T = liin 

£->■0 


|Mi(7’-fe))-|Mi(7’)) 


^.^ e'y\u,{e))-e'y\um) 

£^o e 


= e'ydt\um)- 


(5.43) 


Therefore, the two terms in the square bracket of Eq. ( |5.42| ) acquire the same phase 
factor e'y at the end of the cycle, hence the obtained \u\{t)) solution is stationary. 


5.3 The pumped current is the Berry curvature 

The number of particles pumped through an arbitrary cross section of the one¬ 
dimensional lattice, in the duration T of a quasi-adiabatic cycle, is evaluated com¬ 
bining Eqs. ( |5.25| l and ( |5.42| i. We define the momentum- and time-resolved current 
of the filled band as 


4 +i/ 2(^’0 = ^ {u\{k,t)\dkH{k,t)\ui{kf )), 


N 


(5.44) 


and perform the usual substitution ^ LjteBZ • • • = /bz S ■ ■ ■ ’ yielding the following 
formula for the number of pumped particles: 


J2= dt 

Jo Jbz 2k 




(5.45) 


In the rest of this section, we show that the momentum- and time-resolved current 
is the Berry curvature associated to the filled band, and therefore the number of 
pumped particles is the Chern number, which in turn is indeed an integer. 

To this end, we insert the result ( |5.42| l to the definition ( |5.44| l. The contribution 
that incorporates two lower-band wave functions \ui{k,t)), is finite; however, its 
integral over the BZ vanishes, and therefore we disregard it as it does not contribute 
to particle pumping. Hence the leading relevant contribution is the one incorporating 
one filled-band \ui{k,t)) and one empty-band \u 2 {k,t)) wavefunction: 


.•( 1 ) r,, _A‘^i\[dkH]\u2){u2\d,\ui) 

Jm+l/iy'^d) — I p~ ■■ +C.C. 


E 


(5.46) 
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where the k and t arguments are suppressed for brevity. 
Now we use 


{u\\[di,H]\u2) = (£1 -£2) {dku\\u2) = -E (5iMi|M2), (5.47) 

which has a straightforward proof using the spectral decomposition i7 = £i|Mi)(Mi| + 
£2 \u 2 ) (M2I of the Hamiltonian and the fact that (9^ {u\\ uj) = 0. Therefore, 

7m+i/2 = I M2) {U2\dt |mi) +C.C.. (5.48) 

Since we use the parallel-transport gauge, we can replace the projector | M2 ) (m 2 | with 
unity in the preceding formula, hence the latter can be simplified as 

jm+l/2 = -i{dkui\dtui) +c.c. = -i{{di,ui\d,ui) - {dtui\dkui)) 

= -iidk{ui\dtui) - dt {ui\dkui)). (5.49) 

This testifies that the momentum- and time-resolved current is indeed the Berry 
curvature corresponding to the filled band, and thereby confirms the result promised 
in Eq. ( |5.31 l. 



As a straightforward application of our result, we calculate the time dependence 
of the current and the number of pumped particles through an arbitrary cross section 
in the smoothly modulated Rice-Mele model, see Eq. The results correspond¬ 
ing to four different values of the parameter v are shown in Fig. 5.2 The momentum- 
and time-resolved current filled band can be obtained analytically 

from Eq. ( |5.49| l. Then, the time-resolved current j is defined as the integrand of the t 
integral in Eq. ( |5.45| ). We evaluate j via a numerical k integration, and plot the result 
in Fig. The number of pumped particles as a function of time is then eval¬ 

uated numerically via ^{t) = fQdtj{t)\ the results are shown in Fig. ( |5.2| )b. These 
results confirm that the number of particles pumped through the cross section during 
the complete cycle is an integer, and is given by the Chern number associated to the 
pumping cycle. 
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In this chapter, we have provided a formal description of adiabatic pumping in 
one-dimensional lattices. After identifying the current operator describing particle 
flow at a cross section of the lattice, we discussed the quasi-adiabatic time evolution 
of the lower-band states in a two-band model, and combined these results to express 
the number of pumped particles in the limit of adiabatic pumping. The central result 
is that the momentum- and time-resolved current carried by the lower-band electrons 
is the Berry curvature associated to their band. 


Problems 

5.1. The smooth pump sequence of the Rice-Mele model. 

For the smoothly modulated Rice-Mele pumping cycle, see evaluate the 

momentum- and time-dependent current density, and the number of particles pumped 
through an arbitrary unit cell boundary as the function of time; that is, reproduce Fig. 

EH 


5.2. Parallel-transport time parametrization. 

Specify a parallel-transport time parametrization for the ground state of the two- 
level Hamiltonian dehned by Eqs. ( |5T] ), (1^, and (a) A: = 0 (b) k = 7t. 


5.3. Quasi-adiabatic dynamics with a different boundary condition. 

In Sect. 5.2.2 we described a stationary state of a quasi-adiabatically driven two- 
level system, and used the result to express the number of particles pumped during 
a complete cycle. How does the derivation and the result change, if we describe the 
dynamics not via the stationary state, but by specifying that the initial state is the 
instantaneous ground state of the Hamiltonian at f = 0? Is the hnal result for the 
number of pumped particles different in this case? 


5.4. Adiabatic pumping in multiband models. 

Generalize the central result of this chapter in the following sense. Consider adi¬ 
abatic charge pumping in a one-dimensional multi-band system (n = 1,2,... ,N), 
where the energies of the hrst N- bands {n — 1,2,.. .A^-) are below the Fermi en¬ 
ergy and the energies of the remaining bands (n = N- + 1,...,N) are above the 
Fermi energy, and the bands do not cross each other. Show that the number of par¬ 
ticles adiabatically pumped through an arbitrary cross section of the crystal is the 
sum of the Chem numbers of the hlled bands. 



Chapter 6 

Two-dimensional Chern insulators - the 
Qi-Wu-Zhang model 


The unique physical feature of topological insulators is the guaranteed existence 
of low-energy states at their boundaries. We have seen an example of this for a 
one-dimensional topological insulator, the SSH model: A finite, open, topologically 
nontrivial SSH chain hosts 0 energy bound states at both ends. The bulk-boundary 
correspondence was the way in which the topological invariant of the bulk - in the 
case of the SSH chain, the winding number of the bulk Hamiltonian - can be used 
to predict the number of edge states. 

We will show that the connection between the Chern number and the number of 
edge-state channels is valid in general for two-dimensional insulators. This is the 
statement of bulk-boundary correspondence for Chern insulators. The way we will 
show this inverts the argument above: taking any two-dimensional insulator, we can 
map it to an adiabatic pump sequence in a one-dimensional insulator by demoting 
one of the wavenumbers to time. The connection between the Chern number and 
the number of edge states in the higher dimensional Hamiltonian is a direct conse¬ 
quence of the connection between Chern number and charge pumping in the lower 
dimensional system. 

Chern insulators (two-dimensional band insulators with nonvanishing Chern 
number) were first used to explain the Quantum Hall Effect. There an external mag¬ 
netic field, included in lattice models via a Peierls substitution, is responsible for the 
nonzero value of the Chern number. Peierls substitution, however, breaks the lattice 
translation invariance. This neccessitates extra care, including the use of magnetic 
Brillouin zones whose size depends on the magnetic field. 

The models we construct in this Chapter describe the so-called Quantum Anoma¬ 
lous Hall Effect. Here we have the same connection between edge states and bulk 
Chern number as in the Quantum Hall Effect, however, there is no external mag¬ 
netic field, and thus no complications with magnetic Brillouin zones. The Quantum 
Anomalous Hall Effect has recently been observed in thin films of chromium-doped 
(Bi,Sb)2Te3 0. 

To illustrate the concepts of Chern insulators, we will use a toy model introduced 
by Qi, Wu and Zhang ll24ll . which we call the QWZ model. This model is also impor¬ 
tant because it forms the basic building block of the Bernevig-Hughes-Zhang model 
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for the quantum spin Hall effect (Chapt.j^, and thus it is also sometimes called “half 
BHZ”. 


6.1 Dimensional extension: from an adiabatic pump to a Chern 
insulator 


We want to construct a two-dimensional lattice Hamiltonian H with a nonvanishing 
bulk Chern number. We will do this by first constructing the bulk momentum-space 
Hamiltonian H{kx,ky), from which the real-space Hamiltonian can be obtained by 
Fourier transformation. For the construction we simply take an adiabatic pump se¬ 
quence on a one-dimensional insulator, H{k,t), and reinterpret the cyclic time vari¬ 
able f as a new momentum variable ky. This way of gaining an extra dimension by 
promotion of a cyclic parameter in a continuous ensemble to a momentum is known 
as dimensional extension. This, and the reverse process of dimensional reduction, 
are key tools to construct the general classification of topological insulators 1^ . 


From the Rice-Mele model to the Qi-Wu-Zhang model 


To see how the construction of a Chern insulator works, we take the example of 
the smooth pump sequence on the Rice-Mele model from the previous Chapter, 
Eqs. ( |4.7| l. In addition to the promotion of time to an extra wavenumber, Qt ^ ky, 
we also do an extra unitary rotation in the internal Hilbert space, to arrive at the 
Qi-Wu-Zhang model, 

H{k) = sink;cd’j:-l-sinkj,dv-l-[M + coskj:-|-cosky]a'j. (6.1) 


The mapping is summarized in Table 6.1 


Adiabatic pump in the RM model 

QWZ model (Chern Insulator) 

average intracell hopping v 
wavenumber k 
time t [in units of T/ (2k)] 

tTt, C7y, 

staggered onsite potential u 
wavenumber C 
wavenumber ky 


Table 6.1 Mapping of an adiabatic pump sequence of the Rice-Mele model, H{k,t) to the QWZ 
model for the Anomalous Hall Effect, H{kx,ky). 


The corresponding d(k) vector reads. 


d(kY, ky) 


( sinkx 

sinky 

u + coskx + cos 



( 6 . 2 ) 
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6.1.1 Bulk dispersion relation 

We can find the dispersion relation of the QWZ model using the algebraic prop¬ 
erties of the Pauli matrices, whereby = £(14)12, with I 2 the unit matrix. Thus, 
the spectrum of the QWZ model has two bands, the two eigenstates of H{k), with 
energies 


E±{k,,ky)=±\A{h,ky)\ (6.3) 

= ±\J sin^(^;c) +sin^(^y) -I- {u + cos{k^)+cos{ky))^. (6.4) 


The spectrum of the QWZ model is depicted in Fig |6.1| 

There is an energy gap in the spectrum of the QWZ model, which closes at fine- 
tuned values of M = -1-2,0, —2. This is simple to show, since the gap closing requires 
d(k) = 0 at some k. From Eq. ( |6.2[ ), £4(k) = £/y(k) = 0 restricts us to four inequiv¬ 
alent points in the Brillouin zone; 

• if M = —2: at kx = ky = 0, the F point; 

• if M = 0; at A:;f = 0, ky = n and k^ = 7t,ky = 0, two inequivalent X points 

• if M = +2: at k^ = 7i,ky — n, the M point; note that kx = ±7r, ky = ±n are all 
equivalent 


In the vicinity of a gap closing point, called Dirac point, the dispersion relation has 
the shape of a Dirac cone, as seen in Fig. 6 .1 For all other values of m 7 ^ —2,0,2, the 
spectrum is gapped, and thus it makes sense to investigate the topological properties 
of the system. 


6.1.2 Chern number of the QWZ model 


Although we calculated the Chern number of the corresponding pump sequence in 
the previous chapter, we show the graphical way to calculate the Chern number of 
the QWZ model. We simply count how many times the torus of the image of the 
Brillouin zone in the space of d contains the origin. To get some feeling about the 
not so trivial geometry of the torus, it is instructive to follow a gradual sweep of 


the Brillouin zone in Fig. 6.2 The parameter u shifts the whole torus along the 


direction, thus as we tune it we also control whether the origin is contained inside 
it or not. Generally three situations can occur as it is also depicted in Fig. |6.3| The 
torus either does not contain the origin as in (a) and (d) and the Chern number is 
2 = 0 for \ u\ > 2 , or we can take a straight line to infinity from the origin that pierces 
the torus first from the blue side (outside) of the surface as in (b) with Q = —1 for 
—2 < M < 0 , or piercing the torus from the red side (inside) as in (c) with 2 = 1 for 
0 <u<2. 

To summarize, the Chern number Q of the QWZ model is 
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(a)u = -2 (b)M=0 



Fig. 6.1 The bulk dispersion relation of the QWZ model, for various values of u, as indicated 
in the plots. In (a)-(c), the gapless cases are shown, where the bulk gap closes at so-called Dirac 
points. In (d), a generic value u = —1.8, the system is insulating. 


u < —2 

0 = 0; 

(6.5a) 

-2 < M < 0 

0 = -i; 

(6.5b) 

0 < M < 2 

0 = +i; 

(6.5c) 

2 < u 

0 = 0. 

(6.5d) 


6.1.3 The real-space Hamiltonian 


We obtain the full Hamiltonian of the Qi-Wu-Zhang model by inverse Fourier trans¬ 
form of the bulk momentum-space Hamiltonian, Eq. (|6.1|l, as 
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Fig. 6.2 The surface d(k) for the QWZ model as k sweeps through the whole Brillouin zone. 
To illustrate how this surface is a torus the sweeping is done gradually with m = 0. In (a) the 

image of the ky = —n line is depicted. In (b) the image for the region ky = —n -O.Str, in (c) 

ky = —7t -0.2571:, 'v[i(A)ky = —K -0, in (e) = — TT • • •0.257r, in {f) ky = -K - • -O.Stt, in (g) 

ky = -Tt- ■■ Q.lSn and finally in (h) the image of the whole Brillouin zone is depicted and the torus 
is closed. 



Fig. 6.3 The torus d(k) of the QWZ model for different values of u keeping A = 1. For clarity 
only the image of half of the Brillouin zone is shown. In (a) and (d) u = =f 2.2 and the torus does 
not contain the origin hence (2 = 0. In (b) m = — 1, taking an infinite line from the origin along the 
positive z axis we hit the blue side of the torus once hence Q = —\.\n(c)u = 1, taking the infinite 
line in the negative z direction we hit the red side of the torus thus 2=1. 
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N,-l Ny / 

Y i\mx + l,my){mjc,my 
mjf=lm_y=l V 



N, Ny-l / 

+ Y i\mjc,my + l){mjc,my 

m.v=lmy=l \ 



N, Ny 

+ 1* Y \mx,my) {mjc,my\iSiO^- 

mx=l my=\ 


( 6 . 6 ) 


As sketched in Fig. 6.4 1 , the model describes a particle with two internal states 
hopping on a lattice where the nearest neighbour hopping is accompanied by an 
operation on the internal degree of freedom, and this operation is different for the 
hoppings along the x and y directions. In addition, there is a staggered onsite po¬ 
tential of strength u. Unlike in the case of the SSH model (Chapt.[T]), the real-space 
form of the QWZ Hamiltonian is not intuitive. 


6.2 Edge states 




Fig. 6.4 Sketch of the QWZ model: a particle with two internal states hopping on a square lattice, 
(a): The onsite potential and the hopping amplitudes are operators acting on the internal states, (b): 
A strip, with periodic boundary conditions along y, open boundaries along x (hopping amplitudes 
set to zero). Light blue / dark red highlights left/right edge region, (c): Upon Fourier transformation 
along y, the strip falls apart to an ensemble of one-dimensional Hamiltonians, indexed by ky. 


We constructed a Chern insulator using an adiabatic charge pump. As we saw in 
Chapt.]^ charge pumps also induce energy eigenstates at the edge regions that cross 
from negative energy to positive energy bands, or vice versa. What do these energy 
eigenstates correspond to for the Chern insulator? 
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Dispersion relation of a strip shows the edge states 


To see edge states, consider a strip of a two-dimensional insulator depicted in 


Fig. 6.4 Along y, we take periodic boundary conditions (close the strip to a cylin¬ 
der), and go to the limit Ny —> Along x, the strip is terminated by setting the 
hopping amplitudes to 0 (open boundary condition), and it consists of N sites. Trans¬ 
lation invariance holds along y, so we can partially Fourier transform - only along 
y. After the Fourier transformation, the original Hamiltonian has falls apart to a set 
of one-dimensional lattice Hamiltonians indexed by a continuous parameter ky, the 
wavenumber along y. For the QWZ model, Eq. ( |6.6| ), the ky -dependent Hamiltonian 
reads 


H{ky)= ^ +h.c^ + 

;%=l V ^ / 

N, 

^ |mjc) (mj:| (g) (coskyff^-l-sinkydyM® dj). (6.7) 

m.v=l 

Note that this is the same dimensional reduction argument as before, but for a sys¬ 
tem with edges. Energy eigenstates \'F{ky)') of the strip fall into the categories of 
bulk states and edge states, much as in the one-dimensional case. All states are de¬ 
localized along y, but bulk states are also delocalized along x, while edge states are 
exponentially confined to the left (x = 0) or the right (x = N) edge. If we find energy 
eigenstates with energy deep in the bulk gap, they have to be edge states, and can be 
assigned to the left or the right edge. 

An example for the dispersion relation of a strip is shown in Eig. |6.5[ edge 
states on the left/right edge highlighted using dark red/light blue. We used the QWZ 
model, strip width N = 10, sublattice potential parameter m = —1.5, and the same 
practical definition of edge states as in the Rice-Mele model, 

l'F(kv)) is on the right edge ^ ^ \{'I'{ky)\ mx,a)\ >0.6; (6.8) 

mx=N-l ae{A,B} 

2 , 

l'f'(kv)) is on the left edge ^ ^ |('f'(ky) | >0.6. (6.9) 

nir=l ae{A,B} 


Edge states conduct unidirectionally 

Notice the edge state branches of the dispersion relation of the QWZ strip, Eig. |6.5[ 
which connect the lower and upper band across the bulk gap. They are the edge 
states of the pumped Rice-Mele model, Eig. |4.6| but we now look at them with a 
new eye. Eor the edge states in the QWZ model, dE/dky corresponds to the group 
velocity along the edge. Thus, the dispersion relation tells us that particles in the 
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QWZ model at low energy are confined either to the to left edge and propagate 
towards the right, or to the right edge and propagate towards the left. 

The presence of one-way conducting edge state branches implies that the QWZ 
model is no longer, strictly speaking, an insulator. Because of the bulk energy gap, 
it cannot conduct (at low energies) between the left and right edges. However, it will 
conduct along the edges, but only unidirectionally. 



wavenumber cell index m 

Fig. 6.5 Dispersion relation of a strip of QWZ model, of width N = W, and sublattice potential 
parameter u = —1.5. Because the strip is translation invariant along the edge, the wavenumber 
ky is a good quantum number, and the energy eigenvalues can be plotted (a) as a function of ky, 
forming branches of a dispersion relation. Light blue/dark red highlights energies of edge states, 
whose wavefunction has over 60% weight on unit cells with nix <1! mx>N — \ . These are parts 
of the Mh and W-h 1th branch, which split off from the bulk around —n/A<ky<n/A, and have 
an avoided crossing with an exponentially small gap near ky = 0. We show the marginal position 
probability distribution of the Mh energy eigenstate, PN{inx) = Lae{A,B}L«i. for 

three values of ky. Depending on ky, this state can be an edge state on the right edge (b), on the left 
edge (c), or a bulk state (d). 


6.2.1 Edge states and edge perturbation 

We can use dimensional reduction and translate the discussion about the robustness 
of edge states from Chapt.|^to the edge states of the QWZ model. This treats the 
case where the Hamiltonian is modified in a way that only acts in the edge regions, 
and is translation invariant along the edges. As an example, we introduce an ex¬ 
tra, state-independent next-nearest neighbor hopping, and onsite potentials at the 
left and right edge of the sample. As Fig |6.6| shows, this can modify the existing 
edge state branches, as well as create new edge state branches by deforming bulk 
branches. Including the new local terms the Hamiltonian of Eq. is augmented 
to read 
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H{h)= L 

mx= 


L (6.10) 

mxe{\,N} 

Where /r(0/(A') is the onsite potential on the left/right edge and describes a 

second nearest neighbor hopping on the left/right edge. If not stated otherwise these 
terms are considered to be zero. 

In the top row of Fig |6.6| the spectrum of strips without edge perturbations are 
depicted for Chern number Q = —1 (a.) and 2 = 0 (b) respectively. As we expected, 
a nonzero Chern number results in edge states, one on each edge. In (c) and (d) 
switching on perturbations, we see new edge states moving in to the gap. The onsite 
potential acts as an overall shift in energy on the states around nix = N, the second 
nearest neighbor hopping adds a considerable warping to the states localized around 
nix = 1 The deformations can change the number of edge states at a specific energy, 
but only by adding a pair of edge states with opposite propagation directions. This 
leaves the topological invariant unchanged. 


\mx + l) (otjcI C)- - - \-h.c. j + 

Nx 

^ \mx) {ifixl (S> {cosky.a^ + sin^ydyM® 6’j)+ 

my=l 






Fig. 6.6 Dispersion relation of a strip of the QWZ model, with edge states on the left/right edge 
in highlighted in dark red / light blue. Top row: clean system. Bottom row: with extra next nearest 
neighbor hopping along y, only the left edge, with amplitude =2 and an onsite potential /i A) 
at the right edge. The value of the potential in (c) is = 0.5 and in (d) it is = 1 .5. In (e) the 
schematics of the considered perturbations is shown. The additional potential terms do not affect 
the bulk states but distort the edge modes and bring in new edge modes at energies that are in the 
bulk gap or above/below all bulk energies. 
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Fig. 6.7 Layering sheets of 2-dimensional insulators on top of each other is a way to construct a 
two-dimensional insulator with higher Chem numbers. For uncoupled layers, the Chem numbers 
can simply be summed to give the total Chern number of the 3-layer structure, Gi + 62 + Qs- 
Switching on coupling (hopping terms) between the sheets cannot change the Chem number as 
long as the bulk gap is not closed. 


6.2.2 Constructing models with higher Chern numbers 

A systematic way to construct models with higher Chern numbers is to layer sheets 
of Chern insulators onto each other, as illustrated in Fig. |6.7| The single-particle 
Hilbert space of the composite system of D layers is a direct sum of the Hilbert 
spaces of the layers. 


Jifb =^L2® ■ ■J^LD- ( 6 . 11 ) 

The Hamiltonian, including a state-independent interlayer coupling with ampli¬ 
tude C, is 


D D-l 

Hd=Y, I^) ®HLd+ Y, (I'^+l) (“^1 + 1“^) + (6.12) 

d=l d=l 

with the unit operator on the Hilbert space of a single layer. The operators 

Htd we consider below are of the form of Eq. ( |6.10| l with different values of m, and 
can have an overall real prefactor. In layer d the strength of the local edge potential 
is denoted as As an example, the matrix of the Hamiltonian of a system 

with three coupled layers reads 


7/3 = 


Hl\ Cl 0 
Cl Hl 2 Cl 
0 Cl Hl3 


(6.13) 


Numerical results for two and three coupled layers, with different Chern numbers 


in the layers, are shown in Figs. 6.8 The coupling of copropagating edge modes lifts 
the degeneracies, but cannot open gaps in the spectrum, except in the case of strong 
coupling. This is a simple consequence of the fact that an energy eigenstate has to 
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be a single valued function of momentum. To open a gap, counterpropagating edge 
states have to be coupled. We achieve this by coupling layers of the QWZ model 


with opposite sign of the Chern number Q. For the case of two layers (Fig. 6.8 


second row), this opens a gap in the spectrum. If there are three layers, there is a 
majority direction for the edge states, and so one edge state survives the coupling. 



Fig. 6.8 Dispersion relations of strips of multilayered QWZ model, N = IQ unit cells wide. In 
all cases the bulk Hamiltonian of the first layer Hu is the QWZ Hamiltonian with u = —1.2. 
To elucidate the interplay of the edge states we consider a finite potential acting on the edge 
= 0.2/—0.2. In the left column ((a) and (d)), we have two layers, with bulk Hti = —/Til 
and edge onsite potential In the middle ((b), (e)) and right ((c), (f)) columns, 

we have three layers. The third layer is characterized in the bulk by = Hi^\ and on the edge by 
In (b) and (e) Hi 2 = 2Hii, while in (d) and (e) Hu, = —'^■Hu. For all four 
cases = 0. Coupling in (d), (e) and (f) is uniform with magnitude C = 0.4. 


6.3 Robustness of edge states 

Up to now, we have considered clean edges, i.e., two-dimensional Chern Insulators 
that were terminated by an edge (at lUx = 1 and = N), but translationally invariant 
along the edge, along y. This translational invariance, and the resulting fact that 
the wavenumber ky is a good quantum number, was used for the definition of the 
topologic al invariant N+ — N -, which was the net number of edge bands propagating 
along the edge, equal to the bulk Chern number, Q. With disorder in the edge region 
that breaks translational invariance along y, we no longer have a good quantum 
number ky, and edge state bands are not straightforward to define. However, as we 
show in this section, the edge states must still be there in the presence of disorder, 
since disorder at the edges cannot close the bulk gap. 
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Fig. 6.9 A disordered sample 
of Chem insulator. The dotted 
lines indicate rectangular 
parts of the sample, where 
disorder can be turned off 
adiabatically to reveal edge 
states (indicated in black). 
Since particles cannot turn 
back (unidirectional, or chiral 
channels), and cannot go 
into the bulk (in the gap), 
they have to travel all the 
way on the perimeter of the 
disordered sample, coming 
back to the rectangular, clean 
part. 



Smoothly removing disorder 

Consider a finite, sample of a Chem insulator, with a clean bulk part but disordered 
edge region, as depicted in Fig. |6.9| The bulk gap of the sample decreases due to dis¬ 
order, but we suppose that it is not closed completely (just renormalized). Consider 
now a small part of the sample, containing some of the edge, indicated by the dot¬ 
ted rectangle on the right of Fig. |6.9| Although this is much smaller than the whole 
sample, it is big enough so that part of it can be considered as translation invariant 
“bulk”. Now in this small part of the sample, we adiabatically (and smoothly) de¬ 
form the Hamiltonian in such a way that we set the disorder gradually to 0. This 
includes straightening the part of the open boundary of the sample that falls into the 
dotted rectangle, to a straight line. The deformation is adiabatic in the sense that the 
bulk gap is not closed in the process. Since this small part is a clean Chern insulator, 
with a bulk Chern number of Q, it can be deformed in such a way that the only edge 
states it contains are \Q\ states propagating counterclockwise (if Q>Q, say). 


Unitarity; particles crossing the clean part have to go somewhere 

Consider a particle in an edge state in the small clean part of the sample, with energy 
deep inside the bulk gap (E « 0). What can its future be, as its time evolution follows 
the Schrodinger equation appropriate for this closed system? Since the edge state is 
a chiral mode, the particle has to propagate along the edge until it leaves the clean 
region. Because of unitarity, the particle cannot “stop” at the edge: that would mean 
that at the “stopping point”, the divergence of the particle current in the energy 
eigenstate is nonzero. In other words, the particle current that flows in the edge state 
has to flow somewhere. (Put differently, if the mode describing an edge state particle 
“stopped at the interface”, two particles, initially orthogonal, following each other 
in the mode, would after some time arrive to the same final state. This would break 
unitarity.) After leaving the clean part of the sample, the particle cannot propagate 
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into the bulk, since its energy is deep in the bulk gap. The disorder in the clean part 
was removed adiabatically, and thus there are no edge states at the interface of the 
clean part and the disordered part of the sample, along the dashed line. The particle 
cannot turn back, as there are no edge states running “down” along the edge in the 
clean part. The only thing the particle can do is propagate along the edge, doing a 
full loop around the sample until it comes back to the clean part from below again. 

The argument of the previous paragraph shows that even though the sample is 
disordered, there has to be a low energy mode that conducts perfectly (reflection¬ 
less) along the edge. Since at 0 energy there are Q orthogonal states a particle can 
be in at the edge of the clean part of the sample, unitarity of the dynamics of the 
particles requires that all along the edge of the disordered sample there are Q or¬ 
thogonal modes that conduct counterclockwise. There can be additional low energy 
states, representing trapped particles, or an equal number of extra edge states for 
particles propagating counterclockwise and clockwise. However, the total number 
of counterclockwise propagating edge modes at any part of the edge always has to 
be larger by Q than the number of clockwise propagating edge modes. Because the 
Hamiltonian is short range, our conclusions regarding the number of edge states at 
any point far from the deformed region have to hold independent of the deformation. 

To be precise, in the argument above we have shown the existence of Q edge 
states all along the edge of the sample, except for the small part that was adia¬ 
batically cleaned from disorder. One way to finish the argument is by considering 
another part of the sample. If we now remove the disorder adiabatically only in this 
part, we obtain the existence of the edge modes in parts of the sample including the 
original dotted rectangle, which was not covered by the argument of the previous 
paragraph. 


Problems 

6.1. Phase diagram of the anisotropic QWZ model 

The lattice Hamiltonian of the QWZ model is provided in Eq. ( |6.6| l. Consider the 
anisotropic modification of the Hamiltonian when the first term of Eq. ( |6.6| ), de¬ 
scribing the hopping along the x axis, is multiplied by a real number A. Plot the 
phase diagram of this model, that is, evaluate the Chern number as a function of two 
parameters u and A. 



Chapter 7 

Continuum model of localized states at a domain 
wall 


So far, we have discussed edge states in lattice models, in which the states live on 
discrete lattice sites, and the Hamiltonian governing the physics is a matrix. In this 
chapter, we argue that in certain cases, it is also possible to describe these states via 
a continuum model, in which the states live in continuous space, and the Hamilto¬ 
nian is a differential operator. One benefit of such a continuum description is that it 
allows one to use the vast available toolkit of differential equations for solid-state 
problems in general, including the description of topologically protected states in 
particular. Another interesting aspect of these continuum models is their strong sim¬ 
ilarity with the Dirac equation describing relativistic fermions. A limitation of the 
continuum models is that their validity is restricted to narrow windows in momen¬ 
tum and energy; typically they are applied in the vicinities of band edges. Here, 
we obtain the continuum differential equations for three basic lattice models: the 
one-dimensional monatomic chain, the one-dimensional SSH model, and the two- 
dimensional QWZ model. In the cases of the SSH and QWZ models, the resulting 
equations will be used to analytically characterize the localized states appearing at 
boundaries between regions with different topological invariants. Even though in the 
entire chapter we build our discussion on the three specific lattice models, the appli¬ 
cability of the technique introduced here, called envelope-function approximation, 
is more general and widely used to describe electronic states in various crystalline 
solids. 


7.1 one-dimensional monatomic chain in an electric potential 

We use this minimal lattice model to illustrate the basic concepts of the envelope- 
function approximation (EFA) 0 , the technique that allows us to find the continuum 
versions of our lattice models. Of course, the one-dimensional monatomic chain 
does not host topologically protected states. 
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7 Continuum model of localized states at a domain wall 


7.1.1 The model 

We take a long lattice with ^ 1 unit cells (or sites) without an internal degree 
of freedom, with periodic boundary conditions, and a negative hopping amplitude 
f < 0. We consider the situation when the electrons are subject to an inhomogeneous 
electric potential V{x). This setup is pictured in Fig. HU The lattice Hamiltonian 
describing this inhomogeneous system reads 

H,^H + V, (7.1) 


where 


N 

H = t\m) {m+l\+h.c., (7.2) 

m=l 

N 

V = '£v„,\m){m\, (7.3) 

m=l 

with Vm = V{x = n). Our aim is to construct a continuum model that accurately 
describes the low-energy eigenstates of this lattice Hamiltonian. 


' periodic boundary condition ^ 

'^-o—o—o—o—o—o—o—o—o—o-' 

m = 1 m — N 

t < 0 
hopping 

V{x) 



m = 1 m = N 

Fig. 7.1 one-dimensional monatomic chain in an inhomogeneous potential V (x). 


Before discussing the inhomogeneous case incorporating V (x), focus first on the 
homogeneous system. The bulk momentum-space Hamiltonian is a scalar (1x1 
matrix) in this model, since there is no integral degree of freedom associated to the 
unit cell; it reads 


H{k) = e{k)\k) {k\ 

where e{k) is the electronic dispersion relation: 


(7.4) 















7.1 one-dimensional monatomic chain in an electric potential 


95 


The low-energy part of the dispersion relation is located around zero momentum, 
and is approximated by a parabola; 

t.2 

e(^)«-2|f| + |f|^2^eo + —, (7.6) 

where we introduced the minimum energy of the band ep = —2|f |, and the effective 
mass characterizing the low-energy part of the dispersion relation. (Using 

proper physical units, the effective mass would have the form m* = f?' j (2|f ja^), with 
a being the lattice constant.) For simplicity, we suppress ep in what follows; that is, 
we measure energies with respect to Eq. 


7.1.2 Envelope-function approximation 

Our goal is to find the low-energy eigenstates of the inhomogeneous lattice Hamil¬ 
tonian Hi. The central proposition of the EFA, applied to our specific example of the 
one-dimensional monatomic chain, says that it is possible to complete this goal by 
solving the simple continuum Schrodinger equation 

HuvK<P{x)=E<p{x), (7.7) 

where the envelope-function Hamiltonian Hepa has a very similar form to the free- 
electron Hamiltonian with the electric potential V: 

Hefa = -f V(x). (7.8) 

Here p = —idx is the usual real-space representation of the momentum operator, and 
the function <p(x) is usually called envelope function. Note the very simple relation 
between the low-energy dispersion in Eq. and the kinetic term in the EEA 
Hamiltonian ( |7.8| i: the latter can be obtained from the former by substituting the 
momentum operator p in place of the momentum k. 

Before formulating the EEA proposition more precisely, we introduce the concept 
of a spatially slowly varying envelope function. We say that <p(x) is spatially slowly 
varying, if its Eourier transform 


fN g-iqx 

^{q) = dx-^(p{x) (7.9) 

is localized to the |^| <C 7Z region, i.e., to the vicinity of the center of the BZ. 

With this definition at hand, we can formulate the EEA proposition. Consider the 
inhomogeneous one-dimensional monatomic chain described by Hi. Assume that 
(p{x) is a spatially slowly varying eigenfunction of //efa with eigenvalue E. Then, 
the state | y/) defined on the lattice via 
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N 


\¥) = L <Pix = m)\m). 


(7.10) 


is approximately an eigenstate of the lattice Hamiltonian Hi with eigenvalue E. The 
proof follows below in Sect. [7X3] 

Note that if the envelope function (p{x) fulfils the normalization condition 



(7.11) 


then, for a long lattice N ^ I, the lattice state |y/) will also be normalized to a good 
accuracy; 


(V/|vf>« 1. 


(7.12) 


It is important to point out two possible interpretations of the lattice state | \j/) 
introduced in Eq. ( |7.10[ ). (1) The state 11//) can be interpreted as the zero-momentum 
band-edge eigenstate |A: = 0) = I'”) homogeneous system, modu¬ 

lated by the envelope function (p{x) restricted to the lattice-site positions x = m. (2) 
The state jy/) can also be interpreted as a wave packet, composed of those eigen¬ 
states \k) of the homogeneous lattice Hamiltonian H that have wave numbers k close 
to the band-edge wave number, the latter being zero in this case. To see this, we first 
Fourier-decompose (p{x): 


Jkx / Jkx 


(7.13) 


In the approximate equality, we used the fact that ^(x) is spatially slowly varying, 
i.e., its Fourier transform ^{k) is localized to the central part of the BZ, and intro¬ 
duced the notation Ylk for ^ wave-number sum that goes only for the central part of 
the BZ. By inserting Eq. ( |7.13| ) to Eq. ( |7.10| l, we find 



(7.14) 


That is, I y/) is indeed a packet of plane waves with small wave numbers. 


7.1.3 Envelope-function approximation: the proof 


To prove the EFA proposition, we calculate Hi \ xj/) and utilize the spatially-slowly- 
varying condition on (p{x). Start with the contribution of the bulk Hamiltonian H: 
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H\xif)= ^e(^)|^)(^| £(p(m)|m) . (7.15) 

k€BZ m=l 

Express the envelope function via its Fourier transform: 

1 r ^ / piqm\ 

H\xif)= ^e(fe)|fe)(fe| E E l«) • (7.16) 

itGBZ J [m=l V^GBZ J 

Now use the fact that (p{x) is spatially slowly varying; that implies that the q sum 
can be restricted to the central part of the BZ: 

1 r A' / / piqm\ 

H\w)= EeWI^)(^l E E^(?)^ H • (7.17) 

jtGBZ J [m=l \ q VN J 

Performing the sum for m, we find 


//|V/)= Ee(^)l^)(^l E^(?)l?) • (7-18) 

jtGBZ J L 9 

Performing the scalar product yields 

HW)=Y^e{qmq)\q). ( 7 . 19 ) 

q 

Using the fact that the q sum goes for the central part of the BZ, where the dispersion 
relation e{q) is well approximated by a parabola, we find 

HW) ■ (7-20) 

Utilizing the definition of the plane wave |^), we obtain 



( 7 . 21 ) 


which can be rewritten as 



( 7 . 22 ) 

( 7 . 23 ) 

( 7 . 24 ) 
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Continue with the contribution of the potential V. Using Eqs. ( |73| ) and ( |TT0l ), 
we hnd 


V\w) = 

‘ N 

Y \m) {m\ 

Y (P{m') \m') 


m=l 

N 

_m'=\ 


m=l 


N 

= Y, \m) 

m=l 


(7.25) 


Summing up the contribution ( |7.24[ ) of H and the contribution ( |7.25| l of V, we 
hnd 


(// + y) IV7) « f [//efa(P(^)],=„, \m) = f |m) (7.26) 

m=l m=l 

N 

= £ ^ (p(m)|m)=£|v/), (7.27) 

m=l 


which concludes the proof. 


7.2 The SSH model and the one-dimensional Dirac equation 


To illustrate how the EFA captures topologically protected bound states in one¬ 
dimensional, we use the SSH model described in detail in Chapt. [T] The model is 
visualized in Fig. o The bulk Hamiltonian is characterized by two parameters, the 
intra-cell and inter-cell hopping amplitudes v, w > 0, respectively. The bulk lattice 
Hamiltonian of the SSH model is given in Eq. ( |l.ll i, whereas the bulk momentum- 
space Hamiltonian is given in Eq. ( |1.14| ). We have seen that the (v, w) parameter 
space is separated to two adiabatically connected partitions by the v — w line. In 
Sect. 1.5.1| we have also seen that localized zero-energy states appear at a domain 
wall between two half-inhnite homogeneous regions, if the two regions have dif¬ 
ferent bulk topological invariants; that is, if the sign of v — w is different at the two 
sides of the domain wall. 

This is the phenomenon that we address in this section: we show that an analytical 
description of such localized states can be given using the EFA. First, we discuss the 
electronic dispersion relation of the metallic (v = w) and nearly metallic (| v — w| <C 
|v-|-w|) homogeneous SSH model. Second, we obtain the Dirac-type differential 
equation providing a continuum description for the inhomogeneous SSH model (see 
Fig. |7.2| l for the energy range in the vicinity of the bulk band gap. Finally, we solve 
that differential equation to hnd the localized zero-energy states at a domain wall. 
Remarkably, the analytical treatment remains useful even if the spatial structure of 
the domain wall is rather irregular. 
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-H-^^ ^^^^^^^-►a: 

TO. = 1 m = N 


Fig. 7.2 Inhomogeneous intracell hopping and domain walls in the SSH model. The dashed el¬ 
lipse denotes the unit cell. The dashed line connecting the edges of the chain denotes the periodic 
boundary condition. 


7.2.1 The metallic case 

First, consider the metallic homogeneous SSH model, where v = w. The dispersion 
relation is shown as the blue solid line in Fig. The hlled and empty bands 
touch at the end of the BZ, at k — ko = 7Z. Figure s shows that in the vicin¬ 
ity of that touching point, commonly referred to as a Dirac point, the dispersion 
relations are linear functions of the relative wave vector q = k — ko. The slope of 
these linear functions, corresponding to the group velocity of the electrons, can 
be determined, e.g. by Taylor-expanding the bulk momentum-space Hamiltonian 
H{k) = (v-l- wcosk)d’j: -f wsinkdy, see Eq. ( |1.10| i, to hrst order in q: 

H{ko + q) Pi-wqdy, (v = w). (7.28) 

which indeed has a linear dispersion relation, 

E±{q) = ±wq. (7.29) 

The eigenstates of the linearized Hamiltonian ( |7.28| l are 

*-*(,) = (7.30) 

Note that the dispersion relation of the Dirac equation of fermions with zero mass 
is 

E±{k)=±hkc, (7.31) 

where Ti is the reduced Planck’s constant and c is the speed of light. Comparing 
Eqs. ( |7.29| l and ( |7.31[ ), we conclude that the dispersion of the metallic SSH model 
is analogous to that of massless Dirac fermions, and the hopping amplitude of the 
metallic SSH model plays the role of he. Because of the similarity of the dispersions 
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\1.29) and ( |7.31| l, the linearized Hamiltonian ( |7.28| ) is often called a massless Dirac 
Hamiltonian. 

At this point, the linearization of the bulk momentum-space Hamiltonian of the 
SSH model does not seem to be a particularly fruitful simplihcation; to obtain the 
dispersion relation and the eigenstates, a 2 x 2 matrix has to be diagonalized, no mat¬ 
ter if the linearization has been done or not. However, linearizing the Hamiltonian 
is the hrst step towards the EFA, as discussed below. 


7.2.2 The nearly metallic case 


Now consider a homogeneous, insulating SSH model that is nearly fnetallic; that 
is, the scale of the energy gap |v — w| opened at ko is signihcantly smaller than the 
scale of the band width v + w. An example is are shown in | 1.2) 7, where the dispersion 
relation is plotted for the parameter values v = 1 and w = 0.6. 

We wish to describe the states close to the band gap located around zero energy. 
Hence, again, we can use the approximate bulk momentum-space Hamiltonian ob¬ 
tained via linearization in the relative momentum q: 

H{ko + q) ^ Mdx — wqdy, (7.32) 


where we dehned M = v — w. The dispersion relation reads 

E±{q) = ±s/M^ + w^q^. (7.33) 


The (unnormalized) eigenstates of the linearized Hamiltonian (|7.32|i have the form 



M -f iwq 

E±iq) 


(7.34) 


Note that the dispersion relation of the Dirac equation for fermions with hnite 
mass II ^0 reads 

E±{k) = (7.35) 

Therefore, the parameter M = v — w of the SSH model plays the role of the mass- 
related term jj.c^ of the relativistic dispersion relation ( |7.35| l, and the linearized 
Hamiltonian (|7.32|l is often called a massive Dirac Hamiltonian. 


7.2.3 Continuum description of the nearly metallic case 

We are mostly interested in a continuum description of the zero-energy localized 
states formed at a domain wall between two topologically distinct regions. For sim¬ 
plicity and concreteness, consider the case when the domain wall is created so that 
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the intra-cell hopping amplitude v varies in space while the inter-cell one is constant, 


as shown in Fig. 7.2 


In what follows, we will focus on one of the two domain walls shown in Fig. 7.2 
The inhomogeneous Hamiltonian has the form 


N N 

Hi= Vmi\m,B) {m,A\+h.c.) + w {\m,B) {m+l,A\+h.c.), (7.36) 

m=\ m=\ 


where Vm = v(.x = m) and v{x) > 0 is a continuously varying function of position, 
which takes the constant value v_ (v+) far on the left (right) from the domain wall. 

We also assume that the local Hamiltonian is a nearly metallic SSH Hamiltonian 
everywhere in space. That is, |v(x) — w| ^ v(.r:) +w. This ensures that the local band 
gaps |v± — w\ on the two sides of the domain wall are much smaller than the local 
band widths v± + w. 

Based on our experience with the EFA in the inhomogeneous one-dimensional 
monatomic chain (see Sect. HD’ now we construct the EFA proposition corre¬ 
sponding to this inhomogeneous, nearly metallic SSH model. Recall that in the 
former case, we obtained the EFA Hamiltonian by (i) Taylor-expanding the bulk 
momentum-space Hamiltonian around the wave vector cotTesponding to the band 
extremum (that was ko = 0 in Sect. IlD’ (ii) replacing the relative wave vector q 
with the momentum operator p = —idx, and (iii) incorporating the inhomogeneity 
of the respective parameter, which was the on-site potential V (x) in that case. The 
same procedure, applied now for the SSH model with a hrst-order Taylor expansion, 
yields the following EFA Hamiltonian: 


Hefa = M{x)dx - wpdy. (7.37) 

The EFA proposition is then formulated as follows. Assume that <p(x) = {(Pa{x),(Pb{x)) 
is a spatially slowly varying eigenfunction of Hefa in Eq. 31) , with eigenvalue E. 
Then, the state | yr) dehned on the lattice via 

\V)^Y L (7.38) 

m=\ a=A,B 

is approximately an eigenstate of the lattice Hamiltonian Hi with energy E. 

Note that, in analogy with Eq. ( |7.14| i, the lattice state | xj/) can be reformulated 
as a wave packet formed by those eigenfunctions of the homogeneous (v(x) = w) 
system that are in the vicinity of the band-edge momentum kq: 

/ 

\V)-Y L Mq)\ko + q)<3)\a). 

q a=A,B 


(7.39) 
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7.2.4 Localized states at a domain wall 


Having the envelope-function Schrodinger equation 

[M{x)dx — wpdy] (p{x) =E(p{x) (7.40) 


at hand, we can study the domain wall between the two topologically distinct re¬ 
gions. First, we consider a step-type domain wall, defined via 


M{x) 


Mo if x> 0, 

—Mq if X < 0 ’ 


and Mq > 0, as shown in Fig. |7.3^. 


(7.41) 



Fig. 7.3 (a) Step-like and (b) irregular spatial dependence of the mass parameter M{x) of the 
one-dimensional Dirac equation. 


We wish to use the EFA Schrodinger equation ( |7.40| i to establish the zero-energy 
states localized to the domain wall, which were revealed earlier in the lattice SSH 
model. That is, we look for evanescent solutions of Eq. ( |7.40[ l on both sides of the 
domain wall, and try to match them at the interface x = Q. Eor the .r > 0 region, our 
evanescent-wave Ansatz reads 


(Px>0{x) 



(7.42) 


with K >0. Substituting this to Eq. ( |7.40| l yields a quadratic characteristic equation 
for the energy E, having two solutions 


E± 


= ± 


■sj Mq — f 


(7.43) 


The corresponding unnormalized spinors read 
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a± 




' Mq — WK ' 
E± 


An analogous Ansatz for the x < 0 region is 


^x<q{x) = 



(7.44) 


(7.45) 


with K >0, yielding the same energies as in Eq. (|7.43|l, and the spinors 


c± 

d± 



(7.46) 


Now consider an energy eigenstate with a given energy E. For clarity, set Mq > 
E >0. The (unnormalized) envelope function of the energy eigenstate has the form 

(p{x) = (p^yo{x)©{x)+C(p^<o{x)©{-x), (7.47) 

where (px<o and (Px>o should be evaluated by replacing E± >—>■ E and K i—^ 
and C is a yet unknown parameter to be determined from the boundary conditions 
at the domain wall. The envelope function ( |7.47| l is an eigenstate of the EFA Hamil¬ 
tonian with energy E if the boundary condition that the wave function is continuous 
at X = 0, that is, 

(Px<o{0) = C(px>oiO), (7.48) 

is fulhlled. Note that in our case, the Dirac equation is a hrst-order differential equa¬ 
tion and therefore there is no boundary condition imposed on the derivative of the 
wave function. From the second component of Eq. ( |7.48| l, we have C — 1. From the 
hrst component, we have Mq — wK — —Mq -f wK, implying Mq = wK and thereby 
E — 0. The same result is obtained if the range —Mq < £ < 0 of negative energies is 
considered. Hence we conclude that the zero-energy state at the domain wall does 
appear in the continuum model of the inhomogeneous SSH chain, as expected. 

Let us also determine the coefficients a and c describing this localized state: 

Mq-JmI-e^ 

a = lim-^-= 0, (7.49) 

E^O E 

and similarly, c = 0. These imply that the localized state is completely sublattice- 
polarized, i.e., it lives on the B sublattice, and therefore it is its own chiral partner. 
Considering the a similar mass prohle with negative Mq, we would have found that 
the localized state lives on the A sublattice. These properties are in line with our 
expectations drawn from the lattice SSH model. 

A further characteristic property of the localized state is its localization length; 
from our continuum model, we have an analytical result for that: 
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1 W 

K~ Mo' 


(7.50) 


(In physical units, that is ^ = ^a.) Recall that we are constrained to the nearly 
metallic regime w ^ Mq = v+ — w; together with Eq. ( |7.50| l, this implies that the 
localization length is much larger than one (that is, the lattice constant). This is 
reassuring: it means that the envelope function <p(x) is spatially slowly varying, 
hence is within the range of validity of the EFA. 

The result ( 7.50| l can be compared to the corresponding result for the SSH lattice 
model. Eq. ( |1.50| l provides the localization length ^ of an edge state in a disordered 
SSH model, which corresponds to the localization length ^ obtained above. Taking 
the disorder-free special case of Eq. (|1.50|l, we have 




1 


w 


log^ log 


w+(v—w) 


w — v 


(7.51) 


As we are making a comparison to the nearly metallic v « w case considered in 
this section, we could approximate ^ in Eq. ( |7.51| l using a leading-order Taylor 
expansion in the small quantity (w — v)/w. The approximate result ( |7.51| ) is in line 
with Eq. ( |7.50| ) obtained from the continuum model. 

A further interesting fact is that the existence of the localized state is not con¬ 
strained to the case of a sharp, step-like domain wall described above. The simple 
spinor structure found above also generalizes for less regular domain walls. To see 
this, consider an almost arbitrary one-dimensional spatial dependence M{x) of the 
mass, illustrated in Fig. | 7 . 3 | 7 , with the only condition that M changes sign between 
the half-planes x < 0 and x > 0, i.e., M(x — —o°) < 0 and M(x —>^ oo) > 0. We claim 
that there exists a zero-energy solution of the corresponding one-dimensional Dirac 
equation that is localized to the domain wall and has the envelope function 


(p{x) = 



(7.52) 


To prove this claim, insert this wave function ^(x) to the one-dimensional Dirac 
equation and substitute E = 0 therein. This procedure results in the single differ¬ 
ential equation <3jc/(x) = — implying that Eq. ( |7.52| ) is indeed a zero-energy 
eigenstate of the envelope-function Hamiltonian if the function / has the form 

/(x) = const X » Iodx'M(x') ^ ^7 53 ^ 

Furthermore, the asymptotic conditions of the mass M(x) guarantee that this enve¬ 
lope function decays as x —± 00 , and therefore is localized at the domain wall. 
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7.3 The QWZ model and the two-dimensional Dirac equation 

We have introduced the QWZ model as an example for a two-dimensional Chern 
insulator in Chapt. The lattice Hamiltonian of the model is given in Eq. ( | 6 . 6 | ), 
whereas the bulk momentum-space Hamiltonian is given in Eq. The dispersion 
relation is calculated in Eq. ( |6.4| i, and examples of it are shown in Eig. 

Recall that the model has a single parameter u, and the Chern number of the 
model is determined by the value of u via Eq. ( |6.5[ ). Similarly to the case of the SSH 
model in one dimension, one can consider a domain wall between locally homoge¬ 
neous regions of the QWZ model that have different Chern numbers. Just as the edge 
of a strip, such a domain wall can support topologically protected states that prop¬ 
agate along the domain wall but are localized at the domain wall in the transverse 
direction. The number and propagation direction of those states is determined by 
the magnitude and sign of the difference of the Chern numbers in the two domains, 
respectively. In this section, we use the EEA to provide a continuum description of 
such states. 


6.1 


7.3.1 The metallic case 


Eirst, consider the metallic cases of the QWZ model; that is, when the band structure 
has no energy gap. In particular, we will focus on the u = —2 case. The correspond¬ 


ing band structure is shown in Eig. 6.1 'a). The two bands touch at k = (0,0), and 
form a Dirac cone at that Dirac point. 

To describe excitations in the vicinity of the Dirac point of such a metal, it is 
sufficient to use a linearized approximation of the QWZ Hamiltonian H(k) that is 
obtained via a Taylor expansion of i/(k) up to first order in the k-space location 
q = k — ko measured from the Dirac point kq. In the case u = —2, the Dirac point is 
ko = (0,0), and the linearized Hamiltonian reads 


H{ko + q) « qxdx -f qjdy 


(7.54) 


The dispersion relation is E± (q) = ±^; again, this is analogous to that of the mass¬ 
less Dirac equation Eq. (|7.31|l. 


7.3.2 The nearly metallic case 

Now consider a QWZ insulator that is nearly metallic: u « —2. The dispersion rela¬ 
tion for M = — 1.8 is shown in Eig. |6.1D d). In the vicinity of the metallic state, as seen 
in the figure, a small gap opens in the band structure at the Dirac point ko = (0,0). 

The states and the band structure around ko can again be described by a linearized 
approximation of the QWZ Hamiltonian i/(ko + q) in q: 






106 


7 Continuum model of localized states at a domain wall 


i/(ko + q) + + (7.55) 

where we defined the parameter M = u+ 2. The dispersion relation reads 

E±{q) = ±^/AfiW- (7.56) 

A comparison with the relativistic dispersion relation ( |7.35[ ) reveals that the param¬ 
eter M of the QWZ model plays the role of /rc^; hence M can be called the mass 
parameter. 


7.3.3 Continuum description of the nearly metallic case 

We have discussed that the QWZ lattice with an inhomogeneous u parameter might 
support topologically protected states at boundaries separating locally homogeneous 
regions with different Chern numbers. Similarly to the one-dimensional SSH model 
treated above, these localized states can be described analytically, using the envelope 
function approximation (EFA), also in the two-dimensional QWZ model. In the rest 
of this Chapter, we focus on the nearly metallic case where the inhomogeneous 
u{x,y) is in the vicinity of —2 (i.e., \M{x,y)\ = \u{x^y) -|-2| ^ 1), in which case 
the the low-energy excitations are expected to localize in Fourier space around the 


band extremum point ko = (0,0) (see Figs. 6.1 1 and d). Here we obtain the FFA 
Schrodinger-type equation, which resembles the two-dimensional Dirac equation, 
and in the next subsection we provide its localized solutions for simple domain-wall 
arrangements. 

The considered lattice is inhomogeneous due to the spatial dependence of the 
parameter u{x,y). In the tight-binding lattice model, we denote the value of u in unit 
cell m = {mx,my) as = u{x = mx,y = my), and correspondingly, we introduce 
the local mass parameter via M{x,y) = u{x,y) -f 2 and Mm = M{x = mx,y = my). 

The EFA Hamiltonian can be constructed the same way as in sections |7.1| 
and 7.2.3 The bulk momentum-space Hamiltonian //(ko +q) is Taylor-expanded 


around the band-edge wave vector ko = (0,0), the wave-number components qx and 
qy are replaced by the differential operators px and py, respectively, and the inho¬ 
mogeneous mass parameter M(.r,y) is incorporated. This yields the following result 
in our present case: 

Hepk = M{x,y)d^ + pxdx + Py&y (7.57) 

Then, the familiar EFA proposition is as follows. Assume that the two-component 
envelope function (p{x,y) is a spatially slowly varying solution of the EFA Schrodinger 
equation 

HEFA(pix,y) =E(p{x,y). (7.58) 

Then, the lattice state | y/) associated to the envelope function (p{x,y) is defined as 


IV^) = L ^a(m)|m, 


a) 


(7.59) 
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It is claimed that the lattice state | y/) is approximately an eigenstate of the inhomo¬ 
geneous lattice Hamiltonian with the eigenvalue E. 


7.3.4 Chiral states at a domain wall 


We can now use the EFA Hamiltonian ( |7.57| l to describe the chiral states at a domain 
wall between two topologically distinct regions of the QWZ model. The Dirac-type 
EFA Schrodinger equation reads: 

[M{x,y)a^ + + pdy\<p{x,y) = E<p{x,y). (7.60) 

Consider the homogeneous case hrst: M(x,y) =Mq, where Mq might be positive 
or negative. What is the dispersion relation for propagating waves? What are the 
energy eigenstates? The answers follow from the plane-wave Ansatz 


^(x,y) 


( ^\ pkxX Jqyj 

\b) 


(7.61) 


with qx,qy G K and a,b G C. With this trial wave function, Eq. ( |7.60| ) yields two 
solutions: 


E± = ±JMl + ql + qj, 


(7.62) 


and 


_ qx ^qy 

b± E± — Mq 


(7.63) 


Describe now the states at a domain wall between two locally homogeneous re¬ 
gions where the sign of the mass parameter is different. Remember that the sign 
of the mass parameter in the EFA Hamiltonian is related to the Chern number of 
the corresponding homogeneous half-BHZ lattice: in our case, a positive (negative) 
mass implies a Chern number —1 (0). 

To be specihc, we will consider the case when the two domains are dehned as 
the y < 0 and the y > 0 half-planes, i.e., the mass prohle in Eq. (|7.60|l are 


M{x,y) 


Mq if y > 0, 
—Mo if y < 0 


(7.64) 


Let Mo be positive; the corresponding mass prohle is the same as shown in Fig. |7.3^ , 
with X replaced by y. 

Now we look for solutions of Eq. ( |7.60| ) that reside in the energy range —Mq < 
E < Mo, i.e., in the bulk gap of the two domains, and which propagate along, but 
decay perpendicular to, the domain wall at y = 0. Our wave-function Ansatz for the 
upper half plane y > 0 is 




(Pu{x,y) 


( 7 . 65 ) 
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with qx S K, C and a,b G C. For the lower half plane, (pi{x,y) is defined 

as <Pi,(x,y) but with u gg I interchanged and qy '^ G The wave function <Pi,{x,y) 
does solve the two-dimensional Dirac equation defined by Eqs. ( |7.60| l and ( |7.64| l in 
the upper half plane y>0 provided 


4 “' = iK = i^Ml+ql-E^ 


(7.66) 


and 


“t" It* 

bu E-Mo 


(7.67) 


Similar conditions apply for the ansatz (pi (v,y) for the lower half plane, with the sub¬ 
stitutions ut-G- 1, K i-G —K and Mo i—> —Mo- The complete (unnormalized) envelope 
function has the form 


(p{x,y) = %{x,y)&{y) + c(pi{x,y)©{-y), (7.68) 

where c is a yet unknown complex parameter to be determined from the boundary 
conditions at the domain wall. 

The wave function ( |7.68| l is an eigenstate of the EFA Hamiltonian with energy 
E if the boundary condition that the wave function is continuous on the line y = 0, 
that is, 

(p„{x,0) = (pi{x,0), (7.69) 

is fulfilled for every x. 

The boundary condition ( |7.69| ) determines the value of the parameter c as well as 
the dispersion relation £ of the edge states. First, ( |7.69| l implies 

qx-K = c(qx -f K*) ^ c = ^ , (7.70) 

qx + K 

E -f Mo = c{E-Mo) ^ -qxMo = kE. (7.71) 

Note that K depends on E according to Eq. ( |7.66| ). It is straigthforward to find the 
dispersion relation of the edge states by solving —qxMo = k{E)E for E with the 
condition —Mq < E < Mo'. 

E = -qx. (7.72) 

This simple dispersion relation is shown in Eig. |7.4^ . Together with Eq. ( |7.66| l, this 
dispersion implies that the localization length of edge states is governed by Mq only, 
i.e., is independent of qx. The squared wave function of an edge state is shown in 
Eig.|74]3. 

A remarkable consequence of this simple dispersion relation is that the spinor 
components of the envelope function also have a simple form: 



(7.73) 



















7.3 The QWZ model and the two-dimensional Dirac equation 


109 


(“) E 





Fig. 7.4 Chiral state obtained from the two-dimensional Dirac equation, (a) Dispersion relation 
and (b) squared wave function of a chiral state confined to, and propagating along, a mass domain 
wall. 


Edge states at similar mass domain walls at u{y) « 0 and u{y) « 2 can be derived 
analogously. Note that at u{y) « 0, the low-energy states can reside in two different 
Dirac valleys, around kq = (0, n) or kq = (?t,0), and there is one edge state in each 
valley. The number of edge states obtained in the continuum model, as well as their 
directions of propagation, are in correspondence with those obtained in the lattice 
model; as we have seen for the latter case, the number and direction are given by 
the magnitude and the sign of Chern-number difference across the domain wall, 
respectively. 

An interesting fact is that the existence of the edge state is not constrained to 
case of a sharp, step-like domain wall described above. Moreover, the simple disper¬ 
sion relation and spinor structure found above generalize for more irregular domain 
walls. This generalization is proven in a similar fashion as in the case of the SSH 


model, see Sect. 7.2.4 To see this, consider an almost arbitrary one-dimensional 
spatial dependence of the mass, similar to the one in Fig. 7.3 r. M{x,y) — M{y) with 
the only condition that M changes sign between the half-planes y < 0 and y > 0, i.e., 
M(y ——oo) < 0 and M(y -^ °°) > 0. We claim that there exists a solution of the 
corresponding two-dimensional Dirac equation that propagates along the domain 
wall, has the dispersion relation E = —qx, is confined in the direction perpendicular 
to the domain wall, and has the wave function 


^(x,y) = 


1 

-1 




m- 


(7.74) 


To prove this proposition, insert this wave function (p{x,y) to the two-dimensional 
Dirac equation and substitute E with —qx therein. This procedure results in two 
equivalent equations that are fulfilled if dyf{y) = —M{y)f{y), implying that Eq. 











110 


7 Continuum model of localized states at a domain wall 


( |7.74| l is indeed a normalizable solution with E = —q^ provided that the function / 
has the form 

/(y) =e-/o^>’W). (7.75) 

To summarize: In the preceding chapters, we introduced the topological charac¬ 
terization of lattice models and the corresponding edge states and states bound to 
domain walls between regions of different topological character. In this chapter, we 
demonstrated that a low-energy continuum description (the EFA Schrodinger equa¬ 
tion) can be derived from a lattice model, and can be used to describe those elec¬ 
tronic states. Besides being a convenient analytical tool to describe inhomogeneous 
lattices, the envelope-function approximation also demonstrates that the emergence 
of topologically protected states is not restricted to lattice models. 


Problems 


7.1. SSH model with spatially dependent intracell hopping. 

In Sect. 7.2.3 we provide the EFA proposition for the SSH model with spatially de¬ 
pendent intracell hopping. Prove this proposition, following the procedure detailed 
in Sect. l7.1.3l for the one-dimensional monatomic chain. 


7.2. SSH model with spatially dependent intercell hopping. 

Derive the EFA Hamiltonian for an inhomogeneous SSH model, where w varies in 
space and v is constant. Assume a nearly metallic scenario, w{x) « v. 

7.3. QWZ model. 

Prove the EFA proposition for the QWZ model. The proposition is outlined in 
Sect. |7.3.3] The proof is analogous to that used for the SSH model. 

7.4. QWZ model at m « 2. 

The bulk momentum-space Hamiltonian H{k) of the QWZ model is given in Eq. 
(j^. Starting from this H{k), derive the EFA Hamiltonian describing low-energy 
excitations in the case of an inhomogeneous u parameter for which u « 2. 









Chapter 8 

Time-reversal symmetric two-dimensional 
topological insulators - the 
Bernevig-Hughes-Zhang model 


In the previous chapters, we have seen how two-dimensional insulators can host 
one-way propagating (a.k.a. chiral) edge states, which ensures reflectionless trans¬ 
port along the edge. The existence of chiral edge states precludes time-reversal sym¬ 
metry. Indeed, time-reversed edge states would describe particles propagating back¬ 
wards along the edge. In Chern Insulators (two-dimensional insulators with non¬ 
vanishing Chern number), the absence of these counterpropagating states from the 
spectrum is what ensures the reflectionless propagation of particles along the edges. 

What about time-reversal symmetric or time-reversal invariant two dimensional 
insulators? According to the above, they cannot be Chern insulators. Interestingly 
though, the same time-reversal symmetry that ensures that for every edge state mode 
there is a counterpropagating time reversed partner, can also ensure that no scatter¬ 
ing between these two modes occurs. This means that it is possible for time-reversal 
invariant two-dimensional insulators to host edge states with reflectionless propaga¬ 
tion, in both directions, at both edges. The details of why and how this happens are 
discussed in this and the following chapters. 

We will find that all two-dimensional time-reversal invariant insulators fall into 
two classes: the trivial class, with an even number of pairs of edge states at a single 
edge, and the topological class, with an odd number of pairs of edge states at a single 
edge. We then subsequently show that disorder that breaks translational invariance 
along the edge can destroy edge state conduction in the trivial class, but not in the 
topological class. 

The bulk-boundary correspondence for Chern Insulators stated that the the net 
number of edge states on the edge, N+ —N-, is the same as the Chern number of 
the bulk, Q. We showed this by mapping the 2-dimensional system to a periodically, 
adiabatically pumped one-dimensional chain. After the mapping, the unit of charge 
pumped through the chain during a period could be identified with the net number 
of chiral edge states. 

Unfortunately, identifying and calculating the bulk topological invariant of a 
time-reversal invariant two-dimensional insulator is much more cumbersome than 
for a Chern insulator. We therefore come back to this problem in the next Chapter. 


Ill 
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8.1 Time-Reversal Symmetry 

Before we discuss time-reversal symmetric topological insulators, we first need to 
understand what we mean by time reversal symmetry, and how it leads to Kramers’ 
degeneracy. 


8.1.1 Time Reversal in continuous variable quantum mechanics 
(without spin) 

Take a single particle with no internal degree of freedom, described by a wave- 
function 'f'(r). Its dynamics is presribed by a time independent Hamiltonian H = 
{p — eA(f))^ +V{r), where the functions A and V are the vector and scalar poten¬ 
tials, respectively, and e is the charge of the particle. The corresponding Schrodinger 
equation for the wavefuncion ^{r,t) reads 

idt'F{r,t) = {(-/5r-eA(r))^-fy(r)}'f'(r,t). (8.1) 

Any solution W of the above equation can be complex conjugated, and gives a solu¬ 
tion of the complex conjugate of the Schrodinger equation, 

= {{-idr + eA{r)f+V{r)}^{r,ty. (8.2) 


The operator of complex conjugation in real space basis 

We use K to denote the operator that complex conjugates everything to its right in 
real space basis. 


A/(r)=/(r)*A; K^ = l, (8.3) 

for any complex valued function /(r) of position. The Schrodinger equation above 
can be rewritten using K as 


Kidt'F = KidtKKW = -id,W* = KHKK'P = H*'F*. (8.4) 

Complex conjugation in real space basis conforms to intuitive expectations of 
time reversal: it is local in space, takes x —>■ x, and flips the momenta, idj^ —> —idjc- 


Time reversal 

The above relation shows that for any closed quantum mechanical system, there is a 
simple way to implement time reversal. This requires to change both the wavefunc- 
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lion W to f'* and the Hamiltonian H to H*. The change of the Hamiltonian involves 
A —> —A, i.e., flipping the sign of the vector potential. 


Time reversal symmetry 

In the special case where the Hamiltonian in real space basis is real, H* = H, we 
can implement time reversal by only acting on the wavefunction. In that case, we 
say that the system has time reversal symmetry. For the scalar Schrodinger equation 
above, this happens if there is no vector potential, A = 0. To see this more explicitly, 
consider time evolution for a time f, then apply the antiunitary operator K, then 
continue time evolution for time f, then apply K once more: 

U = (8.5) 

If H* = H, then U = 1, which means that K acts like time reversal. 


8.1.2 Lattice models with an internal degree of freedom 

In these notes we deal with models for solids which are lattice Hamiltonians: the 
position (the external degree of freedom) is discrete, and there can be an internal 
degree of freedom (sublattice, orbital, spin, or other). 


Definition of the operator K 

For the operator of complex conjugation we need to fix not only the external position 
basis, (c^externai = {Im)}, but also an internal basis, (^internal = {|ct)}- The property 
defining K then reads 

Vz G C!, V |m) , |m ^ G (^xtemah ^ 1^) ? ^ internal • 

Kz |m, a) (m', a'\ = z* |m, a) (m', a' \ K, (8.6) 

where z* is the complex conjugate of z. 

We will use the shorthand jT'*) and A* to represent A jf') and KAK, respectively. 
The defining equations are 
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(8.7) 

m a 


'F*)=K\'F)=Y,J^'F,:,„\m)^\a)-, 

(8.8) 

m a 


^ ) (^^1 ^ | OJ ) (OJ | , 

m'ma'a 

(8.9) 

A* =KAK= ^ Y.Kt'.a',m,aW)ifn\(^\a'){a\. 

m'm a'a 

(8.10) 


Time-reversal affects external and internal degrees of freedom 

We look for a representation of time reversal symmetry in terms of a general an¬ 
tiunitary operator. Apart from complex conjugation, which acts on both external and 
internal Hilbert space, we allow for an additional unitary operation on the internal 
degrees of freedom f, that is independent of position, 

^ = tK. ( 8 . 11 ) 

We say that a Hamiltonian H is time reversal invariant (or time reversal symmet¬ 
ric) with respect to time reversal represented by if 

^H^ = H. ( 8 . 12 ) 

In the same sense as for the chiral symmetry (cf. Sect. |1.4| l, when we talk about a 
Hamiltonian, what we really mean is a set of Hamiltonians ), with ^ represent¬ 
ing parameters that are subject to disorder. Thus, Eq. ( |8.12| l should hold for any of 
the H{^), with ^ independent of ^. 


8.1.3 Two types of time-reversal 

We can require that a time reversal operator when squared, should give at most 
a phase; 


XKxK = XX* = e''^Iintemal- (8.13) 

If that was not the case, if the unitary operator xx* was nontrivial, then it would 
represent a unitary symmetry of a time-reversal invariant Hamiltonian, since 

xrH{xrf = mxKHKx^^Kx^^ = mHKx^^ = H. (8.14) 

As explained in Sect. d, when we want to investigate topological phases, the 
usual first step is to get rid of unitary symmetries one by one (except for the lattice 
translation symmetry of the bulk Hamiltonian), by restricting our attention to a sin- 
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gle superselection sector of each symmetry. Thus, the only time reversal symmetries 
that are left are those that fulhl Eq. ( |8.13[ ). 

The phase factor turns out to have only two possible values; -f 1 or —1. 

Multiplying Eq. ( |8.13| l from the left by we get f* = where the 

superscript T denotes transposition. Iterating this last relation once more, we obtain 
f* = which means = ±1, wherefore 

,#^ = ±1. (8.15) 

A Hamiltonian with no unitary symmetries can have only one type of time- 
reversal symmetry: either = -f 1, or = —1, but not both. Assume a Hamil¬ 
tonian had two different time-reversal symmetries, ^ and . Along the lines of 
Eq. ( |8.14| i, the product of the two, the unitary operator would then represent 
a unitary symmetry. The only exception is if !^\ = when they are not re¬ 

ally different symmetries. However, in this case, since ^ is antiunitary, these two 
symmetries square to the same number, ^^ 

An example for a time-reversal operator with = -fl is given by the com¬ 
plex conjugation K. An example for a time-reversal operator with = — 1 is time 
reversal for a spin-1/2 particle. Since time reversal should also flip the spin, it is 
achieved hy ^7 = —iOyK, with K dehned on the basis of the eigenstates of The 
fact that this works can be checked by = —Oj for j = x,y,z. 


The operator f is symmetric or antisymmetric 

Specifying the square of the time-reversal operation constrains the operator f to be 
symmetric or antisymmetric. Consider 

= mtK = = ±1-, r = ±f+= (8.16) 

where the subscript T denotes transpose in the same basis where the complex con¬ 
jugate is defined. Therefore, 

=+l f = symmetric; (8.17) 

=—I f = —antisymmetric. (8.18) 


8.1.4 Time reversal of type = — 1 gives Kramers’ degeneracy 

A dehning property of an antiunitary operator ^ is that for any pair of states jT') 
and |^>), we have 

= (f = (0*1 •?*) = (0| 0)*. 

(8.19) 
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Consider now this relation with |^>) = ^ j'f'): 

= (±'F\^'F) = (8.20) 

where the ± stands for the square of the time reversal operator which is ±1. If 
= -f 1, the above line gives no information, but if = — 1, it leads immediately 
to I f') = 0, which means that for every energy eigenstate, its time-reversed 
partner, which is also an energy eigenstate with the same energy, is orthogonal. This 
is known as Kramers degeneracy. 


8.1.5 Time-Reversal Symmetry of a Bulk Hamiltonian 


We now calculate the effect of time-reversal symmetry ^ = zK on th e bulk 
momentum-space Hamiltonian H{k). This latter is obtained, as in Sect. 1.2 by first 


periodic boundary conditions, and dehning a plane wave basis in the corresponding 
external Hilbert space as 


|k) #|k) = |-k).#. (8.21) 

k 

Next, Hbuik is the part of H in the bulk, with periodic boundary conditions, whose 
components in the plane wave basis dehne the bulk momentum-space Hamiltonian, 

i?(k) = (k|i7b,ik|k); i/bulk = Llk)(k|®-^(k). (8.22) 

k 

The effect of time-reversal symmetry follows, 

>#f/buik^^‘ =l^|-k)(-k|(g)f//(k)*f^ = ^|k)(k|(g)fi/(-k)*f'*'. (8.23) 

k k 

We read off the action of ^ on the bulk momentum-space Hamiltonian, and obtain 
the necessary requirement of time-reversal symmetry as 

fi7(-k)*f^ =.H(k). (8.24) 

Note that time-reversal symmetry of the bulk Hamiltonian is necessary, but not suf- 
hcient, for time-reversal symmetry of the system: perturbations at the edges can 
break time reversal. 

A direct consequence of time-reversal symmetry is that the dispersion relation of 
a time-reversal symmetric Hamiltonian has to be symmetric with respect to inver¬ 
sion in the Brillouin zone, k —> — k. Indeed, take an eigenstate of H(k), with 


H(k)|M(k))=£(k)|M(k)). 
Using time-reversal symmetry, Eq. (|8.24|), we obtain 


(8.25) 
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\u{k)) = £(k) |M(k)) . (8.26) 

Multiplying from the left by and complex conjugating, we have 

|M(k))* = |M(k))* . (8.27) 

This last line tells us that for every eigenstate |M(k)) of H{k.), there is a time-reversed 
partner eigenstate of i7(— k) at the same energy, |M(k))*. This implies inversion 
symmetry of the energies, ^(k) —E (—k). Note, however, that E{k.) — E (—k) is not 
enough to guarantee time-reversal symmetry. 

It is especially interesting to look at points in the Brillouin zone which map 
unto themselves under inversion: the Time-reversal invariant momenta (TRIM). In 
d dimensions there are such points, one of which is at the center of the Brillouin 
zone (so-called F point), and others at the edges. At such momenta, Eq. ( |8.24| i 
implies that 


xH (kTRiivi)* x^' =H (kTRiM) • (8.28) 

If = — 1, then because of Kramers degeneracy, at a TRIM, every eigenvalue of 
the bulk momentum-space Hamiltonian is (at least) twice degenerate. 


8.2 Doubling the Hilbert Space for Time-Reversal Symmetry 


There is a simple way to construct a system with Time-Reversal Symmetry, Htrt 
starting from a lattice Hamiltonian H. We take another two copies of the system, 
and change the Hamiltonian in one of them to H* = KHK. We then couple them, 
much as we did to layer Chern insulators on top of each other in Sect. |6.2.2| 

^TRI = |0) (0| (g)H+ |1) (1| + (|0) (1| (g)Iextemal<8)C + /!-C.) , (8.29) 

where the hopping between the copies is accompanied by a position-independent 
operation C on the internal degree of freedom. In a matrix form, in real-space basis 
(and somewhat simplified notation), this reads 


77tri 


H C 
H* 


(8.30) 


We will use Sx^y^z to denote the Pauli operators acting on the “copy degree of 
freedom”, defined as 


^xjyjz — ^xjyjz ® ^external ® ^internal- 


(8.31) 


Using these operators, the time-reversal invariant Hamiltonian reads 
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^TRI = 


1 + iv 


)i/+- 


1 -Sz 




Sy ~\~ IS-v Sy IS-v 

^^«)Iexternal«)C+-^ ^ TT ■ • ^ 


■ 0 ^external 0 ^ • (8.32) 


The choice of the coupling operator C is important, as it decides which type of 
time-reversal symmetry i/xRi will have. 


Time reversal with = — 1 requires antisymmetric coupling operator C 


If we want a time-reversal symmetry that squares to —1, we can go for 

# = isyK; = K{-i)sy, (8.33) 


with the factor of i is included for convenience, so that the matrix of isy is real. The 
requirement of time-reversal symmetry can be obtained using 


{isyK)HTRi{isyK) * 


0 1 
-1 0 

H 

-H* -C 



H* 

C*' 


0 -T 



H 


1 0 



0 -1' 


■ H 



1 0 


-C* 

H* 


(8.34) 


We have time-reversal symmetry represented hy ^ = iSyK, if 


iSyKHj^yiiSyK)-^ = Ht:^i ^ C =-C^, (8.35) 


where the subscript T denotes transposition in the same fixed internal basis that is 
used to define complex conjugation K. 


Symmetric coupling operator C gives time reversal with = -fl 

If the coupling operator is symmetric, C = C^, then the same derivation as above 
shows that we have time-reversal symmetry represented by = sK, 

SxKHjriKSx = Hjri ^ C = C^. (8.36) 

This time-reversal operator squares to -fl. 

If all we want is a lattice Hamiltonian with a time-reversal symmetry that squares 
to 4-1, we don’t even need to double the Hilbert space. We can just take 

^ = K- Htri = = KH-yriK. (8.37) 

Colloquially, this construction is referred to as taking the real part of the Hamilto¬ 
nian. 
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8.2.1 A concrete example: the Bernevig-Hughes-Zhang model 

To have an example, we use the construction above to build a toy model - called 
BHZ model - for a time-reversal invariant topological insulator starting from the 
QWZ model of Chapt. We follow the construction through using the bulk 
momentum-space Hamilonian, and obtain 

^BHz(k) = so'^[{u + coskx + cosky)d^ + smkydy)] +s^‘^smkxdx + Sx‘^C, 

(8.38) 

where C is a Hermitian coupling operator acting on the internal degree of freedom. 
For C = 0, the Hamiltonian ^bhz reduces to the 4-band toy model for HgTe, intro¬ 
duced by Bemevig, Hughes and Zhang Q. 


Two Time-Reversal Symmetries if there is no coupling 

If there is no coupling between the copies, C = 0, the BHZ model has two differ¬ 
ent time-reversal symmetries, = iSyK and = SxK, due to its block diagonal 
structure reflecting a unitary symmetry, — H^rz- In this situation, the type 

of predictions we can make will depend on which of these symmetries is robust 
against disorder. We will in the following require the = — 1 symmetry. If this 
symmetry is robust, then everything we do will apply to the C = 0 case of Hbhz- 
The extra time reversal symmetry in that case is just a reminder that most features 
could be calculated in a more simple way, by working in the superselection sectors 
of Sj separately. 


8.3 Edge States in 2-dimensional time-reversal invariant 
insulators 

We now consider the situation of edge states in a two-dimensional lattice Hamil¬ 
tonian with time-reversal symmetry, much in the same way as we did for Chern 
insulators in Sect. 16.21 


8.3.1 An example: the BHZ model with different types of coupling 


We start with the concrete example of the BHZ model. We set the sublattice potential 
parameter m = — 1.2, and plot the edge dispersion relation, defined in the same way 
as for the Chern insulators in Section lh^ 

As long as there is no coupling between the two copies, C = 0, the system Hbhz 


is a direct sum of two Chern insulators, with opposite Chern numbers. As Fig. 8.1 (a) 
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Fig. 8.1 Stripe dispersion relations of the BHZ model, with sublattice potential parameter u = 
— 1.2. Right/left edge states (more than 60% weight on the last/first two columns of unit cells) 
marked in dark red/light blue, (a): uncoupled layers, C = 0. (b): Symmetric coupling C = 03(Jx 
gaps the edge states out. (c): Antisymmetric coupling C = 0.3(7,. cannot open a gap in the edge 
spetrum. 


shows, on each edge, there is a pair of edge state branches: a branch on the layer with 
Hamiltonian H, and a counterpropagating branch on the layer with H*. Although 
these two edge state branches cross, this crossing will not turn into an anticrossing: 
the states cannot scatter into each other since they are on different layers. The two 
edge state branches are linked by time-reversal: they occupy the same position, but 
describe propagation in opposite directions. In fact, they are linked by both time- 
reversal symmetries this system has, by S^K and iSyK. 

A co upling between the layers can can gap the edge states out, as shown in 


Fig. 8.1 (b). We here used C = 0.3(7^, which respects the = -fl symmetry but 


breaks the = — 1 one. The crossings between counterpropagating edge states on 
the same edge have turned into anticrossings, as expected, since the coupling allows 
particles to hop between the counterpropagating edge states (on the same edge, but 
in different layers). 

We see something different if we couple the layers whil e respecting the = — 1 
time reversal symmetry, by, e.g., C — 0.3(7y.. As Fig. 8.1 (c) shows, the crossing at 


ky = 0 between the edge state branches now does not turn into an anticrossing. As 
long as the coupling is not strong enough to close the bulk gap, the edge states here 
appear to be protected. 


8.3.2 Edge states in — — I 


The states form one-dimensional edge state bands in the one-dimensional Brillouin 


zone kx = shown schematically in Fig. 8.2 In general, an edge will host 

edge states propagating in both directions. However, due to time-reversal symme¬ 
try, the dispersion relations must be left-right symmetric when plotted against the 
wavenumber kx along the edge direction. This means that the number N+ of right- 
moving edge states (these are plotted with solid lines in Fig. |8.2| i, and A_, the num¬ 
ber of left-moving edge states (dashed lines) have to be equal at any energy. 
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Fig. 8.2 Edge states on a single edge of a 2-dimensional time-reversal invariant topological in¬ 
sulator with = — 1, as the edge region undergoes a continuous deformation, parametrized by 
a, respecting the symmetry and the translational invariance along the edge. In (a)-(d), the edge 
state dispersions relations are shown, in the full edge Brillouin zone and in the 

energy window corresponding to the bulk gap. For clarity, right- (left-) propagating edge states 
are denoted by continuous (dashed) lines. Due to the deformation of the Hamiltonian, the edge 
state branches can move, bend, and couple, while the bulk remains unchanged. From (a) to (b), 
the crossing points between counterpropagating edge states become anticrossings, i.e., gaps open 
in these pairs of dispersion relation branches as a usual consequence of any parameter coupling 
them. Crossings at = 0 and k^ = n cannot be gapped, as this would lead to a violation of the 
Kramers theorem. From (b) to (d), these gaps become so large that at energy £ = 0, the number 
of edge states drops from 6 (3 Kramers pairs) to 2 (1 Kramers pair). The k^ values of the edge 
states at 0 energy are plotted in (e), where this change in the number of edge states shows up as an 
“annihilation” of right-propagating and left-propagating edge states. 


N+{E)=N^{E). 


(8.39) 


As with Chem insulators, we next consider the effect of adiabatic deformations 
of the clean Hamiltonian on edge states. We consider terms in the Hamiltonian that 
conserve translational invariance along the edge, and respect Time Reversal Symme¬ 
try. The whole discussion of Sect. |6.3| applies, and therefore adiabatic deformations 
cannot change the signed sum of of the left- and right-propagating edge states in the 
gap. However, time-reversal symmetry restricts this sum to zero anyway. 

Time reversal symmetry that squares to = — 1, however, provides a further 
restriction: adiabatic deformations can only change the number of edge states by 
integer multiples of four (pairs of pairs). To understand why, consider the adiabatic 


deformation corresponding to Fig. 8.2 (a)-(d). Degeneracies in the dispersion rela¬ 
tion can be lifted by coupling the edge states, as it happens in (b), and this can lead to 
certain edge states disappearing at certain energies, as in (c). This can be visualized 
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by plotting the kx values at £ = 0 of the branches of the edge state dispersion as func¬ 
tions of the deformation parameter (which is some combination of the parameters 
of the Hamiltonian) a, as in Fig. 8.2 (e). Due to the deformation, two counterpropa- 
gating edge states can “annihilate”, when the corresponding modes form an avoided 
crossing. If this happens at a generic momentum value k, as in (c), then, due to the 
time reversal invariance, it also has to happen at —k, and so the number of edge 
states decreases by 4, not by 2. The special momentum values of kx are the TRIM, 
which in this case are kx = 0, ±n. If the edge state momenta meet at a TRIM, as in 
(b) at kx = 0, their “annihilation” would change the number of edge states by 2 and 
not by 4. However, this cannot happen, as it would create a situation that violates the 
Kramers degeneracy: at the TRIM, energy eigenstates have to be doubly degenerate. 
The deformations in Fig. |8.2| can be also read from (d) to (a), and therefore apply to 
the introduction of new edge states as well. 


8.3.3 Z 2 invariant: parity of edge state pairs 


At any energy inside the bulk gap, the parity of the number of edge-states Kramers 
pairs for a given dispersion relation is well dehned. In Fig. |8.2( a), there are 3 
edge-state Kramers pairs for any energy in the bulk gap, i.e., the parity is odd. In 
Fig. |8.2( c), there are 3 of them for every energy except for energies in the mini-gap 
of the bands on the left and right for which the number of edge-state Kramers pairs 
is 1, and for the upper and lower boundaries of the mini-gap [the former depicted 
by the horizontal line in Fig. |8.2( c)], where the number of Kramers pairs is 2. The 
parity is odd at almost every energy, except the two isolated energy values at the 
mini-gap boundaries. 

The general proposition is that the parity of the number of edge-state Kramers 
pairs at a given edge for a given Hamiltonian at a given energy is independent of 
the choice of energy, as long as this energy is inside the bulk gap. Since in a time- 
reversal invariant system, all edge states have counter-propagating partners, we can 
express this number as 


D = 


N{E) 

2 


mod 2 = 


N+{E)+N-{E) 

2 


mod 2, 


(8.40) 


where N{E) = N+{E) +N-{E) is the total number of edge states at an edge. A 
caveat is that there are a few isolated energy values where this quantity is not well 
dehned, e.g., the boundaries of mini-gaps in the above example, but these energies 
form a set of zero measure. 

Since D is a topological invariant, we can classify two-dimensional time-reversal 
invariant lattice models according to it, i.e., the parity of the number of edge-state 
Kramers pairs supported by a single edge of the terminated lattice. Because it can 
take on two values, this ‘label’ D is called the Z 2 invariant, and is represented by a 
bit taking on the value 1 (0) if the parity is odd (even). 
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As a final step, we should next consider disorder that breaks translational in- 
varance along the edge, in the same way as we did for Chem insulators. Due to the 
presence of edge states propagating in both directions along the edge, the treatment 
of disorder is a bit trickier than it was for Chern states. 


8.4 Absence of backscattering 

A remarkable property of Chern insulators is that they support chiral edge states, i.e., 
edge states that have no counter-propagating counterparts. A simple fact implied by 
the chiral nature of these edge states is that impurities are unable to backscatter a 
particle that occupies them. As we argue below, absence of backscattering is also 
characteristic of disordered two-dimensional time-reversal invariant topological in¬ 
sulators, although the robustness is guaranteed only against time-reversal symmetric 
scatterers. 

In this Chapter, we introduce the scattering matrix, a simple concept that allows 
for a formal analysis of scattering at impurites, and discuss the properties of edge 
state scattering in two-dimensional time-reversal invariant topological insulators. 
The scattering matrix will also serve as a basic tool in the subsequent chapter, where 
we give a simple theoretical description of electronic transport of phase-coherent 
electrons, and discuss observable consequences of the existence and robustness of 
edge states. 


8.4.1 The scattering matrix 

Consider a phase-coherent two-dimensional conductor with a finite width in the y 
direction, and discrete translational invariance along the x axis. Think of the system 
as having periodic boundary conditions in the x direction. As earlier, we describe 
the system in terms of a simple lattice model, where the unit cells form a square 
lattice of size Nx x Ny, and there might be an internal degree of freedom associated 
to the unit cells. 

As the system has discrete translational invariance along x, we can also think of it 
as a one-dimensional lattice, whose unit cell incorporates both the internal degree of 
freedom of the two-dimensional lattice and the real-space structure along the y axis. 
Using that picture, it is clear that the electronic energy eigenstates propagating along 
the X axis at energy E have a product structure, as required by the one-dimensional 
Bloch’s theorem: 

|l,±) = |ku)®|0,^±), (8.41) 

On the right hand side, the first ket corresponds to a usual momentum eigenstate 
“ VW ^mj=i l^r) propagating along x, whereas the second ket incorporates 
the shape of the transverse standing mode as well as the internal degree of freedom. 
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The states appearing in Eq. ( |8.41| l are normalized to unity. The integer I = 
labels the propagating modes, also referred to as scattering channels. The -|- and — 
signs correspond to right-moving and left-moving states, respectively; the direction 
of movement is assigned according to the sign of the group velocity; 

dE(k) 

= -IT 


, (8.42) 

k= 


where E{k) is the dispersion relation of the one-dimensional band hosting the state 

ii,±). 

Now we re-normalize the states above, such that different states carry the same 
particle current through an arbitrary vertical cross section of the system. We will use 
these current-normalized wave functions in the definition of the scattering matrix 
below, which guarantees that the latter is a unitary matrix. According to the one¬ 
dimensional relation we obtained between the current and the group velocity in Eq. 
( |5.13| l, the current-normalized states can be defined using the group velocity as 

\l,±), = ^^\l,±). (8.43) 


A 


a 


(in) 

L 


a 


(out) 

R 


y 


(out) 


(in) 

B. 


X 


>> 


Fig. 8.3 Disordered region (gray) obstructing electrons in a two-dimensional phase-coherent con¬ 
ductor. The scattering matrix 5 relates the amplitudes and of incoming waves to the 
amplitudes and of outgoing waves. 


Now consider the situation when the electrons are obstructed by a disordered 


region in the conductor, as shown in Eig. 8.3 A monoenergetic wave incident on the 


scattering region is characterized by a vector of coefficients 


,(“) — 


( (in) 


(in) (in) (in) (in) 

T,2 ’ • ■ ■ ’ ■ 




(8.44) 


The first (second) set of N coefficients correspond to propagating waves ( |8.43| ) in 
the left (right) lead L (R), that is, the clean regions on the left (right) side of the 
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disordered region. The reflected and transmitted parts of the wave are described by 
the vector 


out) _ ^(out) 

— 1 U/^ 1 , Lij^ 2 


L,1 


(out) (out) (out) (out) 


(8.45) 


The corresponding energy eigenstate reads 


\w) = 'La^j\l,+.L)c + atT\l-^^)c + afl\l,-,R)c + a^RT\l,+,R)c(^-^(>) 


1=1 


Here, the notation for the current-normalized states introduced in Eq. ( |8.43| l has 
been expanded by the lead index L/R. 

The scattering matrix S relates the two vectors introduced in Eqs. ( |8.44| l and 
( [8l31 l: 

(8.47) 

The size of the scattering matrix is 2N x 2N, and it has the following block structure: 


S = 



(8.48) 


where r and / are x reflection matrices describing reflection from left to left and 
from right to right, and t and f' are transmission matrices describing transmission 
from left to right and right to left. 

Particle conservation, together with the current normalization Eq. ( |8.43| l, implies 
the unitarity of the scattering matrix S. In turn, its unitary character implies that 
the Hermitian matrices tt\ t't'\ 1 — rr\ and 1 — all have the same set of real 
eigenvalues Ti ,7’2,...7]v, called transmission eigenvalues. 


8.4.2 A single Kramers pair of edge states 


Now we use the scattering matrix S to characterize defect-induced scattering of 
an electron occupying an edge state of a two-dimensional time-reversal invariant 
topological insulator. Consider a half-plane of such a homogeneous lattice which 
supports exactly one Kramers pair of edge states at a given energy E in the bulk gap, 
as shown in Eig. |8.4| Consider the scattering of the electron incident on the defect 


from the left side in Eig. 8.4 The scatterer is characterized by the Hamiltonian V. 
We will show that the impurity cannot backscatter the electron as long as V is time- 
reversal symmetric. 

We choose our propagating modes such that the incoming and outgoing states 
are related by time-reversal symmetry, i.e.. 




(8.49a) 

(8.49b) 
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G> 


<■ 



G> 


<■ 






Fig. 8.4 Scattering of an edge state on a time-reversal symmetric defect V. In a two-dimensional 
time-reversal invariant topological insulator with a single Kramers pair of edge states, the incoming 
electron is transmitted through such a defect region with unit probability. 


Also, recall that = — 1. In the presence of the perturbation V, the edge states of 
the disorder-free system are no longer energy eigenstates of the system. A general 
scattering state |v/) at energy E is characterized by the vector of incoming 
amplitudes. According to Eq. ( |8.46| l and the definition ( |8.47| ) of the scattering matrix 
S, the energy eigenstates outside the scattering region can be expressed as: 




(8.50) 


(Sa^^ 


in) 


L,l 






Using Eq. ( |8.49| ), we find 


(8.51) 


+ ^ |/, +,L), + ^ \l-,R), 


Due to time reversal symmetry, this state is also an energy eigenstate having 

the same energy as | y/). Using the unitary character of the scattering matrix, the state 
— ^ly/) can be rewritten as 


-^\v) 



\ / RJ 


(8.52) 


where denotes the transpose of S, is also an energy eigenstate having the same 
energy as \y/). Comparing Eqs. ( |8.50[ ) and ( |8.51| l, and knowing that the scattering 
matrix at a given energy is uniquely defined, we conclude that S = —S^, that is 


r t' 
t r' 



(8.53) 


implying 


r = r' = 0, 


(8.54) 























8.4 Absence of backscattering 


127 


and hence perfect transmission of each of the two incoming waves. 


Fig. 8.5 Backscattering of 
edge states at a constric¬ 
tion. The states forming the 
edge-state Kramers pairs are 
depicted as solid and dashed 
lines, (a) A time-reversal sym¬ 
metric defect localized to the 
edges, such as a small con¬ 
striction shown here, is unable 
to backscatter the incoming 
electron, (b) Backscattering 
is possible between different 
edges, if the width of the 
constriction is of the order of 
the decay length of the edge 
states, (c) A finite spatial gap 
between the left and right 
part of the wire implies zero 
transmission. 


(«). 


(&). 




(c) 




If the lattice has the geometry of a ribbon, and the TRS scatterer extends to both 
edges, then the absence of backscattering is not guaranteed. This is illustrated in 
Fig. 8.5 where we compare three examples. In (a), the defect is formed as a wide 
constriction on both edges, with a width much larger than the characteristic length 
of the penetration of the edge states to the bulk region of the ribbon. Backscatter¬ 
ing between states at the same edge is forbidden due to TRS, and backscattering 
between states at different edges is forbidden due to a large spatial separation of 
their corresponding wave functions. In (b), a similar but narrower constriction with 
a width comparable to the penetration length of the edge states does allow for scat¬ 
tering between states on the lower and upper edges. In this case, backscattering from 
a right-moving state on one edge to a left-moving state at the other edge is not for¬ 
bidden. In (c), the constriction divides the ribbon to two unconnected parts, resulting 
in zero transmission through the constriction. 

Naturally, backscattering is also allowed if the scatterer is not time-reversal sym¬ 
metric, or if the scattering process is inelastic. Backscattering is not forbidden for 
the ‘unprotected’ edge states of topologically trivial {D = 0) two-dimensional time- 
reversal invariant insulators. 
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8.4.3 An odd number of Kramers pairs of edge states 

The above statement ( |8.54| l implying unit transmission, albeit in a somewhat weak¬ 
ened form, can be generalized for arbitrary two-dimensional time-reversal invariant 
topological insulator lattice models, including those where the number of edge-state 
Kramers pairs N is odd but not one. The proposition is that in such a system, given a 
time-reversal symmetric scatterer V and an arbitrary energy E in the bulk gap, there 
exists at least one linear combination of the incoming states of energy E from each 
side of the defect that is perfectly transmitted through the defect. 

The proof follows that in the preceding section, with the difference that the quan¬ 
tities r and t describing reflection and transmission are N xN matrices, and that the 
antisymmetric nature of the 5-matrix S = —S^ implies the antisymmetry of the re¬ 
flection matrices r = —r^. According to Jacobi’s theorem, every odd-dimensional 
antisymmetric matrix has a vanishing determinant, which is implied by 

det(r) = det(r^) = det(—r) = (—l)^det(r) = —det(r), (8.55) 

where we used the antisymmetry of r in the second step and the oddness of N in 
the last step. As a consequence of Eq. ( |8.55| l we know det(r) = 0, hence det(r^r) = 
det(r^)det(r) = 0. Therefore, at least one eigenvalue of rV is zero, which implies 
that at least one transmission eigenvalue 7/ is unity. 


8.4.4 Robustness against disorder 


Fig. 8.6 A disordered two- 
dimensional time-reversal 
invariant topological insulator 
contacted with two electrodes. 
Disorder is ‘switched off’ and 
the edge is ‘straightened out’ 
within the dashed box, hence 
the edge modes there resem¬ 
ble those of the disorder-free 
lattice. 



The absence-of-backscattering result ( |8.541 i implies a remarkable statement re¬ 
garding the existence of (at least) one perfectly transmitting ‘edge state’ in a finite- 
size disordered sample of a two-dimensional time-reversal invariant topological in¬ 
sulator. 
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(See also the discussion about Fig. |6.9| in the context of Chern Insulators). Such 
a sample with an arbitrarily chosen geometry is shown in Fig. 8.6 Assume that 
the disorder is TRS and localized to the edge of the sample. We claim that any 
chosen segment of the edge of this disordered sample supports, at any energy that 
is deep inside the bulk gap, (at least) one counterpropagating Kramers pair of edge 
states that are delocalized along the edge and able to transmit electrons with unit 
probability. This is a rather surprising feature in light of the fact that in truly one¬ 
dimensional lattices, a small disorder is enough to induce Anderson localization of 
the energy eigenstates, and hence render the system an electronic insulator. 

To demonstrate the above statement, let us choose an edge segment of the disor¬ 
dered sample for consideration, e.g., the edge segment running outside the dashed 
box in Fig. |8.6| Now imagine that we ‘switch off’ disorder in the complementer part 
of the edge of the sample, and ‘straighten out’ the geometry of that complementer 
part, the latter being shown within the dashed box of Fig. 8.6 Furthermore, via 
an appropriate spatial adiabatic deformation of the Hamiltonian of the system in the 
vicinity of the complementer part of the edge (i.e., within the dashed box in Fig. 8 . 61 , 
we make sure that only a single edge-state Kramers pair is present within this com¬ 
plementer part. The existence of such an adiabatic deformation is guaranteed by 


the topologically nontrivial character of the sample, see the discussion of Fig. 8.2 


The disordered edge segment, outside the dashed box in Fig. 8.6 now functions as 
a scattering region for the electrons in the straightened part of the edge. From the 
result ( |8.54| ) we know that such a TRS scatterer is unable to induce backscattering 
between the edge modes of the straightened part of the edge, hence we must con¬ 
clude that the disordered segment must indeed support a perfectly transmitting edge 
state in each of the two propagation directions. 

In the next chapter, we show that the electrical conductance of such a disordered 
sample is hnite and ‘quantized’, if it is measured through a source and a drain con¬ 
tact that couple effectively to the edge states. 










Chapter 9 

The Z2 invariant of two-dimensional topological 
insulators 


In the previous Chapter, we have seen that two-dimensional insulators can host 
topologically protected edge states even if time-reversal symmetry is not broken, 
provided it squares to —1, i.e., = — 1. Such systems fall into two categories: no 

topologically protected edge states (trivial), or one pair of such edge states (topo¬ 
logical). This property defines a Z 2 invariant for these insulators. In the spirit of the 
bulk-boundary correspondence, we expect that the bulk momentum-space Hamilto¬ 
nian H(k) should have a corresponding topological invariant (generalized winding 
number). 

The bulk Z 2 invariant is notoriously difficult to calculate. The original definition 
of the invariant lfT^ [m uses a smooth gauge in the whole Brillouin zone, that is hard 
to construct ED, which makes the invariant difficult to calculate. An altogether dif¬ 
ferent approach, which is robust and calculatable, uses the scattering matrix instead 
of the Hamiltonian 113]. 

In this Chapter we review a definition of the bulk Z 2 invariant lJTl based on the 
dimensional reduction to charge pumps. This is equivalent to the originally defined 
bulk invariants IIJTI . but no smooth gauge is required to calculate it. It can be outlined 
as follows. 

1. Start with a bulk Hamiltonian H{kx,ky). Reinterpret ky as time: H{kx,ky) is a 
bulk one-dimensional Hamiltonian of an adiabatic pump. This is the dimensional 
reduction we used for Chern insulators in Chapt. 

2. Track the motion, with ky playing the role of time, of pumped particles in the 
bulk using Wannier states. We will call this the Wannier center flow. 

3. Time-reversal symmetry restricts the Wannier center flow. As a result, in some 
cases the particle pump cannot be turned off adiabatically - in those cases the 
insulator is topological. If it can be turned off, the insulator is trivial. 

In order to go through this argument, we will first gather the used mathematical 
tools, i.e., generalize the Berry phase and the Wannier states. We will then define 
the Wannier center flow, show that it can be calculated from the Wilson loop. Finally, 
we will use examples to illustrate the Z 2 invariant and argue that it gives the number 
of topologically protected edge states. 
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9.1 Tools: Nonabelian Berry phase, multiband Wannier states 

To proceed to calculate the topological invariants, we need to generalize the tools of 
geometric phases, introduced in Chapt.|^ and of the Wannier states of Chapt. [^to 
manifolds consisting of more bands. 


9.1.1 Preparation: Nonabelian Berry phase 

We defined the Berry phase in Chapt. as the relative phase around a loop L of 
states I'Pj), with j = 1,2,... ,N- Since the Berry phase is gauge independent, it is 
really a property of the loop over N one-dimensional projectors ) In most 
physical applications - in our case as well - the elements of the loop are specihed 
as projectors to a eigenstates of some Hamiltonian for N different settings of some 
parameters. 

As a generalization of the Berry phase, we ask about the relative phase around a 
loop on N projectors, each of which is Np dimensional. The physical motivation is 
that these are eigenspaces of a Hamiltonian, i.e., projectors to subspaces spanned by 
degenerate energy eigenstates iTSll . 


Discrete Wilson loop 

Consider a loop over N >3 sets of states from the Hilbert space, each consisting of 
Np S N orthonormal states, {\u„{k)) |n = 1,... ,Nf}, with k — 1,... ,N- We quan¬ 
tify the overlap between set k and set / by the Np x Np overlap matrix with 

elements 


= {u„{k) I 


(9.1) 


with n,m — I ... ,Nf- 

The discrete Wilson loop is the product of the overlap matrices along the loop, 

w = (9.2) 

We will be interested in the eigenvalues of the Wilson loop, with n — 1,... Nf. 

Wv,=XnYn, (9.3) 

where is the nth eigenvector, with n = 1,... ,Nf- 

Note that we could have started the Wilson loop, Eq. •EH. at the kth group 
instead of the hrst one. 




(9.4) 
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Although the elements of the matrix depend on the starting point k, the eigen¬ 
values A„ do not. Multiplying Eq. ( |9.3| l from the left by , we obtain 

^(k) (^^{k,k +^),, v^=K y„. (9.5) 


U{Nf) gauge invariance of the Wilson loop 

We will now show that the eigenvalues of the discrete Wilson loop over groups of 
Hilbert space vectors only depend on the linear spaces subtended by the vectors of 
each group. Each group can undergo an independent unitary operation to redefine 
the vectors, this cannot affect the eigenvalues of the Wilson loop. This is known as 
the invariance under a U[Nf) gauge transformation. 

A simple route to prove the U{Nf) gauge invariance is via the operator W defined 
by the Wilson loop matrix W in the basis of group 1, 


Np Nf 

E E (9.6) 

n=\m=\ 

This operator has the same eigenvalues as the Wilson loop matrix itself. It can be 
expressed using the projectors to the subspaces spanned by the groups of states, 

A = (9.7) 

n 

The Wilson loop operator reads, 

W =Pip2p^...pNpi. (9.8) 

We show this explicitly for A = 3, 


Np Np Np Np 

^ = E E E E l“«(i)) 

;7=1712 = 1723 = 1 m=l 

(m„(1) I m„,(2)) (m„,(2) I m„3(3)) (m„^(A) I (9.9) 

The generalization to arbitrary A > 3 is straightforward. 

Equation ( |9.81 l makes it explicit that the Wilson loop operator, and hence, the 
eigenvalues of the Wilson loop, are U{Nf) gauge invariant. 
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9.1.2 Wannier states for degenerate multiband one-dimensional 
insulators 

We now generalize the Wannier states of Sect. [13 to a one-dimensional insulator 
with Nf occupied bands. In case of nondegenerate bands, a simple way to go would 
be to debne a set of Wannier states for each band separately. However, time reversal 
symmetry forces degeneracies in the bands, at least at time reversal invariant mo¬ 
menta, and so this is not possible. Moreover, even in the nondegenerate case it could 
be advantegeous to mix states from different bands to create more tightly localized 
Wannier states. 

To be specibc, and to obtain efficient numerical protocols, we take a bnite sam¬ 
ple of a one-dimensional insulator of = 2M unit cells, with periodic boundary 
conditions. An orthonormal set of negative energy bulk eigenstates reads 

\%ik)) = \k) (g) \un{k)) = ^ e""*" \m) (g) \u„{k)) (9.10) 

VN 

with, as before, k G {5^., 25jt,... ,NS/;}, and 5^. = 2n/N. The index n labels the eigen¬ 
states, with n= 1,... for occupied, negative energy states. The \u„{k)) are the 
negative energy eigenstates of the bulk momentum-space Hamiltonian H{k). We 
will not be interested in the positive energy eigenstates. 

Although we took a specibc set of energy eigenstates above, because of degenera¬ 
cies at the time-reversal invariant momenta, we really only care about the projector 
P to the negative energy subspace. This is debned as 

Nf 

p=Lt = Ei^) f\^p{ky, (9.11) 

k n=l k 

Nf 

P(k)=£|M„(k))(M„(k)|. (9.12) 

n=\ 


Defining properties of Wannier states 

We will need a total number NpN of Wannier states to span the occupied subspace, 
\wn{j)), with j = 1,... ,N, and n = 1,... ,Nf- These are defined by the usual prop¬ 
erties: 
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) I = Sjijdn'n Orthoiiormal set (9.13a) 

N Nf 

52 52 |w„(J)) (w„(J)| = P Span occupied subspace (9.13b) 

j=\n=\ 

Vm : (m + 11 w„{i + 1)) = (m | w„{j)) Related by translation (9.13c) 

lim (w„(A^/2)| (i —A^/2)^ |w„(A^/2)) < oo Localization (9.13d) 

7V->oo 


with the addition in Eq. ( 9.13c| l defined modulo N. 

The Ansatz of Sect. 3.1 for the Wannier states, Eq. p.l l[ i, generalizes to the 
multiband case as 


1 Nf 

\Mj)) = ^ L £ Unp{k) I'Fpik)) . (9.14) 

nN P=1 

Thus, each Wannier state can contain contributions from all of the occupied bands, 
the corresponding weights given by a A:-dependent unitary matrix U{k). 


The projected unitary position operator 


As we did in Sect. 3.1 we will specify the set of Wannier states as the eigenstates 
of the unitary position operator restricted the occupied bands. 


Xp = 


(9.15) 


To obtain the Wannier states, we go through the same steps as in Sect. im with an 
extra index n. We outline the derivations and detail some of the steps below. You 
can then check whether the properties required of Wannier states, Eq. ( |9.13| ), are 
fulfilled, in the same way as in the single-band case. 

We note that for finite N, the projected unitary position Xp is not a normal oper¬ 
ator, i.e., it does not commute with its adjoint. As a result, its eigenstates form an 
orthonormal set only in the thermodynamic limit of N ^ Just as in the single¬ 
band case, this can be seen as a discretization error, which disappears in the limit 
A —)■ oo. 

The first step is to rewrite the operator Xp. Eor this, consider 

{%,{k')\X\%{k)) = 5,ps,x {un'ik+5k)\unik)) (9.16) 

where = 1 if kf — k + 5^, and 0 otherwise. Using this, the projected unitary 

position operator can be rewritten as 
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Nf 

k'k n',n=\ 

Nf 

= L L {Un'{k + 5k)\Un{k))-\%,{k + 5k)){^n{k)\. (9.17) 

k «',«=! 


Spectrum of the projected unitary position operator and the Wilson loop 

As in the single-band case, the next step is to consider Xp raised to the Mh power. 
This time, it will not be simply proportional to the projector P, however. Bearing in 
mind the orthonormality of the energy eigenstates, pPiik) \ ^„i{k')) = we 

find 

{Xpf = \'i'm{k)) {'¥n{k)\. (9.18) 

k mn 

The Wilson loop matrices as per Eq. ( |9.4| i, are all unitary equivalent, and have 
the same set of complex eigenvalues, 

= with n = l,...,NF, (9.19) 

|>1«|<1, enG[-7t,7t). (9.20) 

The spectrum of eigenvalues of Xp is therefore composed of the Mh roots of these 
eigenvalues, for j = and n = 1,... ,Nf, 

= e'e«/^+94+iog(|A„|)/iv^ ^ (A„j)^ = A„. (9.21) 


Wannier centers identified through the eigenvalues of the Wilson loop 


As in the single-band case. Sect. 3.1 we identify the phases of the eigenvalues Xnj 


of the projected position operator Xp with the centers of the Wannier states. There 
are Np sets of Wannier states, each set containing states that are spaced by distances 
of 1, 


{x)nj — 

{x)n = 


—ImlogAnj = {x)„+j-, 
On 


(9.22) 

(9.23) 


The phases of the Np eigenvalues of the Wilson loop W are thus identified with 
the Wannier centers, more precisely, with the amount by which the Np sets are 
displaced from the integer positions. 
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9.2 Time-reversal restrictions on Wannier centers 

We will now apply the prescription for Wannier states above to the one-dimensional 
insulators obtained as slices of a two-dimensional = — 1 time-reversal invariant 
insulator at constant ky. We will use the language of dimensional reduction, i.e., talk 
of the bulk Hamiltonian H{kx,ky) as describing an adiabatic particle pump with ky 
playing the role of time. We will use the Wannier center flow, i.e., the quantities 
{x)n = 9„{ky) /(In) to track the motion of the particles in the bulk during a fictitious 
pump cycle, ky = —n ^ n. 

Time-reversal symmetry places constrains the Wannier center flow in two ways: 
it enforces ky o —ky symmetry, i.e., 9„ {ky) — 9„> {—ky), and it ensures that for ky = 0 
and for ky = n, the are doubly degenerate. In this Section we see how these 
constraints arise. 


9.2.1 Eigenstates at k and —k are related 

A consequence of time-reversal symmetry is that energy eigenstates at k can be 
transformed to eigenstates at —k. One might think that because of time-reversal 
symmetry, energy eigenstates come in time-reversed pairs, i.e., that fH(—k)*f^ = 
H(k) would automatically ensure that |m„(— k)) = |M„(k)*). However, be¬ 

cause of possible degeneracies, this is not necessarily the case. The most we can say 
is that the state |m„(— k)) is some linear combination of time-reversed eigenstates, 

Nf Nf 

|M„(-k)) = x^ (B„^(k) |M„,(k)))* = |M„,(k)*). (9.24) 

m=\ m=\ 

The coefficients B„,„(k) define the unitary sewing matrix. An explicit formula for 
its matrix elements is obtained by multiplying the above equation from the left by 
{ua (k)* I V, with some a = 1,... ,Nf. This has the effect on the left- and right-hand 
side of Eq. ^2^ of 

{Ua{k)*\x^ \u„{-k)) = ((M„(-k)| f|M<,(k)*))*; (9.25) 

Nf 

(M«(k)*|r £ Bn^{k)*T\u,n{k)*)=Bn,{k)*, (9.26) 

m=l 

where for the last equation we used the unitarity of f and the orthonormality of the 
set \uynik)). Comparing the two lines above (and relabeling a m), we obtain 

B«m(k) = (M„(-k)|f|M„,(k)*). (9.27) 

Using this definition it is straightforward to show that the sewing matrix is unitary, 
and that B„,„(-k) = -B,„„(k). 
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9.2.2 Wilson loops at ky and —ky have the same eigenvalues 


To see the relation between the Wilson loops at ky and —ky, we first relate the pro¬ 
jectors to the occupied subspace at these momenta. We use a shorthand, 


p jP(^^ + j^k,ky) if;<0; 

\P{j5k,ky), ifj>0. 

Using Eq. (|9.24|i, and the unitarity of the sewing matrix B, we hnd 


(9.28) 


Nf 

P^j{-ky)=P{-k) = 52 |M„(-k)) (M„(-k)| 

n=l 

Np Nf Nf 

= L L L B,™(k)*f|«„,(k)*)B„^,(k)(M,„,(k)*|f+ 

n=lm=l m'=\ 

= Tpj{ky)*T:^ = Tpjikyfz^ 

The consequence of Eq. ( |9.29[ ) for the Wilson loop is 

W {-ky) = tW {kyf 

We write down the proof explicitly for N = 6, 

W{-ky) = P3{-ky)P2{-ky)Pl {-ky)Po{-ky)P^ I {-ky)P^2{-ky) P^ {-ky) 

= tP2{ky)^P_2{ky)^P-l{ky)'^P0{ky)^Pl {ky)'^P2{ky)'^P3{ky)'^f^ 

= TW{kyYt\ (9.31) 

the generalization to arbitrary even N follows the same lines. The set of eigenvalues 
of W is the same as that of its tranpose , as this holds for any matrix. Moreover, 
the unitary transformation of to does not change the eigenvalues either. 

To summarize, we find that the eigenvalues of the Wilson loop at —ky are the same 
as of the Wilson loop at ky, 


(9.29) 

(9.30) 


0n{ky) = 0n{-ky). (9.32) 

Erom Eq. ( |9.32[ ), we have that the Wannier center flow is symmetric around ky = 
0, and hence, also symmetric around ky = n. This means that it is enough to examine 
the Wannier centers from ky = 0 to ky = n. 
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9.2.3 Wilson loops at ky = 0 and ky = n are doubly degenerate 

We now concentrate on the two special values of the y wavenumber, ky =0 and 
ky = n, which are mapped unto themselves by time reversal. The one-dimensional 
Hamiltonians Hbuik(O) and Hbuik(?J^) time-reversal invariant, and thus, their every 
eigenstate must come with its time-reversed partner, obtained by the local time- 
reversal operation Due to Kramers theorem, an energy eigenstate and its time- 
reversed partner are orthogonal, (T'| =0. 

The Wilson loop W ?Aky = 0 and ky = 7Z is doubly degenerate. To show this, take 
an eigenstate of the Wilson loop, W\P) = X jf'). Using Eq. ( |9.31| ), we find 

X\^)=W I'E) = IT') = |T'); (9.33) 

wH|T'*) = A*f |T'*). (9.34) 

We obtained line ( |9.34| ) by multiplication from the left by f and complex conjuga¬ 
tion, and using the antisymmetry of f. In the final line, we have obtained that the 
Wilson loop W has a left eigenvector with eigenvalue X*. Since this is orthogonal to 
|T'), however, the right eigenvalue X must be at least twice degenerate. 


9.3 Two types of Wannier center flow 

We now examine the Wannier center flow, i.e., the functions 9„{ky), in time-reversal 
invariant two-dimensional insulators with = — 1. Due to the restrictions of 
ky o —ky symmetry and degeneracy at ky = 0, 7Z, we will And two classes of Wan¬ 
nier center flow. In the trivial class, the center flow can be adiabatically (i.e., contin¬ 
uously, while respecting the restrictions) deformed to the trivial case, with 9„ =0 
for every n and every ky. The topological class is the set of cases where this is not 
possible. 



Id 







Fig. 9.1 Wilson loop examples. 


To have a concrete example at hand, we examine the Wannier center flow for the 
BHZ model of the previous Chapter, Eq. (|8.38|), in a trivial (a) and in a topological 
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(b) case, and a third, more general topological (c) model. All three cases are covered 
by a modihed BHZ Hamiltonian, 

H{k) = So® [(M-|-cos^;f + cos^_v)^z + sin^_yd’y)] +Sj® sin^j;d;v-l-Sj;®C+ 

gSj 0 Gy (cos kjc + cos Iky — 2). (9.35) 

For case (a), we set g = 0, take the sublattice potential parameter m = 2.1, and cou¬ 
pling operator C = 0.026\.. This is adiabatically connected to the trivial limit of the 
BHZ model at m = +oo. In case (b), we set g — 0, take sublattice potential parameter 
M = 1, and coupling operator C = 0.36\, deep in the topological regime. In case (c) 
we add the extra term to the BHZ model to have a more generic case, with g = 0.1, 
and use m = 1, coupling C = O.ldy. The Wannier center flows for the three cases are 
shown in Fig. |9.1| 

We consider what adiabatic deformations of the Hamiltonian can do to the cen¬ 
ter flow. Focusing to ky = 0 n, the center flow consists of branches 9„{ky), that 
are continuous functions of ky, beginning at 0„(O) and ending at 0„(7r). Due to an 
adiabatic deformation, 

• A branch can bend while 0„(O) and 0„(7r) are fixed; 

• The endpoint at ky = 0 (or ky = ;r) of a branch can shift; in that case, the endpoint 
of the other branch, the Kramers partner at ky = 0 (or ky = n) is shifted with it; 

• Branches 0„ and 9m can recombine: a crossing between them at some ky can turn 
into an avoided crossing. 

Consider the example of Fig. |9.1| Bending of the branches and shifting of the 
endpoints can bring case (a) to a trivial case, where all branches are vertical, 
{ky) = 0 for every n and ky. Case (b) can be deformed to case (c). Notice, however, 
that neither cases (b) nor (c) can be deformed to the trivial case. 


9.3.1 Bulk topological invariant 

We define the bulk topological invariant A^buik^ by choosing some fixed 9 € [—tt, tt), 
and asking for the parity of the number of times the Wannier center flow crosses this 
9. In formulas. 


9&[-TZ,ny, (9.36) 

A), (0) = Number of solutions ky of 0„(ky) = 0; (9.37) 

A)iuik = ^52mod 2 (independent of 0). (9.38) 

The number Abuik is invariant under adiabatic deformations of the bulk Hamiltonian, 
as can be shown by considering the possible changes. Bending of a branch 0„ can 
create or destroy solutions of 9n{ky) — 9, but only pairwise. Shifting of the end- 
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point can create or destroy single solutions of 9n{ky) = 9, but in that case, a single 
solution of 9m{ky) = 9, is also created/destroyed, where 9,„ is the Kramers partner 
of 9n at the endpoint. Finally, recombination of branches cannot change the num¬ 
ber of crossings. The number N{9) is also invariant under a shift of 9, as already 
announced. A shifting of 9 is equivalent to a shifting of the Wannier center flow, 
whose effects we already considered above. 


The bulk topological invariant is the Z 2 invariant of the previous Chapter 

The full proof that the bulk invariant Abuik is the same as the parity D of the number 
of edge state pairs, Eq.( |8.40| l, is quite involved llJTl [161 . We content ourselves with 
just pointing out here that both Aj^uik and D represent obstructions to deform the 
Hamiltonian adiabatically to the so-called atomic limit, when the unit cells are com¬ 
pletely disconnected from each other. Clearly, switching of a charge pump requires 
that there are no edge states present, and therefore, Abuik = 0 requires D = 0. To 
show that the converse is true is more complicated, and we do not discuss it here. 


9.4 The Z 2 invariant for systems with inversion symmetry 


For two-dimensional time-reversal invariant insulators with inversion (i.e., parity) 
symmetry, the Z 2 topological invariant becomes very straightforward. We state the 
result below, and leave the proof as an exercise for the reader. 


Definition of inversion symmetry 


As introduced in Sect. 


3.2 


the operation of inversion, 77, acts on the bulk momentum- 


space Hamiltonian using an operator ft, by 


77H(k)77^^ =M(-k)^^ 


(9.39) 


We now require the operator ft not only to be independent of the wavenumber k, to 
be unitary, Hermitian, but also to commute with time reversal, i.e., 

77^77=1; 77^ = 1; (9.40) 


At a time-reversal invariant momentum, the Kramers pairs have the same 
inversion eigenvalue 

Consider the time-reversal invariant momenta (TRIM), Tj. In the BHZ model, these 
are {kx,ky) = (0,0), (0, k), {7Z,0), {n, k). In general there are such momenta in a 
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£/-dimensional lattice model. Each eigenstate |«(!] )) of the bulk momentum-space 
Hamiltonian at these momenta has an orthogonal Kramers pair 

Hirj)\uirj)} = E\uirj)) H{rj)^\u{rj)) = E^\u{rj)). (9.4i) 

If H is inversion symmetric, \u) can be chosen to be an eigenstate of ft as well, since 

7tH{rj)7t = H{-rj) = HiEj). (9.42) 


Therefore, 


7t\u{rj)) = ±\u{ri)). (9.43) 

The Kramers pair of |M(r))) has to have the same inversion eigenvalue as \u{rj)'), 

ft^\u{rj)) = ^7i\u{rj)) = ±^\u{rj)). (9.44) 

In a system with both time-reversal and inversion symmetry, we get 2^ topo¬ 
logical invariants of the bulk Hamiltonian, one for each time-reversal invariant mo¬ 
mentum Fj. These are the products of the parity eigenvalues of the occupied 

Kramers pairs at Fj. However, inversion symmetry is usually broken at the edges, 
and so these invariants do not give rise to robust edge states. 

The product of the inversion eigenvalues of all occupied Kramers pairs at all the 
time-reversal invariant momenta Fj is the same as the Z 2 invariant, 

( 9 . 45 ) 

i m 

We leave the proof of this useful result as an exercise for the reader. 


9.4.1 Example: the BHZ model 

A concrete example for inversion symmetry is given by the BHZ model of Sec. |8.2[ 
with no coupling C = 0. It can be checked directly that this has inversion symmetry, 
with 


BHZ; jt = SQ® dj. (9.46) 

To calculate the Z 2 invariant of the BHZ model, we take the four time-reversal 
invariant momenta, ki,k 2 ,k 3 ,k 4 , are the combinations of kx,ky with kjc = 0,7r and 
ky = 0,7r. The Hamiltonian EtsHzikxjky) at these momenta is proportional to the 
inversion operator. 
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^^B//z(ki =0,0) = (m + 2)^; //B//z(k 4 = 7r,7r) = (m —2)7r; (9.47) 

Hbhz{^ 2 =0,7t) = uft; HBHzi^s = 7t,0) = uft. (9.48) 

In these cases the Hamiltonian and the inversion operator obviously have the same 
eigenstates. At each TRIM, two of these states form one occupied Kramers pair and 
the two others one empty Kramers pair. If m > 2, at all four TRIM, the occupied 
Kramers pair is the one with inversion eigenvalue (parity) of —1, and so Eq. ( |9.45| l 
gives A^buik = 0. Likewise, if m < —2, the eigenvalues are all +1, and we again obtain 
Abuik = 0. For 0 < M < 2, we have P eigenvalues —1,—1, + 1, — 1 at the four TRIM 
ki,k 2 ,k 3 ,k 4 , respectively, whereas if —2 < m < 0, we have —1,+1,+1,+l.In both 
cases, Eq. ( |9.45| ) gives A^buik = 1- This indeed is the correct result, that we obtained 
via the Chern number earlier. 


Problem s 


9.1. Inversion symmetry and interlayer coupling in the BHZ model 

Consider the BHZ model with layer coupling C = Cf, C) Gy. This breaks the inver¬ 
sion symmetry n = G^. Nevertheless, the Wannier centers of the Kramers pairs 

at ky — 0 and ky — n are stuck to 0 = 0 or 0 = tt, and are only shifted by the extra 
term s^Gy added to the BHZ model in Eq. ( |9.35| ). Can you explain why? (hint; 
extra inversion symmetry) 


9.2. Proof of the formula for the Z 2 invariant of an inversion-symmetric topo¬ 
logical insulator 


Show, using the results of Sect. 3.2 that the Z 2 invariant of a two-dimensional 


time-reversal invariant and inversion symmetric insulator can be expressed using 

Eq. ( |9A31 i. 












Chapter 10 

Electrical conduction of edge states 


It is well known that the electrical conduction of ordinary metallic samples at room 
temperature shows the following two characteristics. First, there is a linear relation 
between the electric current I that flows through the sample and the voltage V that 
drops between the two ends of the sample; I/V = G = R \ where G (R) is the 
conductance (resistance) of the sample. Second, the conductances G,- of different 
samples made of the same metal but with different geometries show the regularity 
GiLijAi — a for V/, where L,- is the length of the sample and A,- is the area of its 
cross section. The material-specific quantity CJ is called the conductivity. Conduc¬ 
tors obeying both of these relations are referred to as Ohmic. 

Microscopic theories describing the above behavior (e.g., Drude model, Boltz¬ 
mann equation) rely on models involving impurities, lattice vibrations, and electron 
scattering within the material. Electrical conduction in clean (impurity-free) nanos¬ 
tructures at low temperature might therefore qualitatively deviate from the Ohmic 
case. Here, we demonstrate such deviations on a simple zero-temperature model of 
a two-dimensional, perfectly clean, constant-cross-section metallic wire, depicted 
in Fig. 10.1^ . Then we describe how scattering at static impurities affects the con¬ 
duction in general. Finally, as a central result in the field of topological insulators, 
we point out that electrical conduction via the edge states of two-dimensional topo¬ 
logical insulators shows a strong robustness against such impurities. 


10.1 Electrical conduction in a clean quantum wire 

As shown in Fig. |10.1^ , take a wire that lies along the x axis and has a finite width 
in the y direction. The latter might be defined by an electric confinement potential 
or the termination of the crystal lattice. Each electronic energy eigenfunction \l,k) 
in such a wire is a product of a standing wave along y, labeled by a positive integer 
/, and a plane wave propagating along x, labeled by a real wave number k (see Eq. 
( |8.41| l). A typical set of dispersion relations F/i- (‘subbands’) for three different / 
indices is shown in Fig.llO.lf). 
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Fig. 10.1 (a) Schematic representation of a clean quantum wire contacted to two electron reser¬ 
voirs (contacts), (b) Occupations of electronic states in the contacts and the quantum wire in the 
nonequilibrium situation when a finite voltage V is applied between the left and right reservoirs. 


We also make assumptions on the two metallic contacts that serve as source and 
drain of electrons. We assume that the electrons in each contact are in thermal equi¬ 
librium, but the Fermi energies in the contacts differ by = |e|y > 0. (Note 

that in this chapter, proper physical units are used, hence constants such as the el¬ 
ementary charge \e\, reduced Planck’s constant h, lattice constant a are reinstated.) 
We consider the linear conductance, that is, the case of an infinitesimal voltage 
|e|y 0. We further assume that both contacts absorb every incident electron with 
unit probability, and that the energy distribution of the electrons they emit is the 
thermal distribution with the respective Fermi energy. 

These assumptions guarantee that the right-moving (left-moving) electronic states 
in the wire are occupied according to the thermal distribution of the left (right) con¬ 
tact, as illustrated in Fig. |10.1| r. Now, we work with electron states normalized to 
the area of the channel. It is a simple fact that with this normalization convention, 
a single occupied state in the /th channel, with wave number k carries an electric 
current of where Nxa is the length of the wire, and v/<. = ^ is the group 

velocity of the considered state. Therefore, the current flowing through the wire is 
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^ = - k| - I^l) -f{Eik - IJ-r)] , (10.1) 

where /(e) = ^exp + 1 j is the Feimi-Dirac distribution. Converting the k 

sum to an integral via ^ g ... yields 

rn/a flh 1 JIT,, 

^=-klE / :^.^[f{Eik-liL)-f{E,t-liR)]- (10.2) 

J-%la2nn dk 


The Fermi-Dirac distribution has a sharp edge at zero temperature, implying 

Icl Isl 6^ 

/=-Um/ dE = -'-^{^L-^R)M^M-V (10.3) 

h J jXf^ h h 

Note that the first equality in ( |10.3| l holds only if the number of subbands intersected 
by jLli and ^Lr are the same, which is indeed the case if the voltage V is small enough. 
The number of these subbands, also called ‘open channels’, is denoted by the integer 
M. From ( |10.3| l it follows that the conductance of the wire is an integer multiple of 
/h (commonly referred to as ‘quantized conductance’): 

G=—M. (10.4) 

h 


The numerical value of /h is approximately 40 /rS (microsiemens), which corre¬ 
sponds to a resistance of approximately 26 ki2. Note the the conductance quantum 
is defined as Go = 2e^/h, i.e., as the conductance of a single open channel with 
twofold spin degeneracy. 

It is instructive to compare the conduction in our clean quantum wire to the ordi¬ 
nary Ohmic conduction summarized above. According to ( |10.3| l, the proportionality 
between voltage and current holds for a clean quantum wire as well as for an ordi¬ 
nary metal. However, the dependence of the conductance on the length of the sample 
differs qualitatively in the two cases: in an ordinary metal, a twofold increase in the 
length of the wire halves the conductance, whereas the conductance of a clean quan¬ 
tum wire is insensitive to length variations. 

Whether the conductance of the clean quantum wire is sensitive to variations 
of the wire width depends on the nature of the transversal modes. Conventional 
quantum wires that are created by a transverse confinement potential have a subband 
dispersion similar to that in Fig. 10.1 1 . There, the energy separation between the 
subbands decreases as the width of the wire is increased, therefore the number of 
subbands available for conduction increases. This leads to an increased conductance 
for an increased width, similarly to the case of ordinary metals. If, however, we 
consider a topological insulator, where the current is carried by states localized to 
the edges of the wire, the conductance of the wire will be insensitive to the width of 
the wire. 
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10.2 Phase-coherent electrical conduction in the presence of 
scatterers 


Having calculated the conductance ( |10.4| l of a clean quantum wire, we now describe 
how this conductance is changed by the presence of impurities. We analyze the 
model shown in Fig. |10.2| where the disordered region, described by a scattering 
matrix S, is connected to the two contacts by two identical clean quantum wires, 
also called ‘leads’ in this context. 



Fig. 10.2 Simple model of a phase-coherent conductor in the presence of scatterers. The ideal 
contacts L and R are connected via ideal leads to the disordered region represented by the scattering 
matrix S. 


First, we consider the case when each lead supports a single open channel. The 
current in the lead connecting contact L and the scattering region consists of a contri¬ 
bution from right-moving states arriving from contact L and partially backscattered 
with probability R — |rp, and from left-moving states arriving from contact R and 
partially transmitted with probability T' = 

I = -\e\^Y.^k[{l-R{Ek))fL{Ek)-T'{Ek)fR{Ek)\ (10.5) 

Converting the k sum to an integral, assuming that the transmission and reflection 
probabilities are independent of energy in the small energy window between jiR and 
jJLi, and using l—R = T = T', we arrive at 

\e\ fp' 

/=-yr/ dE\fR(E)-fR(E)\ = -TY, (10.6) 

n J fii; n 

which implies that the conductance can be expressed through the transmission co¬ 
efficient T: 
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G=^—T. (10.7) 

h 

The result ( |10.7| l can be straightforwardly generalized to the case when the leads 
support more than one open channel. The generalized result for the conductance, 
also known as the Landauer formula, reads; 

2 M 

G=-rY.Tn, ( 10 . 8 ) 

^ n=\ 

where T„ are the transmission eigenvalues of the scattering matrix i.e., the real eigen¬ 
values of the Hermitian matrix tt\ as defined in the preceding chapter. 


10.3 Electrical conduction in two-dimensional topological 
insulators 

After presenting the Landauer formula as a generic tool to describe electrical con¬ 
duction of a phase-coherent metal, we will use it know to characterize the conduc¬ 
tances of various two-dimensional topological insulator samples. 


10.3.1 Chern Insulators 


In Sect. 6.2 we have seen that an impurity-free straight strip of a topologically non¬ 


trivial Chern Insulator supports edge states. The relation between the Chem num¬ 
ber Q of the Chem Insulator and the numbers of edge states at a single edge at a 
given energy E, propagating ‘clockwise’ {N+{E)) and ‘anticlockwise’ (N-{E)), is 
Q — Alf (£) — N-{E). In addition, in Sect. 6.3 it was shown that any segment of 
the edge of a disordered Chern Insulator with Chern number Q and an arbitrary ge¬ 
ometry supports \Q\ chiral edge modes. Here we show that existence of these edge 
modes leads to experimentally detectable effects in the electrical transport through 
Chern Insulator samples. 

We consider a transport setup where the Chern Insulator is contacted with two 


metallic electrodes, as shown in Fig. 10.3 In this discussion, we rely on the usual 


assumptions behind the Landauer formula; phase-coherence of the electrons, good 
contact between contacts and sample, and large spatial separation of the two elec¬ 
trodes ensuring the absence of tunneling contributions to the conductance. 

In the following list, we summarize how the phase-coherent electrical conduc¬ 
tance of a Chern Insulator varies with the sample geometry, absence or presence of 
disorder, and the value of the electronic Fermi energy. 


1. Disorder-free sample with a strip geometry (see Fig. 10.1 1 ) 
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Fig. 10.3 A disordered sam¬ 
ple of Chem insulator, with 
contacts 1 and 2, that can be 
used to pass current through 
the sample in order to detect 
edge states. 



a. Fermi energy lies in a band. In this case, the sample is a clean quantum wire 
(see Sect. |10.1| l with an integer number of open channels. The correspond¬ 
ing transversal wave functions might or might not be localized to the sample 
edges, and therefore the number of channels might be different from any com¬ 
bination of Q, N+ or N-. According to Eq. 10.4[ the conductance of such a 
clean quantum wire is quantized and insensitive to the length of the sample. 
Furthermore, the conductance grows in a step-like fashion if the width of the 
sample is increased. 

b. Fermi energy lies in the gap. The sample is a clean quantum wire with open 
channels that are all localized to the sample edges. The number of those chan¬ 
nels is A^+(£’) +A^- (£), where E is the Fermi energy; in each of the two pos¬ 
sible direction of current flow, there are N+ channels on one edge and N- on 
the other edge that contribute to conduction. Conductance is finite and quan¬ 
tized, a behavior rather unexpected from an insulator. The conduction is not 
Ohmic, as the conductance is insensitive to both the length and the width of 
the sample. 


2. Disordered sample with an irregular shape (see Fig. 10.3 i; 
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a. Fermi energy lies in a band. Because of the presence of disorder, the electrical 
conduction of such a sample might be Ohmic. There are no protected edge 
states at the Fermi energy. 


b. Fermi energy lies well within the gap. According to Sect. 6.3 any edge seg¬ 
ment of such a sample supports Q reflectionless chiral edge modes at the 
Fermi energy. Therefore, conductance is typically quantized, G = \Q\e'^/h, 
although, disorder permitting, it might in principle be larger than this value. 
The quantized conductance is insensitive to changes in the geometry or the 
disorder configuration. This transport property, unexpected for an insulator, 
let alone for one with disorder, is a hallmark of Chern Insulators. 

In the case of two-dimensional samples there is often an experimental possibility of 
tuning the electronic Fermi energy in situ by controlling the voltage applied between 


the sample and a nearby metallic plate (gate electrode), as discussed in Sect. 10.4 


This allows, in principle, to observe the changes in the electrical conduction of the 
sample as the Fermi energy is tuned across the gap. 


10.3.2 Two-dimensional time-reversal invariant topological 
insulators with — 1 

In the following list, we summarize the predictions of the Landauer formalism for 
the conductance of two-dimensional time-reversal invariant topological insulators 
with = — 1 (‘D = 1 insulators’ for short). 

1. Disorder-free sample with a strip geometry: 

a. Fermi energy lies in a band. A simple consequence of the Landauer formula is 
that phase-coherent conductance of an impurity-free D = 1 topological insu¬ 
lator of the strip geometry shown in Fig. |10.1^ is quantized. The conductance 
grows if the width of the strip is increased, but insensitive to change in the 
length. 

b. Fermi energy lies in the gap. Only edge channels are open in this case. These 
also provide conductance quantization. As the number of edge-state Kramers 
pair per edge is odd, the conductance might be 2e^/h, 6e^ jh, lOe^/Zz, etc. 
Conductance is insensitive to width or length changes of the sample. 

2. Disordered sample with an irregular shape and TRS disorder: 

a. Fermi energy lies in a band. The electrical conduction might be Ohmic. 

b. Fermi energy lies in the gap. We have shown in Chapt.j^that a D = 1 insu¬ 
lator supports one protected edge-state Kramers pair per edge, which allows 
for reflectionless electronic transmission if only TRS defects are present. The 
Landauer formula ( |10.8| ) implies, for typical cases, G = 2e^/h for such a sam¬ 
ple, as one edge state per edge contributes to conduction. The conductance 
might also be larger, provided that the number of edge-state Kramers pairs is 
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larger than 1 and disorder is ineffective in reducing the transmission of the 
topologically unprotected pairs. 

We note that in real materials with Z 2 invariant D—\, various mechanisms might 
lead to backscattering and, in turn, to G < 2e^/h. Examples include TRS-breaking 
impurities, TRS impurities that bridge the spatial distance between the edges (see 
Chapt. [^, hybridization of edge states from opposite edges in narrow samples, and 
inelastic scattering on phonons or spinful impurities. 


10.4 An experiment with HgTe quantum wells 

Electrical transport measurements IZD on appropriately designed layers of the semi¬ 
conductor material mercury-telluride (HgTe) show signatures of edge-state conduc¬ 
tion in the absence of magnetic held. These measurements are in line with the the¬ 
oretical prediction that a HgTe layer with a carefully chosen thickness can realize 
a topologically nontrivial (D = 1) two-dimensional time-reversal invariant insulator 
with = — 1. In this section, we outline the main hndings of this experiment, as 
well as its relation to the BHZ model introduced and discussed in chapter]^ 



X 


Fig. 10.4 Schematic representation of a HgTe quantum well of width d, sandwiched between two 
Hg,jCdi_^Te layers. Electrons are confined to the HgTe layer, and their Eermi energy can be tuned 
in situ by adjusting the voltage Vgate of the metallic electrode on the top of the sample (black). Eor 
a more accurate description of the experimental arrangement, see do). 


The experiments are performed on sandwich-like structures formed by a few- 
nanometer thick HgTe layer (quantum well) embedded between two similar layers 
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of the alloy Hg^Cdi_jcTe, as shown in Fig. |10.4| (In the experiment reported in ETIl . 
the alloy composition x = 0.3 was used.) In this structure, the electronic states with 
energies close the Fermi energy are conhned to the HgTe layer that is parallel to the 
x-y plane in Fig. 10.4 The energy corresponding to the conhnement direction z is 
quantized. The carriers are free to move along the HgTe layer, i.e., parallel to the x-y 
plane, therefore two-dimensional subbands are formed in the HgTe quantum well. 
Detailed band-structure calculations of Ref. 0 show that as the the thickness d of 
the HgTe layer is decreased, the lowermost conduction subband and the uppermost 
valence subband touch at a critical thickness d — dc, and the gap is reopened for even 
thinner HgTe layers. (For the alloy composition .r = 0.3 used in the experiment, 
the critical thickness is dc ~ 6.35 nm.) This behavior is illustrated schematically 
in Fig. |10.5| which illustrates the electronic dispersions of the uppermost valence 
subband and the the lowermost conduction subband, at the center of the BZ, for 
three different thicknesses of the quantum well. 



Fig. 10.5 Evolution of the two-dimensional band structure of a HgTe quantum well as a function 
of its thickness d. (a) For a thin quantum well below the critical thickness d <dc, the band structure 
has a gap and the system is a trivial insulator, (b) At a critical thickness d = dc, the band gap closes 
and the system is metallic, (c) For a thick quantum well with d > d^, the band gap reopens and the 
system becomes a two-dimensional topological insulator. 


Band-Structure calculations have also revealed a connection between the sub¬ 


bands depicted in Fig. 10.5 and the BHZ model introduced and discussed in 
Chapt.|gThe4 X 4 effective Hamiltonian describing the two spinful two-dimensional 
subbands around their extremum point at the centre of the HgTe Brillouin zone re¬ 
sembles the low-energy continuum Hamiltonian derived from the BHZ lattice model 
in the vicinity of the u « —2 value. Changing the thickness d of the HgTe layer cor¬ 
responds to a change in the parameter u of the BHZ model, and the critical thickness 
d — dc corresponds to m = —2 and, consequently, a zero mass in the corresponding 
two-dimensional Dirac equation. 

As a consequence of the strong analogy of the band structure of the HgTe quan¬ 
tum well and that of the BHZ model, it is expected that either for d < dc or for 
d>dc the material is a D = 1 insulator with a single Kramers pair of edge states. 
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10 Electrical conduction of edge states 


Arguments presented in 13 suggest that the thick quantum wells with d > dc are 
topologically nontrivial. 

Electrical transport measurements were carried out in HgTe quantum wells pat¬ 
terned in the Hall bar geometry shown in Fig. 10.6 The quantity that has been 


used in this experiment to reveal edge-state transport is the four-terminal resistance 
^ 14,23 = k 23 //i 4 , where V 23 is the voltage between contacts 2 and 3, and I 14 is the 
current flowing between contacts 1 and 4. This quantity was measured for 

various devices with different thicknesses d, below and above the critical thickness 
dc, of the HgTe layer, and for different values of the Fermi energy. The latter can be 
tuned in situ by controlling the voltage between the HgTe layer and a metallic ‘gate’ 
electrode on the top of the layered semiconductor structure, as shown in Fig. |10.4| 



w 


Fig. 10.6 HgTe quantum well patterned in the Hall-bar geometry (gray area). Numbered terminals 
lead to metallic contacts. Solid and dashed lines depict counterpropagating edge states. 


To appreciate the experimental result, let us first derive the four-terminal resis¬ 
tance Ru ,23 for such a device. To this end, we express V 23 with / 14 . Ohm’s law 
implies V 23 = / 23 /G 23 , where 733 is the current flowing through the edge segment 
between contacts 2 and 3, whereas G 23 is the conductance of that edge segment. 
Furthermore, as the current /14 flowing through terminals 1 and 4 is equally divided 
between the upper and lower edges, the relation 123 = / 14/2 holds, implying the 
result 7 ?i4,23 = 1/(2G23). 

If the Fermi energy lies in the bands neighboring the gap, then irrespective of 
the topological invariant of the system, the HgTe quantum well behaves as a good 
conductor with G 23 ^ e^/h, implying 7?i4,23 ^ h/e^. If the Fermi energy is tuned 
to the gap in the topologically nontrivial case d>dc, then G 23 = e^/h and therefore 
f^i 4.23 = h/{2e^). This holds, of course, only at a temperature low enough and a 
sample size small enough such that phase coherence is guaranteed. The presence of 
static time-reversal invariant defects is included. If the system is topologically trivial 
id < dc), then there is no edge transport, and the quantum well is a good insulator 
with 7? 14,23 ^ h/e^. 

The findings of the experiments are consistent with the above expectations. Fur¬ 
thermore, the four-terminal resistance of topologically nontrivial HgTe layers with 
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different widths were measured, with the resistance found to be an approximately 
constant function of the width W of the Hall bar. This is a further indication that the 
current in these samples is carried by edge states. 

To wrap up this chapter, we note that InAs/GaSb bilayer quantum wells are an 
alternative semiconductor material system where two-dimensional topological insu¬ 
lators can be realized ll22l 1^. Graphene is believed to be a two-dimensional topo¬ 
logical insulator as well CD, even though its energy gap between the valence and 
conduction band, induced by spin-orbit interaction and estimated to be of the or¬ 
der of jJ.eV, seems to be too small to allow for the detection of edge-state transport 
even at the lowest available temperatures. The concept of a time-reversal invariant 
topological insulator can be extended to 3D crystals as well, where the role of the 
edge states is played by states localized to the two-dimensional surface of the 3D 
material. The description of such systems is out of the scope of the present course; 
the interested reader might consult, e.g.. Refs. unmi. 
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